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An  experimental  and  analytical  program  was  undertaken  to  lnveati- 
gate  the  behavior  of  frost  formation  under  an  electric  wind 

induced  by  a  positive  wire-plane  corona  in  a  non-uniform  electric  field. 
A  vertical  test  plate  precooled  to  subfreezing  temperatures  was  exposed 
to  the  ati> lent  air  and  the  effects  of  electric  fields  were  studied 
through  the  application  of  a  high  voltage  between  the  grounded  test 
plate  and  a  corona  wire.  Considerable  Increases  in  heat  and  mass  trans¬ 
fer  were  found.  The  changes  were  attributed  to  the  alteration  of  flow 
field  from  free  convection  to  forced  convection. 

A  simple  theoretical  model  was  developed  for  predicting  the  effects 
of  electric  fields  upon  the  heat  and  mass  transfer  processes  associated 
with  frost  formation.  Two  regions  of  the  flow  field  created  by  the 
impinging  electric  wind  are  of  primary  importance  in  beat  and  mass  trans¬ 
fer,  namely  the  s taxation  flow  and  the  wall  Jet  regimes.  In  the  stag" 
nation  flow,  a  linear  relation  for  the  local  heat  and  mass  transfer 
coefficients  was  assisaed  between  the  stagnation  point  and  the  wall  .let 
region.  In  the  wall  Jet,  the  analysis  was  based  on  laminar  flow.  The 
mass  transfer  predicted  by  this  model  agrees  quite  satisfactorily  with 
experimental  results,  but  discrepancy  between  theory  and  experiment 
was  observed  In  heat  transfer. 

An  analysis  was  also  made  of  a  simplified  model  of  transient  heat 
conduction  in  a  frost  layer.  Using  a  method  due  to  Bortnov  the  position 
and  tbs  temps rr';ure  of  the  progressing  frost-edr  interface  was  obtained 
in  a  reries  expansion  in  powers  of  ft.  The  coefficients  up  to  the 
power  n  -  6  were  given. 
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CHAPTER  I 


INTRODUCTION 


1.1  INTRODUCTORY  REMARKS 

In  recent  years,  there  has  been  a  significant  increase  of  interest 
in  electrohydrodynamics .  Problems  in  this  area  are  confronted  in  deve¬ 
loping  controlled  thermal  nuclear  devices,  ion  drag  pumping  arrangements, 
ion  propulsion  systems  as  veil  as  exotic  heat  transfer  schemes,  to 
mention  four  areas . 

In  the  area  of  heat  transfer,  investigations  have  thus  far  mostly 
been  confined  to  the  case  of  two-phase  (liquid-gas)  heat  transfer. 

This  type  of  heat  transfer  is  highly  dependent  on  the  force  field 
experienced  by  the  fluid.  It  has  been  shown  by  Choi1  that  electric 
field  can  enhance  significantly  the  boiling  heat  transfer  properties 
of  a  device.  In  their  recent  work,  Velkoff  and  Miller2,  and  Choi3 
demonstrated  independently  a  similar  improvement  in  condensation  heat 
transfer.  The  instability  of  the  liquid-gas  interface  caused  by  polari¬ 
zation  and  free  charge  interaction  under  the  action  of  an  external 
electric  field  has  direct  bearing  on  these  phenomena. 

The  effects  of  electric  wind  on  heat  transfer  frcm  a  heated  plate 
to  ambient  air  has  been  investigated  by  Velkoff4.  He  showed  that  con¬ 
siderable  increase  in  heat  transfer  rate  could  be  obtained  through  the 
use  of  the  corona  discharge.  Similar  observations  have  been  reported 
by  Godfrey5  for  forced  convection  at  low  free  stream  velocity. 

The  influence  of  electric  fields  on  heat  and  mass  transfer  pro¬ 
cesses  between  humid  air  and  solid  surface  under  frosting  conditions 
has  received  little  attention.  The  problem  is  complicated  because  of 
the  following  facts.  First,  frost  formation  is  essentially  a  transient 
process.  As  the  frost  grows  not  only  the  heat  transfer  rate  and  the 
frost-air  interface  temperature  vary  continuously  with  time,  but  also 
the  frost  properties  change  unpredictably  frcm  one  instant  to  another. 
Second,  the  frost  layer  is  a  heterogeneous  porous  medium  whose  struc¬ 
ture  strongly  depends  on  the  conditions  under  which  it  is  formed.  It 
is  impossible  at  the  present  time  to  delineate  with  any  degree  of 
assurance  the  quantitative  effect  of  the  system  independent  variables 
upon  the  density  of  deposited  frost.  Correlation  of  the  experimental 
frost  thermal  conductivity  data  with  the  analytical  results  evaluated 
on  the  basis  of  the  available  methods  for  predicting  the  conductivities 
of  heterogeneous  materials  is,  as  yet,  not  satisfactory.  Finally,  it 
does  not  appear  that  the  microscopic  study  of  the  physical  process  of 
frosting  has  received  attention  previously.  In  the  absence  of  such 
fundamental  knowledge,  difficulties  are  encountered  in  predicting  the 
possible  effects  of  electric  fields  on  the  formation  of  frost.  Further¬ 
more,  unlike  the  cases  of  condensation  and  boiling  heat  transfer  where 
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the  electric  field  effect  manifests  itself  in  the  form  of  the  instabil¬ 
ity  of  the  deformable  liquid-gas  interface,  the  frost-air  interface 
can  hardly  he  expected  to  play  a  similar  role,  for  such  interface  is 
non-deformable . 

In  spite  of  the  absence  in  the  literature  of  studies  on  frost 
formation  in  an  electric  field,  some  work  on  the  effects  of  electric 
fields  upon  the  growth  of  ice  crystals  from  the  vapor  has  been  reported 
by  meteorologists  in  the  past  few  years.  In  a  study  of  the  effects  of 
electric  fields  on  ice  crystals  growing  from  the  vapor,  Bartlett  et  al6 
observed  that  ice  crystals  grew  rapidly  in  the  form  of  long  thin  needles 
when  an  electric  field  in  excess  of  a  certain  minimum  value  was  applied. 
Their  experiments  were  carried  out  in  a  diffusion  chamber  under  super¬ 
saturation  conditions.  Whether  similar  phenomenon  will  be  observed 
when  frost  layer  grows  from  the  ambient  air  is  still  a  question  mark. 

In  his  study  of  the  electric  deposition  of  ice  in  a  non-uniform 
field,  Arabadzhi7  mounted  two  point  electrodes  right  above  a  container 
with  warm  water  which  released  the  vapor  into  the  surrouniing  air.  A 
positive  potential  of  9  kv  was  applied  to  one  electrode  and  the  other 
was  grounded.  When  the  air  temperature  was  between  -10 °C  to  -12°C,  he 
observed  that  on  the  zero  electrode  first  clear  ice  was  deposited  and 
then  hoar  frost  began  to  form.  When  the  temperature  was  below  -15  °C , 
an  abundance  of  fluffy  frost  deposited  right  away  at  the  zero  electrode. 
He  beleived  that  the  electro-deposition  effects  were  due  to  corona  dis¬ 
charges  between  the  electrodes.  However,  no  detailed  description  of 
this  mechanism  was  given. 

It  is  the  purpose  of  this  study  to  undertake  an  exploratory  inves¬ 
tigation  of  the  influence  of  electric  fields  on  heat  and  mass  transfer 
processes  in  association  with  the  formation  of  frost.  In  the  next 
section,  sane  preliminary  experimental  observations  are  reported.  And 
in  the  section  following,  the  problem  to  be  treated  is  described. 


1.2  PRELIMINARY  EXPERIMENTAL  OBSERVATIONS 

To  determine  the  possible  effects  of  electric  fields  on  the  frost 
formation,  preliminary  exploratory  tests  were  conducted  using  a  variety 
of  electrode  configurations  to  provide  the  fields.  In  all  the  tests, 
the  test  plate  which  served  as  the  forsting  surface  was  grounded,  unless 
otherwise  stated.  The  temperature  of  the  plate  was  maintained  at  about 
10 °F.  The  average  temperature  and  specific  humidity  of  the  ambient 
air  were  about  73 °F  and  0.01  lb  of  vapor  per  lb  of  dry  air  respectively. 
The  experimental  setup  is  to  be  described  in  Chapter  VI.  Here  only 
the  following  observations  are  reported: 

1.  Uniform  field.  A  4-mesh  screen  was  used  as  positive 
electrode  to  provide  uniform  electric  field.  The  screen  was 
parallel  to  and  at  a  distance  £  inch  from  the  test  plate. 
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When  the  potential  applied  to  the  screen  was  increased  to 
12  lev,  no  observable  changes  in  the  frosting  phenomenon  were 
observed.  In  this  case,  no  current  was  measured  across  the 
field  since  there  was  no  corona  discharge  and  the  electric 
wind  was  absent. 

2.  Positive  wire-plane  corona.  A  non-uniform  electric 
field  was  produced  by  applying  a  high  potential  to  a  0.00 6 
inch  corona  wire  also  lacated  at  £  inch  from  the  plate.  When 
the  potential  was  below  the  critical  potential  for  the  corona 
discharge,  approximately  7*6  kv  for  this  case,  there  was  no 
current  and  no  changes  in  frost  formation  could  be  observed. 
However,  when  the  potential  was  increased  above  the  critical 
potential,  the  frost  first  appeared  directly  under  the  wire 
in  the  form  of  a  fine  strip,  then  grew  laterally  to  the  sides 
of  the  wore,  forming  a  wider  frost  strip.  Thus  at  the  initial 
stages  of  frost  formation  the  frost  layer  thickness  decreased 
with  increasing  distance  from  the  centerline  of  the  strip. 

As  the  process  of  frosting  continued,  the  thickness  of  the 
layer  was  more  or  less  evened  out  for  the  whole  plate.  But 
a  definite  non-uniformity  in  the  density  of  the  frost  layer 
was  evidenced  by  the  reduction  in  color  of  the  frost  layer  as 
the  distance  from  the  frost  strip  centerline  was  increased. 

Upon  close  examination  of  the  frost  strip,  it  was  further 
noticed  that  the  frost  thickness  was  not  uniform  along  the 
strip.  In  certain  portions  of  the  strip,  the  deposited  frost 
had  the  appearance  of  periodic  ridges  transverse  to  the  strip. 
These  could  best  be  observed  when  the  current  was  around  150  ua. 
The  scale  of  periodicity  of  the  ridges  is  about  l/8  inch. 
Photographs  of  the  growth  of  frost  and  distribution  of  frost 
thickness  along  the  test  plate  are  shown  in  Appendix  A. 

3.  Positive  wire-screen  corona.  In  this  test,  grounded 
4-mesh  screen  was  inserted  between  the  wire  and  the  frosting 
surface,  the  latter  being  electrically  neutral.  When  the 
corona  discharge  took  place  between  the  wire  and  the  screen 
the  induced  electric  wind  passed  through  the  screen  impinging 
normally  on  the  frosting  surface.  A  frost  strip  appeared  as 
in  the  case  of  positive  wire-plane  corona,  but  transverse 
ridges  could  no  longer  be  observed.  Instead,  transverse 
grooves  appeared  under  each  screen  wire.  Apparently,  the 
grooves  were  due  to  the  obstruction  to  the  passage  of  the  wind 
by  the  screen  wires. 

4.  Negative  and  a-c  wire -plane  coronas.  When  the  wire 
was  subjected  to  a  negative  or  a-c  potential  obove  the  criti¬ 
cal  potential,  the  frost  formed  along  under  the  wire  com¬ 
prised  of  irregular  patches.  This  may  be  accounted  for  by 
the  fact  that  corona  on  negative  or  a-c  wire  is  concentrated 
at  points  along  the  wire.  Observations  indicated  that  for  a 
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given  current,  more  frost  was  depoisted  for  positive  than 
for  negative  or  a-c  corona. 

5.  Positive  point-plane  corona.  When  positive  poten¬ 
tial  higher  than  the  critical  potential  was  applied  to  a 
needle  point,  the  frost  layer  was  approximately  synaetrical 
about  the  point,  as  might  be  expected.  Because  of  the  low 
current  obtainable  for  a  given  applied  potential,  its 
influence  on  frost  deposition  is  not  as  significant  as 
positive  wire  corona. 

On  the  basis  of  the  above  experimental  findings,  the  following 
comments  can  be  made: 

1.  A  uniform  electric  field  does  not  seem  to  have  any 
observable  effects  on  the  frosting  phenomenon,  at  least  at 
the  efirly  stages  of  frost  formation.  Because  of  the  limita¬ 
tions  of  the  experimental  setup,  no  tests  which  last  longer 
than  one  hour  have  been  attempted.  It  remains  to  be  deter¬ 
mined  if  a  unifom  field  will  affect  the  growth  of  frost  at 
a  later  stage  when  loose  frost  begins  to  form. 

2.  A  non-uniform  field  has  significant  effects  on  frost 
formation  only  when  corona  discharges  take  place  in  the  space 
between  the  electrodes.  Of  the  various  corona  discharges  and 
in  so  far  as  the  enhancement  in  the  amount  of  deposited  frost 
is  concerned,  the  positive  wire-plane  corona  is  the  most  effec¬ 
tive. 


3.  Under  the  conditions  of  the  experiments,  it  is  not 
:  ely  that  the  positive  ions,  which  are  presented  in  the 
space  between  the  electrodes  as  a  result  of  corona  discharge, 
will  act  as  nuclei  for  the  frost  crystals.  The  reason  is 
nucleation  requires  the  air  to  be  supersaturated.  The 
increase  in  frost  deposition  under  the  action  of  a  corona 
discharge  thus  cannot  be  attributed  to  nucleations  effected  by 
ions. 


4.  Hie  formation  of  water  vapor  clusters  centered  at 
the  ions  can  neither  account  for  the  increase  in  frost  deposi¬ 
tion.  Calculations  (see  Appendix  B)  indicated  that  a  cluster 
of  the  order  of  105  water  molecules  for  each  ion  is  required 
for  the  postulate  to  be  valid.  Such  clusters  hardly  exist 
under  the  test  conditions. 

5.  The  test  on  positive  wire-screen  corona  strongly 
suggests  that  impinging  electric  wind  produced  by  the  corona 
discharges  is  the  crucial  mechanism  which  effects  the  changes 
in  frost  formation. 


6.  A  possible  mechanism  to  explain  the  ridges  observed 
in  the  frost  strip  is  the  secondary  flow  resulting  from 
electrohydrodynamic  instability.  The  secondary  flow  consists 
of  Goertler  vortices  having  the  shape  shown  in  Fig.  1,  their 
axes  being  parallel  to  the  basic  flow  direction.  The  inflow¬ 
ing  parts  of  Goertler  vortices  would  enhance  the  frost  deposi¬ 
tion  and  outward  flowing  parts  would  retard  it,  so  that  the 
frost  thickness  on  the  wall  would  vary  periodically  in  the 
direction  normal  to  the  plane  of  the  flow. 


1.3  DESCRIPTION  OF  THE  PROBLEM 

In  the  light  of  the  observations  and  consents  made  in  the  preceding 
section,  the  problem  of  interest  can  now  be  described:  the  determina¬ 
tion  of  heat  and  mass  transfer  characteristics  of  impinging  electric 
wind  under  frosting  conditions.  Figure  2  shows  the  overall  schematic 
of  the  physical  model  of  the  present  study.  The  frosting  surface  is 
assumed  to  extend  to  infinity  and  maintained  at  a  constant  subfreezing 
temperature  Tw.  The  wire,  located  at  a  distance  d  from  and  parallel 
to  the  frosting  surface,  is  subject  to  a  positive  potential  V.  The 
temperature  and  the  vapor  mass  fraction  of  the  surrounding  air  are 
I,*  and  w^  respectively.  The  frost  layer  thickness  is  designated  by  Y 
and  the  frost-air  interface  temperature  by  Tf . 

When  the  applied  potential  of  the  wire  exceeds  the  critical  poten¬ 
tial  at  which  corona  starts,  ionization  takes  place  in  the  immediate 
vicinity  of  the  wire  where  an  intense  field  exists.  The  positive  ions 
thus  fonaed  move  away  from  the  wire  under  the  influence  of  the  electric 
fields.  As  they  are  driven  towards  the  cathode  plate,  the  ions  make 
many  collisions  with  the  neutral  molecules  so  that  a  streaming  of  the 
air  results  through  momentvm  transport.  As  a  consequence  of  the  impinge¬ 
ment  of  the  induced  air  stream,  which  is  at  the  same  temperature  and 
specific  humidity  as  the  surrounding  air,  heat  and  mass  are  transported 
between  the  ambient  air  and  the  frost  layer. 

The  problem  of  frost  formation  has,  in  practice,  two  phases.  The 
convective  region  determines  the  convective  neat  and  mass  transfer 
between  the  surface  of  the  deposited  layer  and  the  external  stream. 

On  the  other  hand,  the  frost  properties,  the  frost  layer  thickness  as 
well  as  the  surface  conditions  of  the  frost  layer  such  as  temperature 
and  specific  humidity  are  determined  by  the  transport  mechanisms  in 
the  frost  layer  itself.  Of  course,  the  two  phases  are  interrelate!, 
tor  the  temperature  and  material  distributions  in  the  convective  region 
depend  on  the  surface  conditions  of  the  frost  layer,  whereas  the  deter¬ 
mination  of  the  transport  processes  in  the  frost  layer  requires  a 
specification  of  the  boundary  conditions  imposed  by  the  convective  flow 
on  the  frost-air  interface.  Although  the  main  object  of  this  study, 
as  has  already  been  pointed  out,  is  t..e  convective  transport  processes 
in  connection  with  the  Impinging  electric  wind,  a  chaper  (Chapter  II) 
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Fig*  1  Expected  Goertler  vortices  in  two-dimensional 
stagnation  flow  under  the  impinging  electric 
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Ower-all  schematic  of  the  physical  model 
I  :  Electric  wind  region 
I  :  Stagnation  region 
I  :  Transition  region 
V  :  Wall  jet  region 


is  devoted  to  the  study  of  transient  heat  conduction  in  an  idealized 
frost  layer.  It  is  intended  to  bring  out  the  relative  importance  of 
the  parameters  which  affect  the  frost  layer  thickness  and  the  tempera¬ 
ture  of  the  frost-air  interface.  No  estimate  of  the  variations  of  the 
frost  properties  with  time  is  attempted . 

In  analogy  with  the  corresponding  problem  for  impinging  free 
jet8»9>10». the  flow  field  of  the  impinging  electric  wind  can  be  divided 
into  four  regions,  if  one  follows  the  flow  development  along  the  center- 
line  of  the  wind  to  the  stagnation  point  and  then  from  the  stagnation 
point  along  the  solid  boundary.  The  air  stream  starts  out  from  a  very 
narrow  region  close  to  the  wire  and  is  accelerated  by  the  electric 
field  before  finally  being  retarded  by  the  plate.  This  zone  may  be 
called  the  electric  wind  proper.  At  certain  distance  from  the  plate, 
the  jet  of  electric  wind  begins  to  be  affected  by  the  plate.  The 
velocity  decelerates  rapidly  along  the  centerline  of  the  jet  as  the 
plate  is  approached  and  the  pressure  grows  to  a  maxi  mum  value  at  the 
stagnation  point.  The  jet  is  then  deflected  and  spread  over  the  plate. 
This  is  the  zone  of  the  two-dimensional  stagnation  flow.  In  the  inme- 
diate  neighborhood  of  this  stagnation  zone  lies  a  transition  zone  which 
ends  at  a  distance  Xo  from  the  stagnation  point  along  the  plate.  Beyond 
Xo  follows  the  zone  of  two-dimensional  wall  jet  which  flows  parallel 
to  the  plate. 

In  the  case  of  impinging  free  jet,  the  entire  flow  Held  will  be 
turbulent  in  most  practical  examples  owing  to  the  low  cr.1  tical  Reynolds 
number  for  a  free  mixing  layer.  No  attempt  is  made  in  Tihis  study  to  con¬ 
sider  the  turbulent  flows  created  by  the  impinging  electric  wind,  since 
the  velocities  obtainable  from  the  corona  wind  in  air  are  very  small, 
of  the  order  of  a  few  feet  per  second.  Moreover,  it  is  felt  that  the 
understanding  of  the  ideal  laminar  flow  field  appropriate  to  the  actual 
flow  is  a  necessary  preliminary  to  the  understanding  of  the  actual 
aerodynamic  processes  associated  with  the  impinging  electric  wind. 

The  heat  and  mass  transfer  processes  between  the  impinging  air 
stream  and  the  frost  layer  is  in  reality  non-stationary,  for  the  tem¬ 
perature  of  the  frost-air  interface  is  a  function  of  time.  The  extent 
of  the  deviation  of  the  instantaneous  heat  and  mass  transfer  rates  at 
the  surface  from  the  quasi-steady  values  corresponding  to  the  instan¬ 
taneous  surface  temperature  at  a  given  instant,  depends  on  the  rapidity 
of  the  variation  of  the  surface  temperature  and  the  response  character¬ 
istics  of  the  boundary  layer  to  a  sudden  change  of  the  boundary  condi¬ 
tions.  The  process  of  frost  formation  in  ambient  air  is  in  general 
very  slow,  the  variation  of  the  fro3t  layer  surface  temperature  being 
at  most  of  the  order  of  a  few  *F  per  minute,  and  the  response  of  the 
boundary  layer  in  air  to  a  sudden  variation  of  the  boundary  conditions 
is  generally  very  fast.  Consequently,  it  will  be  assumed  that  quasi¬ 
steady  conditions  prevail  in  the  impinging  flow  of  an  electric  wind. 


Despite  the  fact  that  the  experimental  study  was  carried  out  on  a 
vertical  plate,  the  effect  of  free  convection  is  disregarded  in  the 
analysis.  In  virtue  of  the  smallness  of  the  V'  cities  obtainable  from 
the  corona  wind,  this  assumption  may  not  be  jut  led.  The  aim  of  this 
investigation,  however,  is  to  throw  light  on  the  actions  of  electric 
wind  upon  heat  and  mass  transfer  mechanisms,  it  appears  logical  to  con¬ 
sider  electric  wind  effects  alone. 

The  problem  of  electrohydrodynamic  instability  and  its  influence 
on  heat  and  mass  transfer  processes  is  beyond  the  scope  of  the  present 
investigation,  nevertheless,  it  is  worthwhile  to  make  a  few  remarks 
in  this  regard.  Two  types  of  electrofaydrodynamic  instability  phenomena 
may  be  distinguished:  those  caused  by  surface  forces,  and  those  induced 
by  body  forces.  The  former  take  place  at  the  interface  of  two  media 
having  different  dielectric  properties  (usually  liquid-gas  or  liquid- 
liquid  interface),  while  the  later  occur  in  the  bulk  of  the  medium 
itself.  Another  difference  between  these  two  classes  of  instability 
is  that  in  the  first  case  the  media  are  electrically  neutral  except 
possibly  at  the  interface,  whereas  in  the  second  case  the  presence  of 
space  charges  in  the  medium  is  a  must  for  the  instability  to  occur. 
Theories  for  the  first  category  of  EHD  instability  have  been  developed 
by  Malkus  and  Veronis11,  by  Melcher12,  by  Taylor  and  McEwan13,  and 
more  recently  be  Lee  and  Choi14.  The  theory  for  the  second  category 
of  instability  is  still  in  its  infancy.  Stach15,  in  consideration  of 
flows  in  two-dimensional  flat  channel,  assumed  the  instability  to  be 
in  the  form  of  propagating  waves  (i.e.  the  Tcllmien-Schliching  waves) 
and  jumped  to  the  conclusion  that  the  problem  on  hand  was  analogous  to 
that  of  stratified  flows.  He  fed. led,  however,  to  realize  that  for 
stratified  flows  the  density  distribution  of  the  fluid  was  such  that 
the  arrangement  was  stable,  while  for  the  flows  he  was  considering  the 
arrangement  had  to  be  unstable  in  order  for  the  instability  to  set  in. 

A  different  approach,  proposed  by  Chuang  and  Welkoff1®  in  consideration 
of  quiescent  fluid  between  two  concentric  cylinders,  is  to  assume  that 
the  Instability  takes  the  form  of  stationary  waves  (i.e.  Taylor -Goertler 
vortices).  Although  it  has  been  observed  experimentally  that  stationary 
periodic  waves  exist  in  annular  flows17,  further  experimental  evidence 
is  required  to  verify  the  validity  of  the  postulate. 


CHAPTER  II 


TRANSIENT  HEAT  CONDUCTION  IN  THE  FROST  LAYER 


2.1  FORMULATION  OF  THE  FROBLEM 


Consider  a  flat  plate  precooled  to  a  constant  subfreezing  tempera¬ 
ture  Tw  and  located  at  yO  (Fig.  3) •  When  the  plate  Is  exposed  to  the 
moist  ambient  air  at  constant  temperature  Toe  and  vapor  mass  fraction 
wioo,  the  frost  starts  to  grow  on  the  surface  of  the  plate.  The  frost 
formation  occurs  progressively  at  the  frost-air  interface  as  a  result 
cf  heat  and  mass  transfer  frcm  the  surroundings.  Heat  flows  by  con¬ 
vection  from  the  ambient  air  to  the  frost-air  interface  and  by  con¬ 
duction  through  the  frost  layer  to  the  flat  plate. 

The  transient  process  of  heat  transfer  effected  by  the  growth  of 
the  frost  layer  is  extremely  complicated.  Since  the  temperature  of  the 
frost-air  interface  varies  with  the  time  as  the  frost  grows,  the  heat 
and  mass  transfer  coefficients  sure  not  constant,  but  functions  of  time. 
In  addition,  the  density  and  the  thermal,  conductivity  of  the  frost 
layer  vary  unpredlctably  through  a  wide  range.  It  appears  that  an  exact 
mathematical  solution  to  problem  is  extremely  difficult  if  all  the 
numerous  variables  are  considered. 

In  the  present  study,  the  problem  is  idealized  as  follows:  It  is 
assumed  that  the  frost  forms  uniformly  along  the  surface  so  that  the 
problem  is  one -dimensional.  It  is  also  assumed  that  the  frost  proper¬ 
ties  and  the  heat  and  mass  transfer  coefficients  do  not  vary  with  time. 
Furthermore,  all  water  vapor  transferred  by  convection  from  the  ambient 
air  to  the  frost-air  interface  is  assumed  to  deposit  as  frost  and  is 
responsible  for  the  increment  in  fro  it  thickness.  Special  attention  is 
directed  towards  the  time -dependence  of  the  frost  layer  thickness  and 
of  the  frost-air  interface  temperature,  particularly  the  latter. 

Thus  the  problem  of  frost  formation  can  be  formulated  as  a  bound¬ 
ary  value  problem  in  v  hich  the  governing  equation  is  the  one -dimensional 
heat  conduction  equation  for  the  solid  phase  without  heat  generation 
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subject  to  the  boundary  conditions  that 
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at  y 


T  =  T,. 


(2.2) 
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(2.3) 


at  y  -  Y(t):  k  «  PfLf^J  +  h(T„  -  ?‘f) 

ay  at 

Pf  f£  =  Pha^»  ‘  Wlf  >  (2‘4) 

and  the  Initial  condition  that 

Y(0)  -  0  (2.5) 

where  a  ■  thermal  diffusivity  of  the  frost, 

h  ■  convective  heat  transfer  coefficient, 

*m  ■  convective  mass  transfer  coefficient, 

Lf  =  latent  heat  of  sublimation  of  the  frost, 
t  =  time, 

T  =  temperature, 

Tf  *»  temperature  at  the  frost-air  interface, 
wif  ■  vapor  mass  fraction  at  the  frost-air  interface, 

Y  ■  thickness  of  the  frost  layer, 
p  *  density  of  the  moist  ambient  air 
pf  ■  density  of  the  frost 
and  k  *  thermal  conductivity  of  the  frost. 

In  the  boundary  condition  (2.3),  the  superheat  of  the  vapor  is  neglected. 

To  complete  the  formulation  of  the  problem  it  is  necessary  to 
specify  the  relation  between  wif  and  Tf.  Since  (T-  Tw)  is  usually  small, 
it  can  be  assumed  that  (l8) 

(wi-wiw)  -  b(T-Tw) 

where  b  is  a  constant  of  proportionality  and  wiw  is  the  vapor  mass 
fraction  of  the  saturated  air  at  Ty.  In  the  present  case 


b  =  (wif-viw)/(Tf-Tw) 


Writing 


Eq.  (2.4)  becomes 


Wise  -  Vif  =  (wiarWiy)  -  (wif-Wlw) 

*  (w«-wiw)  -  b^f-T*) 


Pf  H  ■  phm  [(wio^wiv)  -  b^f-T*)] 


or 


—  *  Ci  -  C2Tf 
dt 


(2.6) 


Here  the  constants  Ci  and  C2  ere  defined  as 


Ci  -  ^toUvi^-viw)  +  bTw] 


(2.7) 


and 


C2  -  Ph^b/pf 


(2.8) 


Substituting  Eq.  (2.6)  in  (2.3),  there  is  obtained  an  alternative  form 


for  condition 


;q. 

(2.3). 


k  |E  -  PfCiLf  -  PfC2LfTf  +  h(T«rTf) 

<y 


Or,  after  rearrangement 


k  |2  =  C3  -  C4Tf 
V 


(2.9) 


where 


C3  =  PfCiLf  +  hTfl 


and 


C4  ~  P£C2L£*  h 


The  boundary  value  problem  as  formulated  above  is  written  down 
once  more  for  reference. 

Governing  equation: 


Boundary  conditions: 


to  « .  tor 

St 

(2.10) 

at  yO: 

T  »  Tw 

(2.11) 

at  y*Y: 

k  -  C3 

-  C4,Tf 

(2.12) 

a* -Ci 

dt  1 

-  C2Tf 

(2.13) 

Initial  condition: 

Y(0)  «  0  (2.14) 

The  essential  features  of  the  problem  are  the  existence  of  a 
moving  surface  of  separation  between  the  solid  and  gaseous  phases  and 
the  dependence  on  time  of  the  temperature  of  this  moving  surface.  The 
simplicity  of  Eq.  (2.10)  is  misleading.  The  time-dependent  boundary 
conditions  (2.12)  and  (2.13)  render  the  problem  one  of  considerable 
difficulty.  It  is  only  recently  that  an  exact  solution  of  Eq.  (2.10) 
which  satisfies  any  arbitrary  initial  and  boundary  conditions  ha*  been 
developed.  By  application  of  the  two-sided  Laplace  transform  Portnov19 
was  able  to  obtain  an  exact  solution  in  the  form  of  Poisson  Integral. 


2.2  PORTNOV'S  METHOD  OF  SOLUTION 

In  this  section,  Portnov's  method  of  exact  solution  to  Eq.  (2,10) 
in  the  form  of  Poisson  integral  is  briefly  described.  For  details  and 
rigorous  proofs  of  the  procedure,  reference  should  be  made  to  Portnov's 
original  paper. 

Letting  t  »  t  and  y  ■  ^?Y(t),  it  follows  that 

Y1  (t)  _y_  to  _  3t  _  Y'(t)*3t 

dt  dr  dr  dr  ^(t)  ^  dg  dr  y(t) 
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and 


_  1  c^T 

ay2 "  ¥&) 

Thus  Eq.  (2.10),  in  terms  of  the  new  set  of  variables  t  and  £,  reduces 
to  the  form 

Y2^)  ^  “  Y(t)Y»  (t)  £  ||  -  |p  (2.15) 

where  prime  denotes  differentiation  of  a  function  with  respect  to  its 
argument. 

Now,  applying  to  Eq.  (2.15)  the  two-sided  Laplace  transform20 
with  respect  to  5(t  being  treated  a  parameter)  it  is  obtained 

l*(T)|+KT)r(T)p|.A  (2.16) 

where  p  is  a  new  independent  variable  introduced  in  the  Laplace  trans¬ 
form. 

The  general  solution  to  Eq.  (2.l£)  is 

T(p,x)  =  exp[p2T/Y2(t)]F(p/Y)  (2.17) 

where  F(p/Y)  is  an  arbitrary  function  of  its  argument.  If  the  function 
F  is  represented  as 

F(p/Y)  =  Fi(p/Y)  +  F2(p/Y)  (2.18) 

and  it  is  assumed  that  the  functions  Fx  and  Fa  are  both  Laplace  trans¬ 
form  images  of  the  arbitrary  one-sided  functions  cpi  and  92,  i.e. 

Fi(p/Y)  pi  «Ri(YC)U(Y0 


Fa(p/Y)  <pa(Y0u(-Y0 

where  fi  denotes  the  Laplace  transform  and  U  denotes  the  unit  function, 
inverse  transform  of  Eq.  (2.17)  yields 
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T(£,t)  -  exp  y2^-Z^-]  j  Yz)u(>/a  Yz) 

+  <te(^a  Yz)u(-»/a  Yz)  I 

I  2  vi 


(2.19) 


T({*x)  ■  £  I  /o"-*  [•  ^Er2]  *(/a  *z)  7z 

+  £  «p  ,*(/.  Iz)  m  | 


In  terns  of  the  original  variables,  Eq.  (2.19)  reads 


and  Eq.  (2.20)  becomes 


t  /*  00 

T(y,t)  |/Q  exp  j  q>i(y)dy 


■£exp  [• 


(y-y)2' 


qte(y)dy 


(2.20) 


T(y,t)  "TV  ["  ^^"jf(y)^  (2.21) 


(2.22) 


where  a  »  2  *i/a7,  jr  ■  Y  '/az  and  f(y)  ■  9i(y)U(y)  +  9a(y)u(-y).  Thus 
Eq.  (2.22)  represents  the  temperature  distribution  in  an  infinite  slab 
with  an  initial  temperature  of  f(y).  More  specifically,  the  functions 
<Pi(y)u(+y)  and  qfe(y)U(-y)  represent  the  initial  temperature  for  y  >  0 
and  y  <  0  respectively.  For  the  case  of  frost  formation,  both  9x  and 
qp2  represent  fictitious  initial  temperatures  since  the  frost  layer  does 
not  exist  at  t  ■  0. 

The  unknown  functions  Y,  cpi  and  92  in  Eq.  (2.22)  are  to  be  deter¬ 
mined  by  using  the  initial  and  boundary  conditions  for  the  specific 
problem  under  consideration.  The  general  method  is  to  assume  that 
these  functions  can  be  expanded  in  power  series 
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00 

<Pi(w)  ■  E  cpinw0 

(2.23a) 

n^D 

00 

<P2(“>)  “  E  <P2n^ 

(2.23b) 

n=0 

*(*)  "  E 

(2.23c) 

n=l 

where  Eq.  (2.23c)  satisfies  the  initial  condition  (2.14).  Then  Eq.  (2.22) 
with  Y,  cpi  and  $2  substituted  from  Eqs.  (2.23)  is  used  in  each  of  the 
boundary  conditions,  the  resulting  equations  being  repeatedly  differen¬ 
tiated  with  respect  to  o  and  the  equations  thus  obtained  being  evaluated 
at  cr  **  0.  Each  set  of  the  resulting  stlim.il taneous  algebraic  equations 
may  be  solved  for  the  successive  coefficients  of  the  unknown  series 
(2.23).  Once  these  have  been  found,  Eq.  (2.22)  becomes 


where 


Eq.  (2.24)  can  also  be  written  a u 

T(y,t)  =  i  f  in  erfc(-  +  (-l)n  cpsn  in  erfc(j)] 

nO 

where 

ln  erfo(y)  =  i-  f“  **  bltil  dy 
*'/n  »: 


(2.24) 


(2.25) 


denotes  the  nth  order  repeated  integral  of  the  complementary  error 
function. 
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2.3  SERIES  SOLUTION  OF  THE  FROSTING  PROBLEM 


Portnoy  derived  Eq.  (2.22)  without  using  any  boundary  conditions. 
Since  the  integral  can  be  shown  to  be  a  solution  of  Eq.  (2.10)  tending 
to  f(y)  as  t  approaches  zero,  this  integral  form  of  solution  has  been 
assumed  to  apply  in  problems  involving  melting  and  freezing  of  finite 
slabs  by  Jackson21  and  by  Westphal22.  In  the  present  study  the  same 
method  is  applied  to  frosting  problem. 

To  determine  the  unknown  functions  Y,  <pj.  and  $2,  use  is  made  of 
the  boundary  condition  (2.11)  through  (2.13).  When  Eq.  (2.22)  is 
substituted  into  Eq.  (2.11),  it  is  found 

Tw  "TIL  e"P2  +jf„  (2.26) 


where 


p* 


Substituting  Eq.  (2.23)  into  Eq.  (2.26)  and  letting  o  -»0,  the  following 
relation  is  obtained 


Tw  -  £  (9io  +  920) 


(2.27) 


Differentiation  of  Eq.  (2.22)  with  respect  to  y  gives  the  tempera¬ 
ture  gradient 

I  'VTal  fo  ['  “P 


/.  [-  eicp  [‘^-] 


(2.28) 


Substitution  of  Eqs.  (2.22)  and  (2.28)  into  Eq.  (2.12)  yields 
2h[  jT/^e-72  cpi(7a+Y)dr  *J^^7«T72  q>2(7cr+Y)d7  ] 


=  4  C3-(rC4  [£  e"7  9i(7<r+Y)d7  <P2(7cr+Y)d7] 

(2.29) 
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where 


,2  =  (y-y)2 


Substituting  Eqs.  (2.23)  into  Eq.  (2.29)  and  letting  a  -»0,  it  is  found 


f”  _V2  /*‘Yl-72 

910^  7e  7  d7  +  920^  7e  d7  »  0 


(2.30) 


To  determine  Yi,  substitute  Eq.  (2.22)  into  Eq.  (2.13)  to  yield 


2  -  ^  ■  Cx  -  ^2-f*  /*  e“7  <Pi(70+X)d7  +  f  i"1*  92(70-»Y)d7l 

a  da  fn  I <*  J  J 

Substituting  Eqs.  (2.23)  into  Eq.  (2.31)  and  letting  <r  -»0,  it  is 
obtained 

Yi  -  0 

Using  Eq.  (2.32)  in  Eq.  (2. 30),  it  is  found 


9l0  -  920  ■  0 


(2.31) 


(2.32) 


(2.33) 


Thus,  from  Eqs.  (2.27),  (2.32)  and  (2.33 )»  the  first  coefficients  of 
the  unknown  functions  Y,  cpi  and  92  are  found  to  be 


<Pio  ■  920  -  Tw  ,  Yi  ■  0 


The  next  set  of  coefficients  Y2,  911  and  921  are  obtained  by  taking 
first  derivatives  of  Eqs.  (2.26),  (2.29)  and  (2.31),  substituting 
Eqs.  (2.23)  into  the  equations  thus  obtained,  letting  <r  -*0  and  evaluat¬ 
ing  the  resulting  definite  integrals. 

Ike  first  derivative  of  Eq.  (2.26)  with  respect  to  <7  is 

0e"p2  <pi(l)((3o)de  +  Se"^  92ri)(P*)<tf  J 

°  (2.3*0 

where  superscripts  indicate  the  order  of  differentiation  with  respect 
to  cr.  Substituting  for  91  and  92  from  Eqs.  (2.23)  and  letting  c  -»0, 
Eqs.  (2.34)  becomes 
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Or,  after  evaluation  of  the  definite  integral. 


cpii  -  $21  *  0  (2.36) 


The  first  derivative  of  Eq.  (2.29)  with  respect  to  a  is 


2k 7(74Y(1))e"72  <Rl(1)(7(rfY)d7  +  A^Ae"*2  $i(0) 

+  y~\(7-H f^)®”7*  $2^(7<r>¥)d7  -  $b(0)] 


•Jit  C3-C4[^  e"72  $i(7«r+Y)d7  +  j  e"7*  $2(70^)17} 

-  *  C4[  /  (7+Y^1^)e-72  q>i^1^(7o+Y)d7  +  A^e"**  q*(0) 

L  /-A 

+f  (7’*Y^1^)e“72  $2^1^(70+Y)d7  -  A^e"*2  qte(0)J 


(2.37) 


where 

A  «  Y/a  . 


Substitution  of  Eqt>.  (2.23)  into  Eq.  (2.37)  and  letting  <r  -»0  yields 


i(<Pii+$2i)  -  i(C3-C4Tw)  (2.38) 

k 

In  obtaining  Eqs.  (2.38),  cpio,  $20  and  Yi  have  been  substituted  from 
the  previous  results. 


To  obtain  Ya,  differentiate  Eq.  (2.31 )  with  respect  to  a 


00 

Cl  ~  e”72  <Pl^7CT+Y^d7  +  f  e~7"  <Ba(7«r«T)d7 

(7+Y^^)e"72  9i(1)(7a+Y)d7j 
+  A^e"^  cri(O)  +  /  (7+Y^1))e“72  «P2^1^(7<rfY)d7 

J —to 

-  CP2(0)J 


(2.39) 


Substituting  Eqs.  (2.23)  into  Eq.  (2.39)  and  letting  o  -*0,  it  is 
obtained 

4aY2  «  Ci  -  (cpio  +  920) 

2 

which,  on  substituting  for  <pi0  and  920,  yields 


Ct  “CpTi 


(2.4o) 


Thus  the  second  set  of  coefficients  Y2,  9n  and  921  are 

<Pix  =  <Pax  =  ^(C3-C4Tw) 
k 


Ya  =  (Ci-C2Tw)/4a 

To  obtain  higher  order  coefficients  of  the  unknown  functions  Y,  91  and 
92)  it  is  necessary  to  continue  in  the  manner  described  above,  i.e. 
obtain  higher  order  derivatives  of  Eqs.  (2.26),  (2.29)  and  (2.31),  let 
cr  -» 0  and  evaluate  the  definite  integrals  appearing  in  the  resulting 
equations.  The  first  six  sets  of  coefficients  are  listed  in  Appendix  C. 


2.4  FROST  LAYER  THICKNESS 

Substituting  in  Eq.  (2.24)  the  coefficient  Yn's  from  Appendix  C, 
there  results  an  expression  for  the  frost  layer  thickness. 
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Y(t)  =  JL  bmCwxcorWiwJt  -  ^i)  2  [pfl^m(wi--wiw) 

+  h(To,rTv)J  (wi(*rwiw)t2  +  (p^c)  *  (Wl»"Wlw) 

•  |  (pbLfl^a'Hi)  j^pI*fhjj(wio^wiy )  +  h(T„-T^)J 

•  (wi.-wiw)  ^pLfhaCviorWiw)  +  htT^T*)  j2  1 13  +  ••• 

(2.41) 


It  is  noted  that  the  higher  order  terns  are  due  to  the  unsteadiness  of 
wif.  A  few  remarks  can  be  made  concerning  this  equation.  For  small 
tines,  higher  order  terms  in  the  equation  are  negligible  as  compared 
with  the  first  term.  Hence  the  frost  layer  grows  linearly  with  time. 
The  rate  of  change  of  the  frost  thickness  is  given  by  the  coefficient 
of  the  first  term,  which  represents  the  ratio  between  the  mass  flux 
of  the  water  vapor  to  the  frost  surface  by  convection  and  the  frost 
density.  This  is  in  accordance  with  the  boundary  condition  (2.4),  for 
at  tO,  vif  is  equal  to  wiv.  As  t  increases,  the  higher  order  terms 
became  more  and  more  significant  and  the  rate  of  increase  of  the  frost 
thickness  tapers  off  gradually.  The  frost  thickness  asymptotically 
approaches  its  maximum  value  as  t  tends  to  infinity,  a  direct  conse¬ 
quence  of  the  fact  that  Eq.  (2.24)  can  be  shown  to  converge  uniformly 
for  all  crl9. 

Physically,  the  mass  transfer  process  can  not  go  on  indefinitely, 
since,  in  principle,  as  soon  as  the  concentration  of  water  vapor  at 
the  surface  of  the  deposited  layer  reaches  a  value  equal  to  that  of  the 
ambient  air,  an  equilibrium  state  is  then  established  and  no  further 
mass  transfer  is  possible.  As  for  the  phenomenon  of  frost  formation 
two  cases  may  be  distinguished23.  Where  the  air  dew  point  is  below 
32°F,  deposition  of  frost  will  stop  when  the  frost  surface  reaches  the 
dew  point.  When  the  dew  point  of  the  air  is  above  32 °F,  the  frost 
surface  temperature  can  not  exceed  32 °F.  In  the  latter  case,  it  is 
generally  observed  that  after  a  temperature  of  32  °F  is  attained  at 
frost  surface  a  steady-state  condition  is  reached  in  which  water  vapor 
will  condense  as  liquid.  This  liquid  will,  if  the  frost  is  porous, 
soak  into  the  frost  layer  and  eventually  freeze  resulting  in  a  con¬ 
tinuous  increase  in  the  density  as  well  as  the  thermal  conductivity  of 
the  frost  layer  approching  the  value  of  solid  ice. 

From  the  second  term  of  Eq.  (2.4l),  it  is  noted  that  the  heat 
transfer  rate  h(T«rTw)  between  the  wall  and  the  surroundings  has  the 
effect  of  reducing  the  amount  of  frost  deposited.  If  the  temperature 
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difference  between  the  wall  and  the  free  stream  or  the  surrounding  air 
becomes  extremely  large,  considerable  deviation  of  the  actual  frost 
deposited  from  the  calculated  mass  transfer  by  convection  phm(wi«rwiw) 
is  expectable.  It  has  been  reported  by  Barron24  that  the  magnitude  of 
the  mass  transfer  rate  measured  experimentally  at  cryogenic  temperature 
was  approximately  one  order  of  magnitude  less  than  the  theoretically 
predicted  correlation.  Barron  tried  to  interpret  the  discrepancy 
between  theory  and  experiment  in  terns  of  the  resistance  to  diffusion 
of  water  molecules  by  the  presence  of  frost  particles  which,  supposedly 
formed  by  nucleation,  were  observed  to  exist  in  the  boundary  layer  in 
his  investigation.  Equation  (2.4l)  seems  to  suggest  that,  in  addition 
to  the  appearance  of  the  visibile  frost  particles,  large  values  of 
(T„rTw)  may  also  be  responsible  for  the  discrepancy. 


2.5  FROST-AIR  INTERFACE  TEMPERATURE 

The  expression  for  the  frost-air  interface  temperature  as  a  func¬ 
tion  of  time  up  to  the  fourth  term  may  be  obtained  by  substituting  in 
Eq.  (5»25)  the  coefficients  of  91  and  92  listed  in  Appendix  C  and 
evaluating  the  resulting  expression  at  y=Y,  namely. 


Tf(t)  =  Tw  +  |  [pLfh^wioe-wi*)  +  hCT^-T*)] 

-  i  (Ya2+  |y3  )  (pbIThm-tti)(wlw-wXw) 

•  [pLfhmCwiorWiv)  +  h(T«rTw)  J 

*2*  (j,  I  +  jlt£+4I!) 

3  V 12  cr  12  a3  ifi  a5  ) 


(obLfhm-Hi) 


90k 


(wioo-wiw) 


j  [  (l^t)  ] 

j^PkfhfljCwiflj-Wlw)  +  k(^oo"^-w  )  J 

[!♦£  (&)<-~*>]  *£  (&) 


[pLfhm(wlorW1w)  +  2  J  +  •• 


(2.42) 


Without  recourse  to  mmerlcal  computations,  it  is  seen  directly 
from  the  equation  that  at  t  ■  0,  Tf  »  Tw.  Obviously  this  is  true,  for 
then  the  frost  thickness  is  zero.  To  a  firct  oproxlmatlon,  the  frost- 
air  interface  temperature  is  proportional  to  iue  frost  layer  thickness 
as  well  as  the  total  heat  flux  to  the  frost  surface.  When  Y  and  the 
latent  heat  of  sublimation  of  the  frost  are  small  as  compared  with  con¬ 
vective  heat  flux,  the  latter  is  an  important  parameter  in  deteimining 
the  instantaneous  frost  surface  temperature.  It  is  also  noted  from  Eq. 
(2. 42)  that  the  variation  in  time  of  Tf  is  essentially  effected  by  the 
growth  of  the  frost  thickness  as  expected. 

It  is  of  interest  to  mention  here  the  temperature  distribution  in 
the  frost  layer.  An  expression  for  such  temperature  distribution  is 
readily  reduced  from  Eq.  (2.42),  namely, 

T(y,t)  -  Tw  +  £  ^pLfhntfrieo-viw)  +  Mt^-T*)] 


-  ^(yo^+ly3)  -^-(pbLfhm-HiJCwieo-wiw) 

3  pfak 

•  [pLfh^wxeo-vi*)  +h(T0O-Tw)J  +  •••  (2.43) 

where  y  <  Y.  Clearly  this  equation  indicates  that  at  any  instant,  the 
frost  temperature  increases  from  Tw  at  the  plate  to  Tf  at  the  frost-air 
interface.  But  the  increase  is  non-linear  because  of  the  presence  of 
higher  order  terms  in  powers  of  y.  The  rate  of  increase  in  T  decreases 
with  increasing  y.  On  the  other  hand  the  temperature  at  a  fixed  point 
in  the  frost  layer  varies  with  time.  The  third  term  in  Eq.  (2.43)  is 
the  most  Important  tern  in  the  determination  of  the  temperature  history 
throughout  the  frost  layer.  Consider  the  point  y  -  Yi  and  suppose  the 
frost  starts  to  form  at  t«0.  When  the  frost-air  interface  reaches  the 
point  y  »  Y(tt)  *  yx,  the  temperature  at  yi  takes  the  value  Tf(tj.).  As 
the  frost  continues  to  grow,  i.e.  for  t  >  ti,  T(yi)  decreases  with  time 
until  a  steady  state  is  reached  where  Tf  assumes  a  value  corresponding 
to  the  maximum  physically  possible  frost  surface  temperature. 

2.6  THE  HEAT  TRANSFER  RATE 

To  detexmine  the  heat  transfer  rate  across  a  surface  parallel  to 
the  plate  per  unit  area  of  the  surface,  Eq.  (2.43)  must  be  differentiated 
with  respect  to  y.  Hie  heat  flux  at  a  point  y  and  time  t  is  as  follows: 
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> 


q  =  -  k(|^=  -  [p^hm(wioo-wiw)  +  hfToo-Ty)] 

+  v  (a‘  +2ya)  — ^|~  (pbLflim+h)(wx00-wiw) 
pfk  a 

*  ^PLf>hm(wioo"wiw)  +  h(Tx-Tw^J 

-  i(±  g*4gy+gy\  (gbiyj^tfa) 

3  \12  £  y  16^  /  90 

*  [(^^^2(wiocrWiw)2(pbLfhm-»h)] 

*  |  [pLfhm^ioo-wiw)  +  ^(Too-Tv)] 


+  4  ”  [P4fh,,(»i-«iw)  +  M  W]  | 

(2.1*) 

From  this  equation  it  may  be  concluded  qualitatively,  that  the  instan¬ 
taneous  heat  flux  in  the  frost  layer  decreases  with  y.  3his  is  because 
the  local  heat  storage  in  the  frost  layer  decreases  with  decreasing  y 
as  does  the  local  frost  temperature ,  and  a  certain  amount  of  heat  is 
expected  to  be  released  and  transferred  in  the  negative  y  direction, 
resulting  in  a  heat  flux  increase  in  that  direction.  The  minimum  flux 
occurs  at  the  frost-air  interface  and  is  equal  to  plfh^wiof-wif )  +  h( 

+  h(TorTf)  as  given  by  the  boundary  condition  (2.3).  The  maximum  heat 
flux  is  at  the  plate.  The  second  term  in  Eq.  (2.44)  accounts  for  the 
fact  that  the  frost  layer  has  insulating  effect  upon  heat  transfer. 
Another  point  worth  noting  is  the  unsteadiness  of  the  local  heat  flux. 
It  is  apparent  from  Eq.  (2.44)  that  the  local  heat  flux  decreases  with 
time. 
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2.7  CLOSURE 


The  analysis  developed  In  this  chapter  has,  essentially,  no  direct 
bearing  with  the  remainder  of  this  dissertation  and  is  included  here  as 
an  independent  topic.  As  has  been  emphasized  in  Section  1.3 ,  the  tran¬ 
sient  heat  conduction  in  the  frost  layer  is  an  important  phase  of  the 
problem  of  frost  formation.  It  is  thus  felt  that  this  subject  deserves 
a  discussion  even  in  an  investigation  which  concerns  primarily  with  the 
convective  transport  processes  associated  with  frost  formation. 

Because  of  the  complex  nature  of  the  problem  no  satisfactory  theory 
has  been  developed.  The  present  study  represents  a  new  but  preliminary 
attempt  to  bring  out  the  important  parameters  which  affect  the  heat 
transfer  rate  across  the  frost  layer  and  the  temperature  at  the  frost 
surface.  To  facilitate  the  analysis,  a  simplified  model  of  the  frost 
layer  has  been  considered.  Deviation  of  the  theory  from  the  actual 
frost  formation  is,  therefore,  expected  and  the  results  remain  to  be 
validated. 

In  reality  the  frost  properties  such  as  density  and  thermal  con¬ 
ductivity  are  strongly  dependent  on  time.  A  recent  investigation  in 
this  connection  has  been  reported  by  Yonko  and  Sepsy25.  A  re -examina¬ 
tion  of  the  analysis  discussed  in  this  chapter  indicates  that  the  method 
also  applies  to  the  cue  of  variable  frost  properties  provided  the 
thermal  conductivity  and  density  of  the  frost  are  functions  of  a  (a 
being  two  times  the  square  root  of  the  product  of  frost  thermal  diffu- 
sivity  and  time).  For  this  case  proper  care  has  to  be  exercised  in 
carrying  out  the  differentiations  with  respect  to  a,  for  k,  Ci,  C2,  C3 
and  C«  are  also  functions  of  a. 
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CHAPTER  III 


THE  ELECTRIC  WIND 


3.1  GENERAL  NATURE  OF  THE  ELECTRIC  WIND 

The  phenomenon  of  electric  wind  (also  called  corona  wind)  is  a 
characteristic  feature  of  asymmetric  field  corona  discharge  such  as 
point-plane  or  wire -plane  corona.  Although  the  first  theory  of  elec¬ 
tric  wind  was  given  by  Chat  toe  k20  as  early  as  l899>  very  little  work 
has  been  reported  in  the  literature  ever  since.  As  a  result  the  actual 
flow  processes  of  the  induced  air  stream  is  still  not  well  understood. 
Of  the  scant  work  available,  all  authors  discussed  essentially  only 
one-dimensional  cases.  Typical  examples  are  plane -parallel,  coaxial- 
cylindrical  and  concentric -spherical  electrode  configurations.  An 
exhaustive  historical  review  and  bibliography  has  recently  been  pre¬ 
sented  by  Robinson27  and  no  further  literature  survey  will  be  attempted 
here.  In  this  chapter  the  analysis  will  only  be  developed  to  the  extent 
as  is  required  for  proper  interpretation  of  the  experimental  data  to  be 
presented  in  Chapter  VI. 

Before  proceeding  to  the  discussion  of  the  electric  wind  generated 
by  specific  electrode  configurations,  the  general  nature  of  the  problem 
in  question  will  be  examined  first.  It  appears  that  one  of  the  main 
reasons  why  the  phenomenon  of  corona  wind  has  received  rather  inactive 
attention  is  that  the  phenomenon  is  too  involved  to  permit  analytical 
solutions  for  physically  realistic  models.  To  help  illustrate  the  com¬ 
plexity  of  the  problem  the  governing  equations  are  listed  as  follows: 

Electrodynamic  equations: 


Hydrodynamic  equations: 
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It  is  assumed  here  that  the  flow  is  incompressible  and  that  the  Coulomb’s 
force  is  the  only  extraneous  force.  In  Eq.  (3.4)  the  charge  flux  due 
to  diffusion  is  neglected. 

An  understanding  of  the  actual  aerodynamic  processes  of  the  electric 
wind  necessarily  demands  a  solution  of  the  equations  of  motion  (3 *6), 
which  in  turn  requires  and  expression  for  the  electrostatic  body  force 
(that  is  the  Coulomb’s  force).  The  task  of  solving  the  non-linear 
Navier-Stokes  equations  even  in  the  absence  of  extraneous  forces  is 
already  a  difficult  one;  not  to  mention  that  the  determination  of  the 
Coulomb's  force  is  not  easy  as  well.  In  fact  it  is  evident  from  Eqs. 

Eqs.  (3*1)  to  (3*6)  that  the  electrodynamic  and  hydrodynamic  equations 
are  coupled  when  the  velocity  of  the  fluid  is  comparable  to  the  drift 
velocity  of  the  space  charges.  In  ga§es  it  is  generally  considered 
that  v  is  negligible  as  compared  to  KE,  resulting  in  decoupling  of  the 
two  sets  of  equations.  This,  however,  does  not  resolve  the  difficulties 
involved  in  the  problem,  for  even  then  the  Poisson  equation  (3.1)  can 
only  be  solved  for  a  few  special  cases.  The  common  feature  of  these 
special  cases  is  that  the  geometrical  configuration  of  the  electrodes 
should  be  such  that  the  charge  density  pe  could  be  expressed,  with  the 
aid  of  Eqs.  (3.3)  and  (3.4),  in  terms  of  the  electric  field  intensity, 
the  appropriate  spatial  coordinate,  and  an  experimentally  determinable 
quantity  (usually  the  total  current  at  the  cathode).  The  electrode 
configurations  possessing  this  property  include  parallel  planes,  coaxial 
cylinders  and  concentric  spheres. 

For  the  actual  physical  situation  of  wire-plane  or  point-plane 
corona,  the  electric  field  intensity,  the  charge  density,  and  hence 
the  current  density  are  strong  functions  of  the  relevant  spatial  coor¬ 
dinates.  Lack  of  an  explicit  relation  between  pe  and  E,  and  hence  the 
Coulomb's  force  can  not  be  determined  analytically.  Experimentally 
only  the  current  density  distribution  along  the  plane  electrode  has 
been  measured.  The  data  obtained  by  Chattock26  for  point-plane  corona 
and  those  to  be  presented  in  Chapter  VI  for  wire-plane  corona  exhibit 
a  bell-shaped  distribution  of  current  density,  when  the  latter  is 
plotted  against  the  distance  frem  the  center  or  centerline  of  the  plane 
which  is  taken  to  be  the  point  or  line  of  normal  projection  of  the 
point  or  wire  electrode.  This  Implies  that  the  electric  wind  is  con¬ 
fined  to  a  very  narrow  region  around  the  point  or  wire,  and  its  velo¬ 
city  distribution  is  strongly  non-uniform. 

The  non-uniformity  of  the  jet  of  electric  wind  suggests  that  the 
flow  generated  by  the  non-uniform  electric  field  is  rotational.  It  can 
be  shown  that  the  rotationality  might  be  a  direct  consequence  of  the 
non-conservative  electrostatic  force  and  needs  not  necessarily  be  a 
result  of  viscous  effects.  If  one  takes  the  curl  of  Eq.  (3.6)  and  goes 
through  some  mathematical  manipulations,  there  results  an  expression 
for  the  rate  of  change  of  vorticity 


where 
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Is  the  substantial  derivative.  From  this  equation  it  is  seen  that  in 
the  absence  of  viscous  effects  an  electric  wind  will  remain  irrotational 
(if  it  is  initially  so)  only  either  when  the  charge  density  is  constant 
or  when  the  electric  field  is  colinear  with  the  gradient  of  the  charge 
density.  Neither  condition  is  satisfied  in  a  wire-plane  or  point-plane 
corona  and  the  electric  wind  thereby  induced  is  rotational. 

To  gain  a  further  insight  into  the  nature  of  an  electric  wind,  it 
is  helpful  to  compare  it  with  a  free  jet.  Despite  the  fact  that  the 
two  flows  are  governed  by  different  equations  of  notion  (the  difference 
being  in  the  extraneous  force  term) ,  a  parallel  between  them  may  be 
drawn  in  certain  aspects.  An  electric  wind  is  analogous  to  a  free  jet 
in  the  sense  that  both  discharge  into  a  still  fluid  of  the  sane  kind. 
While  the  total  momentum  of  a  free  jet  in  the  direction  of  jet  axis 
remains  constant  and  Independent  of  distance  from  the  orifice,  the 
quantity  which  possesses  similar  chcracteijistiqsfor  an  electric  wind 
is  the  total  force  (considering  the  cause  v  «  Ke).  Theoretically  a 
free  jet  can  be  regarded  as  emerging  frcm  a  line  source  (the  discussion 
having  been  restricted  to  two-dimensional  case,  similar  statements  apply 
to  axisynmetrlc  situation)  which  issues,  in  one  direction  only,  a  finite 
volume  of  flow  with  a  finite  momentum  (28,  29).  In  the  case  of  electric 
wind,  the  air  stream  is  generated  by  a  finite  force 

F  -  -  /  Jdx  ■  -  (x  being  transverse  to  the  direction  of  the  wind) 
k  k 

at  a  very  thin  region  near  the  emitting  corona  wire  facing  the  collect¬ 
ing  electrode.  It  appears  that  an  electric  wind  could,  therefore,  be 
considered  as  emerging  from  a  line  source  in  a  direction  normal  to  the 
collector  through  the  action  of  a  finite  force  F.  However,  a  discre¬ 
pancy  exists  between  the  two  types  of  line  sources;  whereas  in  a  free 
jet  the  discharging  fluid  is  Introduced  externally,  the  streaming  air 
in  an  electric  wind  must  come  from  the  surroundings  (presumably  frcm 
the  region  behind  the  wire).  Bnerglng  frcm  a  slit  with  uniform  velo¬ 
city,  the  free  jet  carries  with  it  some  of  the  surrounding  fluid  origi¬ 
nally  at  rest  because  of  free  mixing  developed  on  its  periphery  -8  a 
result  of  friction,  an  action  usually  referred  to  as  entrainment.  The 
electric  wind  induced  by  a  wire-plane  corona  is  inherently  non-uniform; 
the  velocity  vanishes  at  Its  periphery  mainly  because  the  space  charge 
density  there  is  zero.  It  is  thus  speculated  that  entrainment  at  the 
boundaries  of  an  electric  wind  is  probably  not  as  significant  as  in  a 
free  jet. 


3.2  ONE-DIMENSIONAL  MODEL 

The  case  considered  first  is  a  one -dimensional  electric  wind  gener¬ 
ated  by  parallel  plane  electrodes  shown  in  Fig.  4.  The  emitting  elec¬ 
trode  (the  anode)  is  considered  to  consist  of  closely  spaced  parallel 
corona  wires  located  at  x=Oj  while  the  collecting  electrode  (the  cathode) 
is  in  the  form  of  a  plane  wire  screen  (at  x=x2)  which  allows  the  induced 
wind  to  pass  through.  A  voltage  of  V  (kv)  impressed  across  the  elec¬ 
trodes  caises  a  current  I  (amp)  to  flow,  producing  a  current  density 
J  -  i/S  (amp/m2).  Since  the  electrodes  can  not  sustain  any  pressure 
gradient,  the  atmospheric  pressure  prevails  throughout  the  entire  space 
between  the  electrodes,  and  the  action  of  the  electrostatic  body  force 
on  the  air  molecules  results  in  the  st .-earning  of  the  electric  wind. 

For  one -dimensional  approximation,  J  E  and  v  have  only  a  component 
in  the  x  direction  and  are  functions  of  x  alone.  In  the  absence  of 
solid  boundaries  and  disregarding  the  free  mixing  region  along  the  edges 
of  the  electrodes,  viscous  effects  on  the  flow  can  be  neglected.  The 
governing  equations  thus  reduce  to,  for  steady  flow, 
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£  5=  -  pe 

dx 


dJ 


dx 


J  ■  pe  (u+KE ) 


du  „ 
pu  ^  -  OeE 


Combination  of  Eqs.  (3*8)  and  (3. 11)  gives 


du  _  dE 
pu  ~  ■  €  E  ^ 
dx  dx 


which  can  be  solved  to  yield 

pu2  -  e  E2  =  constant  =  Ki 


(3.8) 

(3.9) 
(3.10) 
(3.U) 


(3.12) 


From  Eq.  (3.9),  it  is  apparent  that  the  current  density  is  constant. 
Using  Eqs.  (3.10)  and  (3.12)  in  Eq.  (3.8),  there  is  obtained  an  equation 
determining  the  electric  field  intensity  E. 
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€  [jf n/e2  +  Ei/ €  +  Ke]  g  =  J 


Upon  integration,  it  is  found 

|e  n/E2  +  (Ki/e)  +  (Ki/e)ln  [e  +  n/e2  +  (Ki/e)  ]  J  +  ^  ■  Jx  +  K2 

(3.13) 

where  Ks  is  a  constant  of  integration. 

It  remains  to  determine  the  two  constants  Ki  and  K2.  The  fact 
that  the  highly  ionized  region  constituting  the  corona  envelope  occupies 
only  a  very  small  volume  (say  between  x  «  0  and  x  *  xx)  near  the  emitter 
provides  the  appropriate  boundary  conditions  at  the  emitter.  In  this 
region  the  positive  and  negative  ions  are  approximately  equal  in  number. 
Therefore,  it  may  be  assumed  that  pe  *  0  there  and  hence  the  absence  of 
electric  force.  At  x  =  xx  the  field  strength  is  equal  to  the  spark-break- 
down  field  strength  Ei.  Disregarding  the  entrainment  of  air  molecules 
from  the  surroundings  (x  <  0)  the  velocity  of  the  air  atx*xx  may  be 
assumed  zero.  There  conditions  permit  the  determination  of  Ki  and  fe. 
They  are 

Ki  -  -  €  Ei2 


and 


K2  -  |  Ei2(K  -  InEi)  -  Jxi 

Thus  Eqs.  (3.12)  and  (3.13)  becomes 

pu2  -  c  (E2-Ex2)  (3.14) 

and 

€  {Hz  £e  >/e2  -  El2  -  Ei2ln(E  We2  -  Ex)]  +  |  E2  j 

-  J(X-xx)  +  (k  -  'JTTp  InEi)  (3.15) 

For  the  sake  of  simplicity,  it  will  be  assumed  that  Ex  is  zero  and  xx 
is  negligible  against  xg.  With  these  assumptions  the  velocity  of  the 
induced  electric  wind  at  the  screen  electrode  is  found  to  be 
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02  ■  *^JX2/p(^€7p+K)  (3.16) 


and  the  velocity  at  x  is 


u  *  >/2J(x-xi)/p(>/e7p+K)  (3.17) 

That  the  assumption  of  zero  field  at  emitter  is  far  from  true  has  been 
emphasized  by  Stuetzer30.  Fortunately  a  similar  expression  can  be 
derived  in  an  alternative  manner.  When  u  is  negligible  as  compared  to 
KE,  Eq.  (3*H)  can  be  written 


(3.18) 


which,  upon  integration,  yields 


u  =  V2J(x-xx)/pK  .  >/2l(x-x1)/pKS 


(3.19) 


At  x=xs,  the  velocity  is  approximately 

us  -  >/2Jxs/pK  =  V2IX2/pKS  (3.20) 

A  great  disadvantage  of  the  one -dimensional  model  of  electric  wind 
is  that  the  continuity  equation  is  not  satisfied.  It  is  therefore  con¬ 
cluded  that  continuous  streaming  of  the  electric  wind  requires  replenish¬ 
ment  from  the  surroundings  by  entrainment. 


3.3  TWO-DIMENSIONAL  MODEL 

The  electric  field  associated  with  wire-plane  electrode  configura¬ 
tion  is  two  dimensional  and  the  electric  wind  thereby  induced  is  two- 
dimensional  as  well.  Referring  to  Fig.  4,  instead  of  a  plane  of  parallel 
wires,  it  is  now  considered  that  the  emitting  electrode  is  a  single 
wire.  The  physical  model  and  the  coordinate  system  cure  depicted  in 
Fig.  5.  The  asymmetry  of  the  electric  field  and  the  absence  of  a  theory 
of  breakdown  in  non-uniform  fields  make  it  not  possible  to  determine 
analytically  the  distribution  of  the  electric  wind,  certain  assumptions 
have  to  be  made: 

1.  The  transverse  velocity  component  u  of  the  jet  of  electric 
wind  i 8  small  compared  to  the  longitudinal  velocity  com¬ 
ponent  v. 
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Fig.  5  Two-dimensional  model  of  electric  wind 


2.  ^x/cbc  is  small  compared  to  80  that  the 

Poisson*!  equation  (3*l)  may  be  approximated  by 


c 


3.  For  a  first  approximation,  the  value  of  Ey  in  the 

absence  of  space  charges  will  be  used  to  replace 
the  actual  Ey. 

On  the  basis  of  the  first  two  assumptions,  the  equations  of  motion 
reduce  to  a  single  equation 


pv  |  .  e  Ey  (3.2!) 

which  upon  integration  becomes 

pv2  -  e  Ey2  -  f0(x)  (3.22) 

where  fo(x),  a  function  of  x  only,  is  introduced  through  Integration. 

Once  again  the  wind  will  be  considered  to  start  blowing  from  at 
rest  near  the  wire  so  that  the  boundary  condition  which  must  be  satis¬ 
fied  by  Eq.  (3.22)  is,  assuming  the  radius  of  the  wire  is  small  against 
d. 


u  =  0  when  [  Ey  ®  Ey]y^  »  Ey(x,d) 

Using  this  the  function  fo(x)  is  found  to  be 

f0(x)  -  -  e  [Ey(x,d)]2 
and  Eq.  (3*22)  can  be  rewritten 

PV2  «  €  -  [ Ey (x,d)]  2|  (3.23) 

Assuming  again  the  field  vanished  at  the  edge  of  ionization  zone  near 
the  emitting  electrode,  Eq.  (3.23)  takes  a  simple  form 


v  « 'fe/p  Ey 


(3.24) 


An  analogous  expression  relating  the  pressure  to  the  field  strength  ha 
been  deduced  by  Velkoff*  by  means  of  a  different  argument.  Neglecting 
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space  charge  distortion  the  electric  field,  as  noted  by  Velkoff ,  can 
be  expressed  as 


%  = 

In  r— 


(y-d): 


(y^)J 


when  the  radius  ra  of  the  wire  electrode  is  small  as  compared  to  the 
distance  d  of  the  wire  from  the  plane  electrode.  Following  assumption 
3,  Eq.  (3 -24)  becomes 


y  -  d 

c2  +  (y-d)2 


y  + d 

+  (j*Y 


(3.25) 


As  in  one -dimensional  case,  the  assumption  of  negligible  field  near  the 
emitting  electrode  is  unrealistic.  It  has  been  pointed  out  by  Kapzow31 
that  at  the  emitter  the  field  is  the  highest  in  the  system.  The  valid¬ 
ity  of  Eq.  (3.25)  thus  needs  further  justification. 

A  different  approach  makes  use  of  the  relation  PgEy  =>  Jy/K  and 
assumption  1.  the  other  two  assumptions  being  discarded.  Thus  in  place 
of  Eq.  (3.21),  an  expression  similar  to  Eq.  (3.18)  is  obtained 


Jy/K 


(3.26' 


Integration  from  the  emitter  y  =  ye,  where  v  =  0,  to  the  collector  y  =  yc 
gives  the  following  relation  for  the  velocity  distribution  at  the  col¬ 
lecting  electrode 


1  iCH" 


(3.27) 


This  equation  is  similar  to  the  one  obtained  by  Robinson27,  who  further 
considered  that  Jy  is  proportional  to  the  total  current  I  and  arrived 
at  a  relation  between  vc  and  I,  for  coaxial  cylindrical  and  concentric 
spherical  electrode  configuration,  namely, 


VolC/ 


(3.28) 


where  K'  is  a  constant  of  proportionality  to  be  determined  experimen¬ 
tally.  Since  two-dimensional  situation  is  being  considered  here,  vc 
is  a  function  of  x,  it  follows  that  K'  must  also  depends  on  x  for 
Eq.  (3.28)  to  be  valid. 
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3.4  INVISCID  ELECTRIC  WIND  WITH  STAGNATION 


The  generation  of  the  electric  wind  is  usually  associated  with  a 
solid  plane  which  acts  as  the  collecting  electrode.  The  presence  of 
this  solid  plane  tends  to  make  the  situation  more  intractable.  When 
the  wind  discharges  onto  the  plane,  a  static  pressure  gradient  is 
expected  to  build  up  near  the  stagnation  point.  For  the  electrode  con¬ 
figuration  of  a  wire-plane  corona,  the  electric  field  component  E* 
parallel  to  the  plane  is  in  practice  negligible  as  compared  to  the  com¬ 
ponent  Ey  normal  to  the  plane  near  the  plane  surface.  Therefore,  the 
equations  of  motion  for  a  two-dimensional  invlscid  flow  in  the  neighbor¬ 
hood  of  a  stagnation  point  can  be  written  as 


pu 


du 

£ 


(3.29) 


&  %  ay 


+  PeEy 


(3.30) 


Instead  of  solving  these  two  equations,  which  is  no  easy  task,  a 
simple  relation  between  the  velocity,  pressure  and  electric  quantity 
will  be  derived.  When  multiplied  by  dx  and  dy  respectively  and  use  is 
made  of  the  relation  along  a  stream  line 


dx  .  dy 
u  v 


Eqs.  (3.29)  and  (3*30)  combine  to  give  an  equation  of  the  form 

O^+^J  +  dp  -  PgEydy  -  0  (3.31) 

Upon  integration,  there  results  a  Bernoulli  equation  for  the  problem 
of  interest 

|  (u2-^)  +  p  -L  1  PeEydy  ■  constant  (3.32) 


where  s.l.  indicates  the  integration  is  along  a  streamline.  In  general 
the  constant,  the  so-called  Bernoulli  constant,  takes  different  values 
for  different  streamlines  when  the  flow  is  rotational.  When 


Eq/  (3.32)  can  be  written 


9 
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^  (u2+\r)  +  p  -  \  e  Ey2  ■  constant.  (3*33) 

If  on  the  other  hand  the  relation  peEy  =  Jy/K  is  used,  it  follows 


f  (U2V)  +  p  -  [  5l  dy 


=  constant. 


(3.34) 


The  x  component  of  the  pressure  gradient  is  of  utmost  importance  in 
considering  viscous  stagnation  flow.  Differentiation  of  Eq.  (3*34) 
with  respect  to  x  gives 


(Bernoulli  constant) 


(3.35) 

It  is  noted  from  this  equation  that  dp/dx  is  determined  by  the  velocity 
field  as  well  as  the  electric  field. 


For  one  -dimens  tonal  electric  wind  Jy  =  constant  and  the  Bernoulli 
constant  takes  an  identical  value  throughout  the  flow  field.  It  follows 
that  for  this  particular  case 


(3.36) 


The  relation  (3.36)  also  holds  true  for  two-dimensional  potential  flows 
in  the  absence  of  extraneous  forces.  Hence  the  flow  produced  by  a  one¬ 
dimensional  electric  wind  Impinging  on  a  plane  electrode  can  be  treated 
as  a  classical  two-dimensional  stagnation  flow.  For  the  latter  case 
the  velocity  distribution  in  frictionless  potential  flow  in  the  neigh¬ 
borhood  of  the  stagnation  point  is  given  by20 

u  -  uix;  v  »  -  uiy  (3.37) 

where  ux  is  a  constant  determined  by  the  oncoming  stream.  Thus 

^  ■  pu  7  =  ■  puj/Sc  (3.38) 

oc  ox 


CHAPTER  IV 


HEAT  AND  MASS  TRANSFER  IN  TWO-DIMENSIONAL  STAGNATION  FLCW 


4.1  GENERAL  CONSIDERATIONS 

In  view  of  the  non-uniformity  of  the  impinging  tvo -dimensional  jet 
of  electric  wind,  the  flow  field  in  the  neighborhood  of  the  stagnation 
point  is  rather  complex.  At  the  present  state  of  our  knowledge,  it  is 
impossible  to  obtain  an  exact  solution  for  the  actual  stagnation  flow 
resulted  from  the  impingement  of  a  non-uniform  stream  with  arbitrary 
velocity  distribution.  Information  on  the  heat  and  mass  transfer  pro¬ 
cesses  at  the  stagnation  point,  however,  may  be  secured  by  considera¬ 
tion  of  the  stagnation  flow  under  an  impinging  uniform  stream. 

Consider  a  uniform  electric  wind  induced  by  positive  corona  dis¬ 
charges  taking  place  between  a  plane  of  closely-spaced  parallel  wires 
and  a  solid  plane  wall  electrode.  Moving  in  the  direction  of  the 
electric  field,  the  wind  arrives  from  the  y-axis,  strikes  normally  on 
the  plane  electrode  placed  at  yO,  divides  into  two  streams  on  the  wall 
and  leaves  in  both  directions,  Fig.  6.  The  temperature  and  the  mass 
fraction  of  water  vapor  of  the  impinging  air  are  considered  to  be  the 
same  as  the  surrounding  air.  The  facts  that  the  streaming  air,  though 
induced  electrically,  is  actually  replenished  from  the  surroundings 
through  entrainment,  and  that  the  corona  discharge  does  not  have  signi¬ 
ficant  effects  on  the  temperature  and  moisture  content  of  the  air, 
seem  to  lend  strong  support  to  this  assumption. 

The  viscous  flow  theory  for  two-dimensional  stagnation  flow  has 
already  been  well  developed.  It  has  long  been  recognized  that  for  this 
type  of  flow  the  solution  based  on  Prandtl's  boundary  layer  approxima¬ 
tion  is  also  an  exact  solution  to  the  complete  Navler-Stokes  equations 
for  two  -dimen  s  ional  flows.  For  rapid  calculations  of  the  heat  and  mass 
transfer  values  near  the  stagnation  point,  it  is  proposed  to  make  use 
of  the  integral  method  based  on  the  boundary  layer  equations. 

The  complete  set  of  the  differential  equations  which  describe  the 
combined  momentum,  energy  and  mass  transport  in  a  two-dimensional  incom¬ 
pressible  steady  laminar  boundary  layer  of  a  binary  mixture,  which  is 
humid  air  for  this  case,  are,  as  given  by  Eckert  and  Drake32, 

Continuity: 

!♦!-»  <m> 

Momentum: 

u2“  +  v2<i.uSK  +  vSfa  (i,.2) 

dx  dy  dx  dy2 


4l 


Energy: 


Diffusion : 


a 


a2? 


(4.3) 


(4.4) 


The  pressure  gradient  term  in  the  momentum  equation  has  been  substituted 
in  favor  of  the  free  stream  velocity.  In  the  energy  equation,  the 
viscous  dissipation,  the  thermal  diffusion  and  the  diffusion  thermo- 
effect  are  neglected.  The  diffusion  equation  is  written  in  teims  of 
the  mass  fraction  wi  of  the  water  vapor. 

The  boundary  conditions  appropriate  to  the  flow  over  a  frosting 
surface  are 

At  y  =  0  :  U“0  ,  v=v^  ,  T  ®  T^  ,  wi«  wij.  (4.5) 

At  y  =  00  :  u  ■  U  ,  T  ■  T»  ,  vi  =  wi^  (4.6) 


The  mass  fraction  wi-p  and  temperature  Tf  at  the  frost  surface  will  be 
assumed  constant.  The  condition  v  =  vf  accounts  for  the  fact  that  a 
convective  flow  is  generally  connected  with  mass  transfer  from  or  to  a 
surface.  Such  a  convective  flow  is  always  present  when  the  surface  is 
impermeable  for  one  component  of  the  mixture.  In  order  to  complete 
the  statement  of  the  problem  under  consideration,  the  normal  velocity 
at  the  frost  surface  Vf  is  to  be  supplied  by  application  of  the  con¬ 
dition  that  the  mass  transfer  of  air  to  the  frost  surface  is  considered 
to  be  zero24.  The  mass  flux  of  air  consists  of  two  parts:  a  diffusive 
mass  flux  and  a  convective  flux.  In  the  absence  of  diffusion  thermo- 
effect  and  thermal  diffusion,  the  total  mass  flux  of  air  at  the  frost 
surface  is 

dwp 

-pD  +  PaVf  "  0 

where  the  subscript  2  denotes  the  second  species  of  the  mixture,  which 
is  dry  air  in  this  case.  The  mass  fraction  of  the  air  can  be  expressed 
in  terms  of  the  mass  fraction  of  the  water  vapor  by 


■  W2  ,  W2  a  1  —  wi  ,  _  _  dwi 

p  dy  ’  dy 

This  gives 

D  ^ 

Vf*  =  -  — - - 

1  -  wi  dy 


(4.7) 
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where  w*  and  its  derivative  are  to  be  evaluated  at  the  frost  surface. 


The  integral  relations  which  form  the  basis  of  the  integral  method 
are  most  simply  obtained  by  integration  from  the  boundary  layer  equations 
(4.1)  through  (4.4).  Equations  (4.1)  and  (4.2)  may  be  written  as 

-v^i-U^-u^-v^ 

Sy2  dx  &  2y 

o  -  (u-u)  +  (u-u)  ^ 

«  cy 

where  Eq.  (4.1)  has  been  multiplied  by  (U-u).  By  addition  of  these, 

•  ^~m  4:  (Uu"u2)  +  (u“u)  —  +  (vu-vu) 

3y*  «  dx  ay 

Upon  Integration  with  respect  to  y  from  0  to  •  this  yields,  since 
au/ay  and  v(U-u)  tend  to  zero  as  y  -» •,  the  equation  expressing  the 
balance  of  momentum 

v|y-)'s/o  (U-u)uay  +  I  fQ  (u-u)dy  *  vfu 

(4.8) 

The  same  process  of  integration  applies  to  Eqs.  (4.3)  and  (4.4) 
leads  to  similar  equations  expressing  the  balance  of  heat  and  material. 

In  terms  of  the  new  variables  9  and  q>  defined  by 


e  ,  I-r-JBC.  <p  ,  Vl.7  v^. 

T.  -  Tf  wiw  -  wij 

the  resulting  integral  equations  have  the  following  forms 

£  (!)  y*>  ■  5  r  '9)uay  ■  Vf  (M) 

£  (f)  iro  '  X  J  (1'<p)udjr  •  rt  {kM) 

VQ 

where  Pr  *  v/a  is  the  Prandtl  number  and  Sc  ■  v/d  is  the  Schmidt  number. 


4.2  SOLUTION  FOR  Vf  NEGLIGIBLY  SMALL 

The  simpler  cue  where  vf  is  negligible  will  be  considered  first. 
The  problem  is  simplified  considerably  u  a  result  of  the  decoupling  of 
the  momentim  equation  from  the  diffusion  equation  when  the  boundary 
condition  v  •=  Vf  ■  0  is  used.  Thus  the  momentum  integral  equation  may 
be  solved  independently  in  exactly  the  same  manner  as  Bohlhausen's 
approximate  method.  An  excellent  polynomial  representation  of  the 
velocity  profile  takes  the  fallowing  font 

where  U  *  uix,  A  ■  B^ui/v,  &  is  the  boundary  layer  thickness  and  ui'  is 
a  proportionality  constant  having  the  dimension  sec*1.  It  has  been 
shown  in  Schlichting2®  that  the  profile  for  stagnation  point  flow 
corresponds  to  X  -  7*052,  so  that 

8  -  (7.052  (4.13) 

'  Me  ' 


In  order  to  solve  the  energy  and  the  diffusion  integral  equations, 
suitable  functions  have  to  be  chosen  for  9  and  9.  On  account  of  the 
similarity  in  form  of  Eqs.  (4.9)  and  (4.10)  in  the  absence  of  Vf  and 
the  identity  of  the  boundary  conditions  (4.5)  and  (4.6)  for  0  and  9,  it 
is  only  necessary  to  carry  through  the  calculation  for  one  of  the  two 
equations.  The  solution  thus  obtained  is  also  a  solution  of  the  other 
equation  provided  Pr  ■  Sc.  If  Pr  j*  Sc,  it  is  necessary  to  replace  Pr 
by  Sc  or  vice  versa,  as  cue  may  be,  in  the  process  of  calculation. 

As  the  temperature  profile,  the  function 


(4.14) 


which  satisfies  the  following  boundary  conditions 


y  -  0  , 

y  -  Bt  , 


0-0  , 

0-1  , 


£°  -o 
—  -  0 


(4.15) 


will  be  used.  Although  strictly  the  conditions  at  Infinity  are  only 
approached  asymptotically,  it  is  assumed  that  these  conditions  can  be 
transferred  from  Infinity  to  y  -  ftj  without  appreciable  error  (siailar 
assumption  applies  to  the  velocity  and  material  boundary  layers).  As 
seen  from  Eq.  (4.3),  the  conditions  (4.15)  are  all  satisfied  by  the 
exact  solution.  Here  the  thermal  boundary  layer  thickness  ft^  is  inde¬ 
pendent  of  x,  since  Tf  is  assumed  to  be  spatially  uniform". 

Introducing  Eqs.  (4.12)  and  (4.14)  in  Eq.  (4.9)  and  carrying  out 
the  indicated  integrations  and  differentiations ,  Vf  being  treated  zero, 
there  results 


This  expression  is  only  strictly  valid  of  ftr  £  ft,  which  is  not  usually 
the  case,  since  Eq.  (4.12)  is  only  valid  for  y  £  ft.  However,  the 
algebra  would  become  more  cumbersome  with  little  gain  in  accuracy  if  the 
strict  expression  were  used. 

Writing  dp  -  — ,  Eq.  (4.l£)  becomes 
ft 

-  14.9315  V  +  67.0067  £%  -  144.991)  -  -  97-1503  (^) 

(4.17) 

An  approximation  to  the  solution  of  this  equation  when  Pr  does  not  differ 
too  much  from  unity,  i.e.  dp  1,  can  be  obtained  by  setting  dp  ■  1  in 
the  parenthesis.  Hence,  approximately, 

Ap  -  (Pr)"1/5  (4. IS) 

An  accurate  solution  of  Eq.  (4.17)  for  Pr  has  been  calculated  by 
Sparrow88.  Hs  gave  the  value  dp  -  1.177*  Calculation  based  on  Eq.  (4.18) 
gives  dp  -  1.14.  Thus  the  expression  (4.10)  is  in  error  by  less  than 
5%  as  compared  with  the  accurate  solution. 


♦  r  >  « 


Similarly,  the  material  profile  may  he  assumed  to  hare  the  form 


1  JL  -  i  (JLY 

2  8*  2  \5»/ 


(4.19) 


which  satisfies  the  boundary  conditions 


7  -  0  ,  flP  - 


7  ■  8m  •  <P  -  1 


I- 


(4.20) 


where  8*  is  the  material  boundary  layer  thickness. 

Writing  4i  ■  &•/&(  Eq.  (4.17)  msy  now  be  modified  to  read 

Ama(dma  -  14.9315  Am2  ♦  67.0067  A*  -  144.991)  -  -  97.1503/Sc 

(4.21) 

The  approximate  solution  corresponding  to  Eq.  (4.13)  now  beccsms 

(Sc)"1/3  (4.2 

The  heat  flow  to  the  plate  near  the  stagnation  point  per  unit  area  is 
given  by.  In  the  absence  of  mess  transfer, 


-*«) 


dr\ 


Upon  substitution  from  Eq.  (4.14),  (4.10)  end  (4.13),  it  is  obtained 


«  ■  ik(7^f(Pr)l/‘(I'Tf) 


(4.23 


On  the  other  hand  a  heat  transfer  coefficient  h  can  be  introduced 
through  "Eewton's  law  of  cooling" 


» 


Equating  this  .aist  two  expressions  gives 

h  «  2k/_Sl_\,A(Pr)x/a  (4.24) 

2  \7.052v  /  '  ' 

The  mass  flux  to  the  plate  near  the  stagnation  point  can  be  obtained 
in  a  similar  fashion,  namely. 

» -  -  (*••«) 

And  the  mass  transfer  coefficient  defined  as 


m  -  -  P^(wi„-wif) 

can  be  readily  deduced  from  Eq.  (4.24) 


- 1 D  (?^)  *  (sc>v“  (4-26) 


4.3  SOLUTION  FOR  vf  NOT  ZERO 

To  account  for  the  connection  of  the  convective  flow  at  the  surface 
with  the  ovmr-all  transport  process,  it  is  required  to  deal  with  the 
coupled  system  of  integral  equations  (4.8)  through  (4.10)  simultaneously. 
The  method  which  will  be  used  is  essentially  similar  to  that  of  the 
preceding  section.  The  only  differences  are  that  in  selecting  the  suit¬ 
able  forms  of  the  velocity,  temperature  and  material  profiles  the  bound¬ 
ary  condition  that  v  ■  Vf  has  to  be  taken  into  account  and  that  terms 
involving  vf  are  retained  in  the  integral  equations.  Eckert  and 
Llebleln*4  first  applied  the  method  to  solve  Eqs.  (4.8)  and  (4.10)  for 
the  particular  case  of  the  vaporization  of  water  from  a  flat  plate  in 
forced  convection.  An  improved  version  of  the  method  was  used  by 
Spalding"  to  find  solutions  of  the  same  two  equations  for  three  types 
of  flows:  forced  convection  from  a  flat  plate,  free  convection  from  a 
vertical  plate  and  the  forward  stagnation  point  of  a  sphere  suspended 
in  a  fluid  stream.  All  three  equations  were  solved  in  connection  with 
the  frosting  plate  in  forced  convection  by  Sugawara  et  al3®. 

The  integral  equations  (4.8)  through  (4.10)  may  be  reduced  to 
simplified  forms  by  introduction  of  new  quantities  defined  in  terms  of 
the  Integrals  appearing  in  the  equations .  Considering  Eq.  (4.8)  first, 
it  is  a  well-known  fact  that,  when  expressed  in  terms  of  the  displace¬ 
ment  and  momentum  thicknesses,  which  are  defined,  respectively, 
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a*  ■  y6  (a.  -  g  ay)  (4.27) 

ud 

*•*  *  C  ■  §  *)  (i,-28) 

the  moment  im  integral  equation  reduces  to 

(|)^ + Tf°  •  5 +  **  ui  c*-^) 


Similarly,  Bq«.  (4.9)  and  (4*10)  can  be  rewritten  as 

s(|U+^-s(°T)  <k-*» 

and 


♦  Tf  .  A(wj») 


(4.31) 


where  &*»  and  5**  are  defined  by 

T  a 

e*».jp[i-e)a<ijr  (k.32) 


and 


•  MM 

Qwsr 


-/"(l-qOgdy 


(4.33) 


Following  Spalding,  it  will  be  assamd  that  the  Telocity,  tempera¬ 
ture  and  Material  profiles  can  be  expressed  in  cubic  polynomial  font, 
nsMly, 


ir  * ,w  *>? ♦  ♦  ‘K#)’  (4-34) 


8  -  0(y)  -  bxX  ♦  b2(i)*  ♦  b.(X)»  (4.35) 

and 

9  -  H(y)  -  cij  ♦  ca^)2  +  ca(^)a  (4.36) 


49 


f 


The  coefficients  of  these  functions  may  be  determined  by  consideration 
of  the  boundary  conditions,  which  are 


u  ■  0  , 

(4.37) 

0*0  , 

de  v 

f  v  Pr  v 

(4.38) 

9 . 0  , 

dy  Sc  cfcr2 

(4.39) 

y  *  B  : 

u  *  U  , 

o 

■ 

(4.40) 

y  ■  * 

#  -1  , 

o 

■ 

*1^ 

(4.41) 

y  -  Cm  : 

<p  -  1  , 

■ 

o 

(4.42) 

The  second  conditions  of  (4.37)  through  (4.39),  which  are  direct  con¬ 
sequences  of  the  differential  equations  (4.2)  through  (4.1),  are  of 
particular  Importance.  They  determine  the  curvatures  of  the  velocity, 
temperature  and  material  profiles  near  the  surface.  It  is  through  them 
that  the  shapes  of  the  profiles  are  affected  by  the  convective  flow  at 
the  surface,  thus  resulting  in  changes  in  the  over-all  transport  process. 

Inserting  the  value  of  Vf  from  Eq.  (4.11),  there  is  dbtained  the 
following  relations  determining  the  coefficients  of  the  functions  (4.34) 
through  (4.36) 

A  ■  4ai  -  6  -  B"ai2  , 

•a  ■  3  -  2ai  ,  a3  -  ax  -  2  , 

[l  -  "'I  -  <3/2)B']  ,  (4.43) 

b2  ■  3  -  2bi  ,  b3  »  bi  -  2  , 
ci  -  |  [l  -  (3.2)B  ]  , 

c2  -  3  -  2ci  ,  c3  -  Cl  -  2 
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where 


di  ,,  PrB  ?£ 

B  sc  e*bi 


B «  .  J.  ±  £1 

8c  8m  M 

It  should  be  noted  that  the  functions  (4.3*0  through  (4.36)  repre¬ 
sent  physically  Meaningful  profiles  only  when  the  coefficients  are  all 
real.  Upon  exsMl nation  of  the  shore  expressions!  it  Is  evident  that 
when  B  >  2/3  no  real  coefficients  are  obtainable.  More  specifically! 
the  Method  used  here  is  strictly  applicable  when  and  only  when  B  £  2/3. 
Fortunately!  in  frosting  problsn  (vx^erx? )  is  usually  very  •mil,  of 
the  order  of  1CT  ,  and  the  condition  B  £  2/3  is  fulfilled. 

With  the  aid  of  Eqs.  (4.34).  (4.35),  (4.36)  and  (4.43),  the 
integral  equations  (4.29),  (4.30)  and  (4.31)  can  now  be  rewritten  as 


(1-B")ai  -  —  —  (u*8**)  +  8*U  —  (4.44) 

vu  dx  dx 

(1-B*  )bx  -  A  (ubj##)  (4.45) 

(i-b)cx  -  5iia  A  (uBm**)  (4.46) 

V  dx 

where,  after  carrying  out  the  indicated  integrations  in  Eqs.  (4.27), 
(4.28),  (4.32)  and  (4.33), 

T  ■£  "  i  *  5 

¥ *  T  £  *  ft)  +  *«*>(&  ■  ft) 

+  V(*x-2)(^  -  ±)]  (4.49) 
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and 


^■t/01(1-h)w© 

♦  A.a(i-2)(^  - 


-  S) + -  is) 


_c 
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)] 


(U.50) 


Thelaat  two  expressions  are  strictly  ralid  for  8g  8  and  8*  £  8. 

Solution  of  the  simultaneous  ordinary  differential  equation* 

(4.44),  (4.45)  and  (4.146)  for  the  general  boundary  layer  flow  with 
preaaure  gradient  preaenta  cooalderable  difficulty  and  will  not  be 
attempted.  Tor  the  particular  caae  of  stagnation  flow  under  considera¬ 
tion,  advantage  nay  be  taken  of  the  information  available  In  the  litera¬ 
ture.  m  the  absence  of  mass  transfer,  i.e.  when  Vf  -  0,  the  boundary 
layer  thickness  &  In  a  stagnation  flow  is  constant,  a  well-established 
fact  which  can  be  readily  seen  from  Eq.  (4.13).  An  inspection  of  the 
existing  solutions  for  mass  transfer  in  lM« flows  over  flat  plate 
and  in  free  convection  from  a  vertical  plate  indicates  that  the  presence 
of  vf  does  not  effect  the  functional  relation  between  5  and  x.  On  the 
basis  of  these  above  observations,  it  can  be  considered  that  6  Timmins 
constant  in  the  stagnation  flow  when  Vf  ■  0.  Since  dg  -  8g/8  and 
^  ■  &■/&  are  independent  of  x,  it  follows  that  6*,  8**,  &£**  and  Q^** 
are  all  constant.  Hance  Eqs.  (4.4)  through  (4.46)  reduce  to 


a*.*. 


(**.») 


(l-B*  )bi  -  Mnst  3- 
o 


(l-B)ci  -  >BcAn  — 
5 


(4.52) 

(4.53) 


therein  8*/®  and  8**/6  are  known  in  terms  of  %i  from  Eqs.  (4.47)  and 
(4.48),  while  8j /8  and  8j*/8  cone  from  Eqs.  (*>  >49)  and  (4.50). 

A  relation  between  B"  and  ai  may  be  found  from  (4.51)  by  elimina¬ 
tion  of  A  with  the  aid  of  the  first  relation  of  Eq.  (4.43),  namely. 


B"  ■ 


1  -  ai(S. 


aiUL 

O  "  " 


8 *\ 
8  ' 


(4.54) 
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It  has  been  demonstrated  by  Spalding  that  the  approximation  ^Mli 
be  substituted  in  the  brackets  in  the  shore  equation  without  causing 
appreciable  error.  Thus 

£  .  a*  Hi  fi  ,(,1'c>)  (kb  *Sao)l  (i».. 


Dividing  Fq.  (4.51)  by  (4.53)  using  Eq.  (4.56),  there  results  an 
expression  for  uarsely, 


Le 


(^f  .  d-B'  ftl  Sc  .  (1-B')ti 
VAi/  (l-B)ci  Pr  (l-B)oi 

(Le  is  the  Lewis  number) 

and  then  substituting  this  in  the  equation  defining  B' ,  that  is 
B'  ■  ~  The  resulting  expression  is 


(1-BQb!  (l-B)Cl 
Le^B'bj*  .  (Bex)4 


(4.59) 


4.4  BElERMniAXI<H  OF  THE  HEAT  AND  MASS 
TRANSFER  COEFFICIENTS 


The  determination  of  the  heat  and  mass  transfer  coefficients  in 
the  presence  of  mass  transfer  is  rendered  more  difficult  than  in  the 
case  where  mass  transfer  is  absent.  As  will  be  seen  below,  the  heat 
and  mass  transfer  rates  per  unit  surface  area  depend  on  the  first 
coefficients  of  the  temperature  and  material  profiles  as  well  as  the 
thicknesses  of  the  temperature  and  material  boundary  layers.  The  latter 
in  turn  are  related  to  the  Telocity  boundary  layer  thickness.  Since 
it  is  the  quantity  B  which  is  normally  specified,  it  is  desirable  to 
express  the  heat  and  mass  transfer  coefficients  in  terms  of  B.  However, 
this  can  only  be  accomplished  when  it  is  possible  to  express  explicitly 
the  coefficients  of  the  assumed  profiles  (4.34)  through  (4.36)  in  terms 
of  B.  It  is  evident  from  the  preceding  section  that  such  explicit 
relations  are  difficult  to  find.  Consequently,  the  heat  and  mass  trans¬ 
fer  coefficients  can  only  be  obtained  either  numerically  or  graphically 
for  given  values  of  B. 

The  rate  of  heat  transfer  (sensible  heat  transfer)  per  unit  surface 
area  is,  as  before,  given  by 


q  * 

Substituting  the  expression  for  the  temperature  gradient  from  Eq.  (4.35 ) > 
this  leads  to 

q  -  -kbiCr.-T^/br 

Using  the  equations  defining  X,  B',  and  B",  the  above  equation  may  be 
put  in  the  following  form 

q  "  ’  §7  ai(T--Tf )  -  -H*  § 7  *i(T«r Xf )  (4.60) 
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In  terns  of  the  heat  transfer  coefficient  defined  previously. 


h  ■  kPr  J ^  ai  (4.6l) 

1  \v  B' 

The  procedure  for  calculating  heat  transfer  rates  or  coefficients 
at  the  two  -dimens  tonal  stagnation  point  is  as  follows: 

1.  The  quantity  B  is  ascertained  from  the  data  of  the 
problem. 

2.  The  coefficient  c,  is  determined  by  the  seventh  rela¬ 
tion  of  Eq.  (4.43). 

3.  B'  is  obtained  from  Eq.  (4.59)  with  the  help  of  the 
fourth  relation  of  Eq.  (4.43). 

4.  B"  and  ai  are  found  by  solving  simultanelusly  Eqs.  (4.45) 
and  (4.58). 

5.  Hie  dimensionless  quantity  A  follows  from  the  first 
relation  of  Eq.  (4.43). 

6.  Finally,  the  heat  transfer  rates  of  coefficients  are 
calculated  from  Eq.  (4.6o)  or  (4.6l). 

For  given  values  of  Pr  and  Sc,  approximate  values  of  the  heat  and 
mass  transfer  coefficients  or  rates  may  be  found  rapidly  by  the  use  of 
a  graphical  method.  The  expression  on  the  right-hand  side  of  Eq.  (4.59) 
is  plotted  in  Fig.  7*  Because  of  the  relation  (4.59),  Interpolating  in 
Fig.  7  enables  bi  and  hence  B*  to  be  determined.  To  determine  Bn  and 
ai,  it  is  convenient  to  plot  both  sides  of  Eq.  (4.58).  The  left-hand 
side  is  plotted,  with  the  aid  of  Eq.  (4.54),  in  Fig.  8  for  various 
values  of  ax;  while  the  right-hand  side,  excluding  the  quantity  Sc2,  is 
plotted  in  Fig.  9  for  various  values  of  ex  with  ax  as  parameter.  Since 
I  is  dependent  on  both  ax  and  ci,  Eq.  (4.58)  has  to  be  solved  by  a 
method  of  successive  approximation.  Let  the  quantity  on  the  left-hand 
side  of  Eq.  (4.58)  be  designated  by  M.  For  any  given  value  of  B  and 
hence  Cx,  a  value  of  M  may  be  obtained  by  means  of  Fig.  9  for  an  arbi¬ 
trarily  selected  value  of  ax.  Corresponding  to  this  value  of  M,  an  ax 
is  read  from  Fig.  8.  If  this  value  of  ax  differs  from  the  selected 
one,  go  back  to  Fig.  9  and  obtain  another  value  of  M,  corresponding  to 
the  new  value  of  ax.  Repeat  the  procedure  until  identical  value  of  ax 
is  obtained  for  successive  iterations,  nils  value  of  ax  is  then  the 
solution  of  Eq.  (4.58)  for  the  given  value  of  B.  With  ax  thus  found, 

B"  follows  from  Fig.  10,  in  which  B"  is  plotted  versus  ax  using  the 
relation  (4.54).  Inserting  the  values  of  ax,  B*  and  B"  in  Eq.  (4.6l), 
the  heat  transfer  coefficient  is  readily  calculated. 
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Fig.  8  Diagram  for  uaa  with  aquation  (4.58) 
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Similarly,  the  mass  transfer  rate  is  given  by 

m  *  “  PD(-^-)yo  “  “  PDci(vl00-vXf) 

Using  the  definition  of  B",  this  becomes 

m  =  -p  2|£  |_  aiCwi^-wif)  «  -  pDSc  |-  aiCwic-wif) 

(4.62) 

The  mass  transfer  coefficient  defined  in  Section  4.2  assumes  the  form 

*b-firr‘l -“'VIr*1  <4-63) 

The  technique  described  above  for  finding  the  heat  transfer  coefficient, 
omitting  the  step  of  determining  B',  applies  here. 
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CHAPTER  V 


HEAT  AND  MASS  TRANSFER  IN  A  TWO-DIMENSIONAL  WALL  JET 


5.1  INTRODUCTORY 

The  term  "well  jet"  was  first  Introduced  by  Qlauert37  to  describe 
the  flow  in  a  viscous  Jet  bounded  on  one  side  by  a  wall  and  on  the  other 
by  the  saae  fluid  at  rest.  Such  a  Jet  arises,  for  example,  when  a  free 
Jet  is  allowed  to  impinge  no  really  aganist  a  plate  and  spread  out  over 
it*  The  situation  envisaged  is  the  final  state  of  the  Jet  when  distances 
from  the  stagnation  point  (i.e.  the  point  of  intersection  between  the 
axis  of  the  1  wringing  Jet  and  the  plate)  are  sufficiently  large. 

The  lwpingaamnt  of  a  Jet  of  electric  wind  onto  a  plane  surface 
provides  another  example  of  a  wall  Jet.  The  measured  current  density 
distributions  over  the  plane  surface  in  a  positive  wire-plane  corona 
reveal  a  sharp  concentration  near  the  stagnation  point.  At  large 
distances  from  this  point,  the  current  density  is  essentially  zero  and 
the  pressure  gradient  caused  by  the  electro-static  force  is  negligible. 
Thus  the  conditions  of  a  wall  Jet  prevail  in  this  region. 

The  main  feature  of  a  wall  Jet  is  its  hybrid  structure  common  to 
both  the  free  Jet  and  the  ordinary  boundary  layer.  Thus  the  motion  of 
the  fluid  is  retarded  by  frictional  forces  near  the  wall  and  the  inner 
layer  of  the  flow  nay  be  expected  to  show  a  certain  structural  similar¬ 
ity  to  a  boundary  layer.  Boar  the  outer  edge  of  the  flow,  entrainment 
of  quiet  fluid  occurs  and  accordingly  the  outer  layer  of  the  flow  is 
likely  to  resenfcle  a  free  Jet  in  character. 

An  Important  contribution  to  the  theory  of  i— rfiw  vail  Jet  has 
been  made  by  Qlauert^  »*®.  He  found  a  similarity  solution  of  the  bound¬ 
ary  liyer  equations  governing  such  flow  in  an  Incompressible  fluid. 
Riley**  extended  the  theory  to  the  case  of  compressible  flow.  He  also 
obtained  similarity  solutions  for  the  energy  equation.  In  this  chapter, 
the  similarity  solution  of  the  momentum  equation  first  found  by  Qlauert 
is  briefly  discussed.  A  technique  used  by  Riley  in  obtaining  similarity 
solution  of  the  energy  equation  for  compressible  i— vail  Jet  is 
employed  to  seek  solutions  of  the  energy  and  diffusion  equations  for  two- 
dimensional  incompressible  laminar  wall  Jets  over  a  frosting  surface, 
under  the  simplified  assumption  that  the  normal  velocity  at  the  frost 
surface  is  negligibly  small.  The  effects  of  non-vanishing  normal  velo¬ 
city  at  the  frost  surface  on  the  flow  and  the  heat  and  mass  transfer  are 
investigated  using  an  approximate  method  based  on  the  integral  equations. 
The  coordinate  system  shown  in  Fig.  6  will  be  used  in  this  chapter. 


5.2  BASIC  EQUATIONS 

The  equations  governing  a  two-dimensional  laminar  jet  flowing  over 
a  plane  surface  will  now  be  considered.  In  view  of  the  fact  that  a  wall 
Jet  is  in  essence  of  combination  of  a  free  jet  and  a  boundary  layer  flow 
over  a  flat  plate ,  the  boundary  layer  approximations  may  be  used.  In 
the  absence  of  a  pressure  gradient,  the  momentum  equation  for  incompres¬ 
sible  flow  may  be  written 


u 


du 


v 

V2 


The  energy  equation  including  viscous  heating  is 


u 


& 


+  v  — 


d^T 
Pr  cy2 


(5.1) 


(5.2) 


Using  the  mass  fraction  wx  of  water  vapor  as  dependent  variable,  the 
diffusion  equation  is 


u 


_v_ 

Sc  by2 


(5.3) 


Here  the  thermal  diffusion  and  diffusion-thermo  effect  are  not  considered. 
They  are  important  only  when  the  temperature  and  mass  fraction  gradients 
are  extremely  large.  The  equation  of  continuity  for  two-dimensional 
incompressible  flow  is 


du  . 

¥  by 


(5.4) 


Turning  next  to  tbs  boundary  conditions,  it  should  be  noted  that 
in  connection  with  the  problem  of  frost  formation,  v  does  not  vanish 
at  the  surface  of  the  porous  frost  layer.  However,  the  magnitude  of 
Vf  is  in  practice  very  small  and  can  be  disregarded  in  forced  convection 
flow  without  causing  appreciable  error.  (The  case  where  Vf  differs 
from  zero  will  be  examined  later.)  With  this  in  view,  the  boundary 
conditions  appropriate  to  the  flow  of  a  wall  jet  over  a  frosting  surface 
are  exactly  the  same  as  those  for  the  same  flow  over  an  Impermeable  sur¬ 
face.  These  conditions  are 


aty-0  :  u  ■  v  ■  0 


T  ■  Tf 


wi  ■  Wlf 

(5.5) 


at  y  00  :  u  =  0 


T  =  T*  »  wl  ■  vi* 

(5.6) 
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In  order  to  seek  similarity  solution  of  the  momentum  equation, 
with  u  a  x®,  and  the  jet  width  6  a  xn,  two  relations  are  required  to 
determine  the  exponential  constants,  m  and  n.  One  of  the  relations  is 
obtained  from  the  momentum  equation  itself.  The  second  relation  is  to 
be  provided  by  the  general  nature  of  the  flow.  For  a  free  jet  con¬ 
stancy  of  the  total  momentum  in  the  x-dlrection  gives  the  second  rela¬ 
tion,  and  for  a  boundary  layer  flow,  one  of  the  constants  is  determined 
by  the  distribution  of  the  external  flow.  For  a  wall  jet,  however, 
there  is  no  such  obvious  quantity  to  consider.  In  the  connection  a 
valuable  relation  has  been  deduced  by  Glauert37.  Integrating  the  momen¬ 
tum  equation  with  respect  to  y  between  the  limits  y  and  ®,  multiplying 
the  resulting  equation  by  u,  integrating  again  with  respect  to  y  between 
limits  0  and  ®,  and  making  use  of  the  continuity  equation  (5.4)  and  the 
boundary  conditions  (5.5)  and  (.56),  the  following  relation  is  obtained. 


HU  u2dy  dy  «  constant J 


(5.7) 


5.3  VELOCITY  DISTRIBUTION 


Glauert  has  shown  that  a  similarity  solution  of  the  momentun  equa¬ 
tion  (5.l)  in  which  the  form  of  the  velocity  profile  does  not  vary  with 
x  is  available.  This  ary  be  written 


u  -  U(x)f '  (ij)  where  t|  -  -^y 


(5.8) 


Here  the  prime  indicates  differentiation  of  a  function  with  respect  to 
its  argument.  Equation  (5.4)  is  satisfied  by  the  introduction  of  the 
stream  function  f  such  that 


U-h.  ,  y  .  .5* 

is  a* 

With  the  aid  of  Eq.  (5*8),  this  may  be  integrated  to  give 

♦  -  U  6  f (tj) 


and  Eq.  (5.7)  shows  that 


U3  B2  ■  constant. 


Making  use  of  Eqs.  (5.10)  and  (5.H),  Eq.  (5.1)  then  reduces  to 

f'*»  -A  U'b2?"  -  U,ft2  f'2  a  0 

2v  v 

The  condition  for  the  existence  of  similarity  solution  then  is 

U*  82  *»  constant 


(5.9) 


(5.10) 


(5.1 1) 


(5.12) 


(5.13) 
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(5.14) 


If  the  constant  in  Eq.  (5.13)  is  chosen  to  have  the  value  of  -2v, 

Eq.  (5.12)  becomes 

ft*  •  +  ff"  +  2f’2  =  0 

The  solutions  to  the  simultaneous  equations  (5.11)  and  (5.13)  are 
U(x)  =  2A  [v(x+i)]  ,  8(x)  -(I)1*  [v(x+i)]  «/♦ 

(5.15) 

where  A  and  /  are  arbitrary  constants.  The  boundary  conditions  (5*5) 
and  (5*6)  require 

f(0)«f’(0)=0  ,  f'(«)  ■  0  (5.16) 

Multiplying  by  f  and  using  Eq.  (5.16),  Eq.  (5.14)  may  be  integrated 
to  yield 

f"  -  £  f ,2  +  f^ '  »  0  (5.17) 

Multiply  by  f  and  integrate  again.  Thus 

f  f  *  f '  +£  fa/2  ,  constant  (5.18) 

3 

Without  loss  of  generality.  It  may  be  required  that  f(°°)  ■  1  and  the 
constant  in  Eq.  (5 .18)  may  be  taken  to  have  the  value  of  2/3. 

It  is  convenient  to  make  the  transformation  f  ■  g2.  Then  f'  ■  2gg' 
and  Eq.  (5*18)  becomes 

g'  -i  (1-g3) 

Upon  integration,  it  is  found 

q  +  V3  tan'1  £-6  (5.19) 

1-g  2-tg 

The  arbitrary  constant  A  is  to  be  determined  by  Eq.  (5*7).  It 
can  be  verified  that 

R  "  A*f0  f ’[  /,  f ,2<iTl]  dTl  "  |  A2  (5.20) 
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Hence  we  nay  write 


A  -  ±(5R/2)* 


t  ■  £40Rv(x+l)J1/4  f(i)) 


(5.21) 


n 


i 


5R 


32va(x+J)a 


y 


The  quantity  R  and  the  constant  i  nay  he  estimated  roughly  as 
follows.  If  the  wall  jet  is  the  result  of  a  free  jet  Impinging  on  a 
flat  plate,  the  magnitude  of  R  nay  he  obtained  from  the  conditions  in 
the  free  jet  as 

R  *  £( typical  velocity )(volune  flow/unit  length)2  (5.22) 


The  constant  length  i  can  he  determined  by  equating  the  maaimum  velocity 

at  x  ■  0  as  given  hy  Eq.  (5.21  to  the  marl  mm  velocity  in  the  imping¬ 
ing  jet. 


5.4  SIMILARITY  SOLUTION  OF  THE  ENERGY  EQUATION 

Having  solved  the  nomentvm  equation,  the  energy  equation  nay  now 
he  considered.  Making  use  of  the  results  of  the  previous  section,  the 
following  relations  are  obtained: 


u  ■ 


■J* 

2v(x+i) 

v  * 


]*  '■<”  •  &  ■ 

-  4<40Rv)l/‘  (*+!)-»/.  f(t|)  +  1(1|0H»)V4  (x+i)-5/* 


f’W  » 


dr 

dx 


^  +  (x+t)~xn  ^ 

abc  dq  &  &  4  1  ^ 


&  -  r  5R  1 1/4  dr 


32va(x+l) 


and 


<fT  ,  r_5R_t  &T 

^2  |_32v3(x+i)3  J  V 

Thus,  with  (x,tj)  as  independent  variables,  the  energy  equation  (5«2) 
may  be  written  as 


^  -  4Pr(x+/)f '  —  +  Prf—  » 
drj2  dx  dr) 


Pr 

CP 


[ 


-2—1 

2v(x+l) J 


^if  2 


(5.23) 


with  the  boundary  conditions 

T  -  Tf  at  y  -  0  ,  T  =  Ta  at  y  »  « 


(5-24) 


The  general  solution  of  Eq.  (5*23)  consists  of  a  particular  solution 
and  suitable  complementary  solutions.  The  particular  solution  may  be 
found  in  the  form 

T  -  T0  +  Ko(x+-er1e0(4)  (5-25) 

where  To  and  Kb  are  constants.  If  Ko  is  chosen  as  Ko  e  ■  5R/2cp,  then 
00(11)  satisfies 

ej  +  Pr(f0o+4f  0O)  =  Prf”2  (5.26) 

Eq.  (5.23)  also  has  an  infinity  of  complementary  functions  of  the  form 

T  =  Ki  +  K2(x+i)“a  0n(Ti)  (5.27) 

with  0n(il)  satisfying 

0J1  +  Pr(fO^ +4af ’0n)  =  0  (5.28) 

Here  a  is  a  constant  to  be  determined  by  appropriate  boundary  conditions. 
A  complete  picture  of  the  temperature  distribution  in  a  particular  prob¬ 
lem  can  thus  be  obtained  by  adding  to  the  particular  solution  (5.25) 
appropriate  supplementary  solutions  of  the  type  (5 .27). 

To  examine  the  effects  on  temperature  distribution  of  viscous 
heating  in  the  wall  jet,  it  is  required  to  seek  a  solution  of  Eq.  (5*23) 
in  the  form  of  Eq.  (5*25)  where  00(4)  satisfies  Eq.  (5*23)  with 


0O(«)  -  0 

and  0o(O)  “  0  for  constant  wall  temperature  Tw 


(5.29) 
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since  T  ■  T„  at  t)  ■  To  ■  T*.  When  the  Prandtl  nisaber  Pr  is  unity, 
Eq.  (5*26)  subject  to  (5.29)  has  the  solution 


8o(<i)  -  i  t’*w 


(5.30) 


For  arbitrary  Pr,  no  solution  in  the  closed  fora  has  been  found.  The 
equation  has  been  integrated  numerically  by  Riley38  for  Pr  -  0.72,  the 
value  appropriate  for  air.  He  gave  0o(O)  m  0.0036.  Thus,  the  rate  of 
heat  transfer  across  the  vail,  per  unit  area,  is 


-  Kok(x+irx^  e'o(0) 


-  0.0036  (5/4 )  (k/cpV  )  [5R8/2v3(x+J)7  y/* 


(5.31) 


To  study  the  effect  on  the  temperature  distribution  of  maintaining 
the  vail  at  a  constant  teaqwrature  Tf,  a  complementary  solution  (5.27) 
vith  a  -  0  is  required.  To  satisfy  the  boundary  condition  (5.24), 

Eq.  (5.27)  nay  be  put 


where  0i(t))  satisfies 


T  -  Tf  +  (T^-Tf  )0i(tj) 


♦  Prf  0x  ■  0 


(5.32) 


(5.33) 


0l(o)  -  0 

fll(-)  - 1 


(5.34) 


The  solution  of  Eq.  (5*33)  subject  to  (5-34)  may  be  Immediately  obtained. 
Direct  Integration  of  Eq.  (5.33)  gives 

®x(n)  ■  0i(O)axp  [“PrJ^  fdqj 
Integrating  once  more,  it  is  found 

0i(tj)  -  0i(«)  -  ©1(0 )f  exp^-Pr^  fdtjjdi] 

Using  Eq.  (5.34)  to  obtain 

ei(o)  -  vf-M H  dn  (5.35) 


(5.35) 
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the  final  solution  of  Eq.  (5.33)  becomes 


0i in)  =  1  ~  j  exp[-pr Jq  fdT,]dTJ//0  jo  fdT)]dT) 

In  Section  5*3  we  have  f  =  g2  whore  g'  ■  (l-g3)/3,  hence 

fdq  «  -f  &  =  ln(l-g3) 

J0  J0  (l-g3) 

Thus  in  terms  of  g,  Eq.  (5*36)  may  be  written 

Si  -  1  -  j£  (l-g3)1*-1^/ /  ( l-g 3)Pr-1dg  (5.37) 

When  Pr  =  1,  Eq.  (5*37)  may  be  readily  integrated  to  give 

0i  -  B  (5-38) 

When  r-  4  1,  writing  s  =  (l-g3),  Eq.  (5*37)  then  becomes 

Ox  ~  -J  s^>r"1(l-s)“2/3d s//  sPr_1(l-s)_  ^  da 

where  the  denominator  is  the  beta  function  B(Pr,l/3)  and  the  numerator 
is  the  incomplete  beta  function  Bs(Pr,l/3).  Therefore, 

0i  -  1  -  Bs(Pr,l/3)/B(Pr,l/3)  (5.39) 

The  rate  of  heat  transfer  per  unit  area  across  the  wall  is  given 
by 

q  =  -  Kfjpy^Q  ’  ■  k(T» ' rf )  [5R/32v3(x+/)3]  A/4  ei(0) 

(5.4o) 


From  Eq.  (5.35), 

ei(C)  =  l/B(Pr,l/3)  -  r(Pr+l/3)/r(Pr)r(l/3) 

where  the  relation  between  the  beta  and  gazana  functions  has  been  used. 
For  Pr  ■  0.71,  9i(0)  =  0.273*  The  local  heat  transfer  coefficient 
defined  as 

h(x)  =  -  q(x)/(TC0-Tf) 
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(5.41) 


Is  readily  found  from  Eq.  (5.40)  to  be 

h(x)  -  k6x(0)(5R/32va)1/4(x+l)"3/4 

Defining  the  local  Nusselt  number  as 

Nu(x)  -  hx/k 

it  follows  from  Eq.  (5.4l)  that 

Nu(x)  -  01(O)(5R/32va)V4x(x+l)"*/4  (5.42) 


5.5  SOLUTION  OF  THE  DIFFUSION  EQUATION 

Upon  comparison  of  Eq.  (5*3)  and  (5.2),  it  is  evident  that  complete 
analogy  exists  between  the  two  equations  when  the  viscous  dissipation 
term  in  the  energy  equation  is  absent.  By  '’complete  analogy"  it  is 
meant  that  the  two  equations  are  similar  in  form.  Perfect  analogy 
between  the  transport  processes  through  the  boundary  layers  requires 
also  that  the  boundary  conditions  be  identical.  Fortunately,  this  is 
the  case  for  the  problem  under  consideration,  as  can  be  seen  from  the 
boundary  conditions  (5.5)  and  (5.6).  Therefore  the  solution  of  Eq.  (5-3) 
can  be  obtained  directly  from  the  results  of  Section  5*4  by  replacing 
the  Prandtl  number  Pr  by  the  Schmidt  number  Sc. 

In  terms  of  the  independent  variables  x  and  q,  the  diffusion 
equation  (5*3)  assumes  the  fora 

-  4Sc(x+A)f '  ^1  +  Scf  -  0  (5.43) 

dq2  dx  dtj 

and  the  boundary  conditions  become 

wi  -  wif  at  q  ■  0  ,  wi  ■  wx  at  q  ■  *  (5.44) 

Introduction  of  the  dimensionless  variable 


dq)  - 

"l.-Wlf 

in  Eqs.  (5*43)  and  (5.44)  leads  to 


with 


9"  <f  Scccp*  =  0 

Cp(O)  ■  0 
q>(»)  ■  1 


(5.45) 
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n 


•-  «.  i.  r  ■■■  ■  >  ..  .  1 

r 

n 

n 

The  solution  is,  from  Eq.  (5.39) » 

1 

9*1“  Bs(Sc,l/3)/B(Sc,l/3) 

(5.46) 

0 

The  rate  of  mass  transfer  per  unit  area  across  the  wall  is  given 
by 

n 

m  =  -  -  -  pD(wi08-wxf, )  [  5H/32v3(x+i)3 

7y=o  L 

(5.*tT) 

0 

j  where 

n 

fl  9'(0)  =  l/B(Sc,l/3) 

1 

For  air  in  the  temperature  range  between  0°F  and  100 °F,  Sc  *  0.53  and 
qp'(O)  *  0.258.  In  analogous  to  the  case  of  heat  transfer,  the  local 
mass  transfer  coefficient  is  defined  as 

0 

i*a(x)  =  -  m(x)/p(\ixar^Lf) 

and  the  local  Sherwood  number  is  defined  by 

Sh(x)  =  hox/D 

j 

It  follows  from  Eq.  (5.47)  that 

hjx)  -  Dq)'(0)(5R/32v3)1A(x+i)-3/4 

(5.48) 

-1 

and 

Sh(x)  =  qp'(0)(5R/32v3)1/4x(x+/)”3/4 

(5.49) 

— 1 

j 

5.6  TWO-DIMENSIONAL  WALL  JET  WITH  Vf  ^  0 

] 

In  the  foregoing  sections,  it  has  been  demonstrated  that  closed 
fora  solutions  are  available  for  the  momentum,  energy  and  diffusion 
equations  when  the  transverse  velocity  v  vanishes  at  the  frost  surface, 
i.e.  when  Vf  *  0.  Ifce  situation  is  quite  different  when  the  condition 

Vf  ^  0  is  taken  into  consideration.  In  order  to  clarify  the  difficulties 
that  arise,  the  governing  equations  (5*1)  through  (5.4)  will  be  re¬ 
examined  subject  to  the  new  boundary  conditions  that 

1 

n 

J 

at  y  »  0  :  u*v,  v=vf,  T=Tf, 

Wi  -  wif 

□ 

(5.50) 

| 

at  y  =  *  :  u  =  0  ,  T  =  T  ,  wi 

■  Wloo 

(5.5D 

LJ 

n 

10 

] 

D 

The  relation  (5. 8)  will  he  retained,  but  the  stream  function  ♦ 
defined  by  Eq.  (5.9)  will  now  be  assumed  to  take  a  different  form, 
namely, 


t  -  U(x)6(x)f(n)  +  to(x)  (5.52) 

The  last  term  of  Eq.  (5.52)  has  been  introduced  to  take  care  of  the 
condition  v  ■  Vf  at  y  ■  0.  It  follows  from  Eq.  (5*52)  that 

t  -  -  ffc  -  -  (U6)'f  +  US’  f  -  *0  (5.53) 


and 


vf  -  to  00  (5»5*0 

With  the  help  of  Eqs.  (5.8),  (5-53)  and  (5-54),  Eq.  (5.1)  leads  to 

f***  4-  Ofcff"  -  Qfef*2  +  dbf"  -  0  (5.55) 

where 

,  Cfe  .22*  and  (5.56) 

V  V  V 

are  constants.  Because  of  the  relation  between  and  the  material 
profile  as  indicated  by  Eq.  (4.1l)  the  momentum  and  diffusion  equations 
axe  now  coupled.  Undoubtedly  this  makes  the  exact  solution  much  more 
difficult  to  be  obtained . 

What  makes  the  situation  worse  is  the  fact  that  the  relation  (5.7) 
is  not  valid  any  more.  Instead  the  same  procedure  yields 

Scjf  U[j£  u2dyjdy  =  vf  ■  vfj^  u2dy  (5*57) 

which  does  not  provide  the  additional  relation  required  for  the  deter* 
mination  of  the  dependence  on  x  and  U  and  8.  That  this  is  so  can 
readily  be  seen  by  making  the  substitutions  U  a  x®  and  6  a  x®  in  Eqs. 
(5.56)  and  (5.57).  It  appears  that  further  investigation  is  necessary 
in  order  to  determine  whether  similar  solution  still  exists  in  this 
case.  NO  efforts  will  be  made  along  this  direction  in  the  present 
study. 

An  approximate  estimate  of  the  effects  of  vf  on  the  heat  and  mass 
transfer  processes  can  be  made  using  the  integral  method.  Two  cases, 

Vf  differing  considerably  from  zero,  will  be  discussed  separately  in 
Sections  5.7  and  5*8. 
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5.7  TWO-DIMENSIONAL  WALL  JET  WITH  Vf 
SMALL  BUT  NOT  ZERO 


To  examine  the  influence  of  small  vf  on  the  heat  and  mass  transfer 
processes  in  a  two-dimensional  laminar  wall  jet,  a  method  used  by 
Whitehurst40  and  by  Barron24  in  consideration  of  frost  formation  under 
free  convection  conditions  will  be  employed.  The  integral  equations 
may  be  obtained  in  the  same  manner  as  in  Chapter  IV  from  Eqs.  (5.1) 
through  (5.4).  They  are: 


d 

5 


(5.58) 


•0)udy  =  — 
Pr 


if«- 


y*C 


Vf 


(5.59) 

(5.60) 


The  boundary  conditions  which  are  to  be  satisfied  by  the  velocity, 
temperature  and  material  profiles  are 


at  y  =  0  , 

at  y  =  8  , 


u  »  9  =  <p 

8u 


u  ■  =  0 


(5.61a) 

(5.6lb) 


aty-Bj  ,  0=1  ,  —  =  0 

v 


at  y  »  8a  ,  <p  =  1 


£*E  =  0 


(5.61c) 

(5.62d) 


It  should  be  noted  that  the  method  ignores  the  boundary  conditions  that 

8©  v  d2©  ,,  dcp  _  v  d^q 


at  y 


du  d^ 


v  dg0 
Vf  "  Pr  dp  * 


v  era 
8c  d^2 

(5.62) 


The  velocity,  temperature  and  material  profiles  which  satisfy  the 
conditions  (5.61)  axe 


u 

u(x) 


(5-63) 
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fl-  2-L-(.L) 

(5.64) 

*  =  2  f  -  (*-) 

(5.65) 

Substituting  Eqs.  (5-63)  through  (5*65)  in  Eqs.  (5.58)  through  (5.60), 
carrying  out  the  indicated  integrations  and  making  use  of  Eq.  (4.11), 
there  results 

-i-  -i  (u=8)  .  .  JS 

105  dx  '  '  8 

(5 .66) 

d  _  v  2  v  2B 

(U&T  '  PF 

(5.67) 

S  (U8"">  -  B  k  - B  g  *  +  (1-B) 

(5.68) 

where 

^  ■  g  C {i-*'  s .  ■  -  &  * 

(5.69) 

and 

sf- iC 

♦  i  A,2) 

60  411  f 

(5.70) 

As  noted  earlier  the  relation  (5*7)  no  longer  holds  true  if 
differs  from  zero.  However  when  vf  is  very  small,  it  may  be  assumed 
that  Eq.  (5.7)  still  applies.  Using  the  velocity  profile  (5*63)  in 

Eq.  (5*7)  there  results 

U*  62  =  2500R 

(5.71) 

Equations  (5 .66)  and  (5«7l)  may  be  solved  by  means  of  the  substitutions 


U(x)  ■  DiCx+Z)01  and  8(x)  =  Dz(xH)n 
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Eqs.  (5.66)  and  (5.71)  then  become 


jir  (iWnjD^Cx+i)®*-1  =  -  j£  (x^f- 


Dx^Cx+Z)3114211  =  2500R 


Since  these  two  equations  must  be  valid  for  any  value  of  x,  each  term 
must  be  in  the  same  power  of  x.  Therefore 


2m  +  n-  l=  m-  n 


3m  +  2n  =  0 


This  gives 


m 


—  and  n  ^ 
2  4 


Thus  there  results  from  Eqs.  (5/$)  and  (5*7l) 


D1D22  «  420 


eu, 


D;.^2  =  2500R 

These  two  equations  may  be  solved  for  Dx  and  D2 


D2 


r  5R  1-1/4 

Li6(21v)3J 

.[  125R  ]* 

(5.72) 

Ulv(x+/)J 

=  r  5r  ■  -1-1/4 

(5.73) 

Ll£(21v)3(x+Z)3  J 
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And  hence 


The  linked  equations  (5.67)  and  (5 .68)  may  now  be  solved  for  Ap 
and  Aaj  which  will  be  assumed  to  be  independent  of  x.  With  the  aid  of 
Eqs.  (5*72)  and  (5*73),  Eqs.  (5.67)  and  (5.68)  may  be  rewritten 

Md  a**  (i _  i  4*  +  JLa  (5.74) 

4  1  Vl2  15  ^  60  r  /  Pr  Ar  Sc  Aa 


M2?  a1 2  (  L  .  i  a  +  _L  a  2)  =  -§_  -L  (l-B) 

i.  ^  V12  15  411  '  Sp  a_  '  ' 


Sc  ^ 


(5-75) 


When  Aa  «*  1,  a  very  accurate  approximate  solution  for  An  nay  be  obtained 
from  Eq.  (5*75)  by  letting  4  =  1  in  the  parenthesis,  namely, 

An3  =  24°  a  (l-B) 

ScDiD22 


a  aikB  r 

411  7  Sc 


Similarly,  Eq.  (5-74)  leads  to 


Af3  +  ^  -A-  Aj  -  -  ■  0 

1  7  Sc^a  1  7  Pr 


The  only  real  root  of  this  equation  is 

Ar=/_2.+f(JL_L.)3+(_2_)2l,Al1/3 

1  l7Pr  L'21  Sca/  '7Pr'  J  / 


2_  %  .  l'3 


(5.76) 


(5.77) 


The  heat  flux  at  the  frost  surface  is 

q  =  -kf§)  _  »  -2kia-^  -  -2k(Te#-Tf)/Ar5 


l6(21v)3(x+i) 


When  B  u  (that  is  vf  =  0),  Eq.  (5*78)  reduces  to 


q  »  -2k(T„-Tf )  £ 
For  Pr  =0.71 


5R 


i6(21v)3(x+i) 


1/3 


(5.79) 


,  =  -0.259k(T„-Tf)  [-^3] 


.32 v3(x+i)3. 

which  is  within  6 %  of  the  exact  solution  (5*39)  with  0i(O)  =  O.273. 
The  mass  flux  at  the  frost  surface  is 


m 


**  -<*>(&)  *  -2pD(wi00-wif  5/6^ 

/y=0 


*  -2pD(wi00-wif )  £ 


5*, 


I1/4  r  7Sc  I 

J 


l£(21v)3(x+i)3  J  L4(i-b) 


i/s 


(5.80) 


When  B  =  0, 


“  ’  -2CD(WI-W1f 5  L(21v)*(x.f)»  ]  ^  (T) 
For  Sc  =  0.53» 


1/3 


,.-o.»(V>f)p[J2iSfrt)!] 


‘/ 


1/4 


which  is  a  good  approximation  as  compared  with  the  exact  solution 
obtained  in  Section  5.5  with  <p'(0)  “  O.258. 

Substitution  of  Eqs.  (5.78)  and  (5*80)  in  the  definitions  of  the 
local  heat  and  mass  transfer  coefficients  and  the  local  Nusselt  and 
Sherwood  number  leads  to,  for  B  ^  0, 


hm(x)  *  2D  [ 


~lt  r  5r  ] 

1/4 

/at 

(5.81) 

Ll£(21v)3(x+/)3 J 

5R  I1/4 

f  7Sc  ]l/3 

(5.82) 

l6(21v)3(x+/)3-* 

Lu(i-b)  j 

!  tm3]  x(lt+/) 

(5.83) 
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and 


a(x)-2h<!i^]l/4  x(x+i)"s/4  [fel 


1/3 


(5.84) 


5.8  SOLUTIONS  FOR  IARGER  VALUES  OF  Vf 

In  Section  5.7,  it  has  been  assumed  that  Vf  is  so  small  that  the 
relation  (5*7)  applies  and  that  the  forms  of  the  velocity,  temperature 
and  material  profiles  are  not  affected  by  Vf .  These  two  assumptions  do 
not  apply  for  the  case  where  Vf  deviates  considerably  from  zero.  To 
account  for  the  effects  of  large  Vf ,  a  method  similar  to  that  used  in 
Section  4.4  will  be  employed. 

The  governing  equations  are  Eqs.  (5 .58)  through  (5.60)  and  the 
appropriate  boundary  conditions  are  Eqs.  (5.6l)  and  (5.62).  As  in 
Section  4.3,  the  velocity,  temperature  and  material  profiles  will  be 


assumed  to  be  of  cubic  form,  namely, 

* -“I 

(5.85) 

(5.86) 

and 

*  ■ ci  &  *  +  cs(£)3 

(5.87) 

By  consideration  of  the  boundary  conditions  (5.6l)  and  (5.62),  and 
inserting  the  value  of  Vf  from  Eq.  (4.11),  there  is  obtained  the  follow¬ 
ing  relations  determining  the  coefficients. 


1 

*  -2ai 

• 

» 

*3 

"  ax 

(5.88) 

bi 

■f  < 

(i-vr 

-  i  B' 
2 

) 

b2 

*  3  - 

2bi  , 

b3  ■ 

bi 

-  2 

(5.89) 

Cl 

-l(: 

L  -  >/lL~— 

Tb) 

B  V 

2  ' 

C2 

-  3  - 

2ci  , 

C3  = 

Cl 

-  2 

(5.90) 
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where 

B'  =  and  Bn  =  ■—£*-  (5.91) 

Sc8mbx  ScSjnai 

Equations  (5*88)  through  (5*90)  indicate  that  the  method  is  valid  when 
B  <  §-,  B*  <3  and.  B"  +  0. 

Substituting  Eqs.  (5-88)  through  (5-90)  in  Eqs.  (5*58),  (5*59)  and 
(5  •  90)  >  there  results  the  following  ordinary  differential  equations 


.m 


-2-  (Ua6) 

Uv  dx 


v  5  dx 


(5.92) 

(5.93) 


(1-B)cx  -  h£l  A  (us)  (5.94) 

v  6  dx 


where 


k*  Is)**  +(^'^)^2] 

(5.95) 


and 

0m*» 

a 


2_ 

15 


(5.96) 


Hie  relation  between  B'  and  B"  will  be  determined  first.  Dividing 
Eq.  (5.93)  by  Eq.  (5.94)  and  using  the  approximations  63  1  and  A«  1 
as  in  Section  4.3,  there  results 


,  (l-B')bi 
V^/  (l-B)cx 


Le 


Noting  Ar/An  »  ST/bm*  1*  follows  from  the  first  relation  of  Eq.  (5.91) 
that 


(l-B’ )bi  .  (l-B)ci 
(B’lli)3!*3  (Bci)3 

This  relation  is  the  same  as  Eq.  (4.59)  and  Fig.  7  can  be  used. 
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Turning  next  to  the  determination  of  the  relation  between  B"  and 
B,  rewrite  Eq.  (5.92)  as 

.  10^  .  8f  dU  .  8  A  (U8) 

ax  v  dx  v  dx 

and  note  that  Eq.  (5*96)  reduces,  when  ^  »  1,  to 

^  -  aAx(J=L  - 
5  ^  '210  105/ 


Dividing  Eq.  (5.92)  by  Eq.  (5*9*0  and  using  the  above  two  relations, 
there  yields 


(l-B)ci 


Scai(2TO  “  E&)  "(^  “  V 


(5.98) 


Substitution 
leads  to 


of  this  expression  in  the  second  relation  of  Eq.  (5. 91) 


(Bci)3 

(l-B)ci 


Sca(B"aj )' 


(— 

- 

/105  +8f  dU\ 

'210 

105' 

'ax  v  dx 

To  solve  Eq.  (5.99)  the  functions  U  and  6  must  be  known.  In 
principle  two  equations  are  required  to  determine  these  two  functions 
and  Eqs.  (5.57)  and  (5-92)  may  be  used  for  this  purpose.  As  noted  in 
Section  5.6,  substitution  of  the  relations  U  a  x?1  and  8  a  x11  in  these 
two  equations  fails  to  determine  uniquely  the  exponential  constants 
m  and  n.  Therefore,  a  different  approach  must  be  adopted.  In  view  of 
the  similarity  solution  obtained  in  Section  5o>  U  and  8  may  be  expressed 
in  the  following  form 


U(x)  -  diR(x)'A(x+i)"'A  (5.100) 

and 

6(x)  =  d2R(xr1/4(x+i)s/4  (5.101) 

where  R  is  now  assumed  to  be  a  function  of  x  rather  than  a  constant 
and  dx  and  d2  are  constants 


Defining,  as  in  Section  5.2, 


and  using  Eqs.  (5*85),  (5 -HO)  and  (5*101),  Eq.  (5.57)  becomes 

dR  m  _  vBai2cidi2  R 


105Sc4a  (x+i) 


Integration  of  this  equation  gives 


R  =  Ro(x+/rvB»i2cidi2A05ScAa 


(5.103) 


(5.104) 


Here  the  constant  Ro  can  be  interpreted  as  the  value  of  R  for  B  =  0. 
Because  of  its  physical  meaning,  Ro  may  be  roughly  estimated  by  means 
of  the  relation  (5.22).  From  Eqs.  (5.100),  (5.101)  and  (5.IO3),  the 
quantity  B2/v  dU/dx  is  found  to  be 


.  1  L  + 

v  dx  2  v  ' 


vBai2cidi2 

105ScAa 


(5.105) 


It  remains  to  determine  the  constant  di  and  ds.  Two  algebraic 
equations  relating  these  two  constants  may  be  obtained  by  inserting 
Eqs.  (5.100)  and  (5.101)  in  Eqs.  (5.92)  and  (5.102),  and  making  use  of 
Eq.  (5.103),  namely 

vBaiZci  diW  +  dife2  .  ^  (5.106 

3580!^  ai 


(5.106) 


Equations  (5.107)  and  (5.108)  can  be  solved  to  give 


83160 


V  [13020.x2  -  249480(g|ii; 


(5.109) 


_ (31ai3)2(83l60) _ 

v3  £1302 Oai2  -  24948o(^c^)]  3 


(5.110) 


^3^2  .  §2160 

31ai3 

Substituting  Eq.  (5.107)  in  Eq.  (5.106),  it  is  obtained 

dld22  =  (vBaj.2Cl/35Sc^) 

ai  31ai 


(5.107) 


(5.108) 
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The  heat  and  mass  transfer  rates  across  the  frost  surface  per  unit 
aree  are,  respectively, 


q  =  =  -kPrBciCTc-TfJ/B'Sc^  =  -kPrB"ai(T00-Tf  )/B’ 

'cfcr'y'O 


(5.111) 


and 


m 


-pd( 


aw 


^  )y=5  -Dci(wi-W1f)/Ai.6 


-pDScBnai(wioo-  wif)/B 


(5.112) 


The  local  heat  and  mass  transfer  coefficients  are,  respectively, 


and 


h  «  kPrBoi/B'Sc^  =  kPrBnai/B,6 


=  Dci/^B  *  DScB"ai/B6 


(5.H3) 

(5.114) 


The  solutions  obtained  above  are  undersirable  when  B  is  very  small, 
for  it  is  evident  from  Eqs.  (5*88)  and  (5.91)  that 


B"ai  =  4  and  An  =  Bci/4Sc 


When  B  =0,  An  =  0.  This  means  that  the  material  boundary  layer  does 
not  exist  when  B  =  0,  in  contrast  with  the  results  of  the  foregoing 
sections.  To  remedy  this  disadvantage  the  first  condition  of  Eq.  (5.62) 
may  be  released  and  the  second  relation  of  Eq.  (5.91)  defining  B"  must 
not  be  discarded.  Since  the  coefficient  ai  may  be  absorbed  in  the  con¬ 
stant  di,  there  is  no  loss  of  generality  by  putting  ax  =  1.  As  a 
result,  Eqs.  (5.98) ,  (5.IO9)  and  (5. 110)  become 


-Sc(-==- 

'•210 


(l-B)c 1 

H)  (105  +  T  |) 


105' 


83160 


v  £13020  -  24948o(^-i^  ) j 
(31)2(83360) 

v3  13020  -  24948o(fei  ci) 
'35ScAn' 


(5.115) 

(5.116) 


(5.117) 
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For  given  value  of  B,  B’  may  be  found  by  means  of  Eqs.  (5.97)  and 
(5.89).  while  Eq.  (5.H5),  with  the  help  of  Eqs.  (5.IO5),  (5.116)  and 
(5.117)  may  be  solved  for  4a*  With  4n  so  found,  it  follows  that  da., 
de,  R,  lt  q,  m,  h  as  well  as  ha  may  all  be  determined. 


82 


CHAPTER  VI 


EXPERIMENTAL  PROGRAM 


6.1  EXPERIMENTAL  APPARATUS 

The  initial  objective  of  the  experimental  investigation  was  two¬ 
fold.  In  the  first  place,  it  was  intended  to  determine  the  possible 
effects  of  electric  fields  upon  frosting.  Hie  second  objective  was 
to  provide  an  answer  to  the  question:  how  effective  the  electric 
field  effects,  if  there  were  any,  would  be  as  compared  to  other  mass 
transfer  processes  such  as  forced  or  free  convection?  Because  of  the 
exploratory  nature  of  the  current  study,  it  is  logical  to  use  as  simple 
a  test  setup  as  possible.  With  this  principle  in  mind,  it  was  decided 
to  use  free  convection  conditions  as  reference  state  to  which  the 
experimental  results  would  be  compared. 

A  schematic  of  the  test  setup  is  shown  in  Fig.  11  and  photographs 
of  the  apparatus  are  given  in  Fig.  12.  Hie  setup  consisted  of  three 
parts: 

a.  Test  plate 

b.  Refrigerating  system 

c .  Electrical  equipment 


6.1.1  Test  Plate 

The  test  plate  used  in  this  investigation  has  been  used,  prior  to 
the  frosting  work,  by  Velkoff  and  Miller2  to  study  the  condensation  of 
vapor  in  an  electrostatic  field.  The  plate  was  made  of  inch-thick 
sheet  copper  and  had  a  frosting  surface  measuring  6  by  9  inches.  To 
the  back  of  the  plate  was  soldered  £  inch  copper  tubing.  Refrigeration 
of  the  plate  was  accomplished  by  pumping  refrigerant  through  this  tubing. 

The  test  plate  was  installed  in  a  hole  cut  in  a  vertical  plexiglas 
wall  which  formed  one  side  of  a  wooden  box  enclosing  the  tubing  and  the 
back  of  the  plate.  The  wooden  box  was  filled  with  insulating  material 
to  insulate  the  plate  so  that  only  the  exterior  surface  of  the  plate, 
when  exposed  to  the  ambient  air,  was  available  for  frosting  and  heat 
transfer. 

Nine  iron-cons tantan  thermocouples  were  constructed  and  inserted 
through  the  back  of  the  plate  to  within  l/32  inch  of  its  front  to 
measure  the  temperature  of  the  frosting  surface.  To  eliminate  the 
effects  of  strong  electric  fields  during  the  tests,  a  small  amount  of 
insulation  was  provided  to  insulate  the  thermocouples  electrically 
from  the  plate.  These  thermocouples  were  located  along  the  sides  and 
at  the  center  of  4  3/4  by  5  inch  rectangle  centered  at  the  center  of 
the  plate  as  shown  in  Fig.  13. 
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6.1.2  Refrigerating  System 

The  refrigerating  system  used  to  refrigerate  the  test  plate  con¬ 
sisted  of  a  freezer,  a  refrigerant  reservoir,  a  pump  and  a  rotameter. 
Refrigeration  of  the  frosting  surface  was  accomplished  by  passing  the 
precooled  liquid  through  the  copper  tubing  on  the  back  of  the  test 
plate.  A  mixture  of  anit-freeze  (ethylene  glycol)  and  water  was  used 
as  refrigerant.  It  was  felt  that  such  a  mixture  was  adequate  to  fur¬ 
nish  the  desired  test  plate  temperatures .  In  addition,  it  is  not  haz¬ 
ardous  and  can  provide  lubrication  for  the  pump. 

A  £~ton  commercial  freezer  was  used  to  refrigerate  the  refrigerant. 
The  latter  was  contained  in  a  10  gallon  tank  installed  in  the  freezer 
to  serve  as  refrigerant  reservoir.  Liquid  was  drawn  from  the  tank, 
passed  through  the  copper  tubing  and  the  rotameter,  and  then  pumped 
back  to  the  tank.  Circulation  of  the  liquid  was  achieved  by  means  of 
a  rotary  pump.  Liquid  flow  rate  was  measured  with  the  rotameter. 

Tygon  plastic  hose  was  used  to  connect  the  copper  tubing  on  the 
back  of  the  test  plate  to  the  rest  of  the  flow  circuits,  thus  insulat¬ 
ing  the  plate  electrically  from  the  refrigerating  system.  Two  four- 
junction  copper -cons  tan  tan  thezmocouples  were  constructed  and  inserted 
in  the  plastic  hoses  at  the  inlet  and  outlet  to  the  copper  tubing  to 
measure  the  temperature  difference  between  the  inlet  and  the  outlet 
refrigerant.  Leads  of  the  thezmocouples  were  led  through  holes  drilled 
in  the  plastic  tubing  walls. 


6.1.3  Electrical  Equipment 

The  electrical  fields  used  during  the  tests  were  provided  by  a 
Sorensen  high  voltage  power  source  with  a  capacity  of  0-30  kv  and  up 
to  a  4  milliamperes .  The  output  from  the  high  voltage  supply  was  first 
connected  to  a  series  of  resistors  of  the  order  of  5  meg  ohms  total  and 
then  to  the  positive  electrode.  The  test  plate  which  served  as  cathode 
was  always  grounded  except  for  one  occasion  when  the  screen  electrode 
inserted  between  the  plate  and  anode  was  used  as  the  cathode.  In  that 
case,  the  test  plate  was  kept  neutral  electrically.  The  voltage  across 
the  field  was  measured  by  a  Singer  electrostatic  voltmeter.  Two  milli- 
ammeters  were  used  to  measure  the  current.  The  one  on  the  hot  side  was 
incorporated  as  part  of  the  power  supply.  The  other  was  connected  from 
the  test  plate  to  the  ground. 

The  electric  field  for  negative  corona  was  supplied  by  the  same 
power  source  by  reversing  the  polarity  of  the  output.  To  provide  an 
a.c.  field  for  the  a.c.  corona,  another  high  voltage  power  supply  was 
used.  A  0.006  inch  diameter  stainless  steel  wire  was  employed  as  the 
discharging  electrode.  The  wire  was  stretched  acioss  a  plexiglas  frame 
which  in  turn  was  mounted  on  a  movable  plexiglas  support.  Distances  of 
the  wire  from  the  test  plate  were  adjusted  by  moviig  the  support  or  the 


frame.  The  support  was  designed  so  that  the  vertical  position  of  the 
wire  could  be  adjusted  by  lifting  or  lowering  the  frame  relative  to  the 
support.  Several  parallel  wires  may  be  stretched  on  the  plexiglas  frame 
when  electrode  with  more  than  one  wire  is  desired. 


6.2  EXPERIMENTAL  PROCEDURE  AND  MEASUREMENT 
6.2.1  Mass  Transfer  to  the  Test  Plate 

Because  of  the  non-uniformity  of  the  frost  layer  thickness  and  of 
the  frost  density  along  the  entire  plate,  measurement  of  frost  thickness 
distribution  is  not  necessarily  meaningful  in  the  determination  of  mass 
transfer.  In  a  situation  like  this,  it  is  always  necessary  to  measure 
local  mass  transfer  rate.  Such  measurements,  however,  require  either 
a  measurement  of  the  local  concentration  profile  or  a  segmented  test 
plate  which  facilitates  the  determination  of  the  local  frost  deposition. 
The  former  requires  an  instrument  which  can  measure  accurately  the  local 
dry  and  wet  bulb  temperatures  in  the  thin  boundary  layer  in  a  strong 
electric  field.  The  latter  needs  a  special  design  of  the  test  plate  so 
that  each  segment  of  the  plate  can  be  maintained  at  the  same  temperature 
and  can  also  be  removed  easily.  Unfortunately,  technical  difficulties 
prevent  such  measurements  to  be  accomplished  in  the  limited  time  avail¬ 
able.  In  the  absence  of  a  convenient  method  to  deteraine  the  mass  trans¬ 
fer  rate  locally,  it  was  decided  to  measure  the  total  amount  of  frost 
deposited  on  the  test  plate  in  a  given  period  of  time. 

To  facilitate  the  comparison  of  the  total  mass  transfer  for  dif¬ 
ferent  electrical  conditions,  it  is  of  necessity  that  the  tests  be  con¬ 
ducted  under  identical  test  conditions,  i.e.  same  test  plate  tempera¬ 
ture  and  same  specific  humidity  and  temperature  of  the  surrounding  air. 
Furthermore,  to  determine  the  variation  in  time  of  the  mass  transfer 
under  the  influence  of  electric  fields,  it  is  also  necessary  to  maintain 
the  same  test  conditions  over  a  long  period  of  time.  Limitations  set 
by  the  present  test  facilities  render  this  impossible.  In  spite  of  the 
difficulty,  a  qualitative  trend  of  the  electric  field  effects  on  frost 
formation  can  be  obtained  by  shortening  the  test  period  so  that  it  may 
be  regarded  that  approximately  the  same  test  conditions  prevail  in  each 
test. 

The  test  period  selected  in  accord  with  the  capability  of  the  test 
setup  was  five  minutes.  Two  sets  of  tests  were  run  in  this  relation. 

In  the  first  set  of  tests,  the  distance  between  the  plate  and  the  wire 
was  kept  constant,  while  the  current  was  varied.  The  second  set  of 
tests  was  aimed  to  measure  the  variation  of  frost  deposition  with 
distances  of  the  wire  from  the  plate  for  a  given  current. 

The  test  plate  temperature  was  measured  with  the  nine  thermocouples 
inserted  to  within  l/22  inch  from  frosting  suriace  and  was  recorded  by 
a  self-blancing  potentiometer  on  a  continuous  recorder.  The  nine 
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thermocouples  were  connected  in  parallel  and  then  connected  to  the  re¬ 
corder.  A  reference  junction  in  the  open  air  was  used.  Temperature 
of  the  air  was  read  frcm  the  dry  bulb  thenncmeter  of  sling  type  psy- 
chrcmeter.  Specific  hunidity  of  the  ambient  air  was  determined  from 
the  dry  and  wet  bulb  temperatures  read  frcm  the  p&ychrcmeter.  Barom¬ 
eter  and  psychrameter  readings  were  taken  immediately  prior  to  the 
beginning  and  at  the  conclusion  of  the  test  run.  Average  value  was 
used  in  the  reduction  of  the  data.  The  variation  of  the  readings  be¬ 
fore  and  after  a  test  run  was  usually  small,  less  than  1°F. 

In  general,  it  took  about  eight  minutes  to  cool  the  test  plate 
down  to  a  steady  constant  temperature.  As  soon  as  the  test  plate  tem¬ 
perature  became  steady,  as  could  be  judged  frcm  the  temperature- time 
curve  on  the  recorder,  the  frost  already  formed  on  the  plate  during 
the  cool-down  period  was  scraped  off  and  the  electric  field  turned  on 
and  the  current  adjusted  to  a  predetermined  value.  After  five  minutes 
the  electric  field  was  turned  off  and  the  frost  scraped  off  for  measure¬ 
ment.  Another  run  of  a  given  teat  was  then  started  by  turning  on  the 
field  again  and  adjusting  the  current  to  a  new  value. 

The  frost  scraped  off  for  measurement  was  contained  in  a  light 
rectangular  tray  made  of  aluminum  foil.  The  tray,  being  weighed  prior 
to  containing  frost,  together  with  the  frost  was  then  weighed  on  a 
balance  and  measured  to  within  0.1  gram.  The  net  weight  of  the  frost 
was  obtained  from  the  gross  weight  by  subtracting  the  weight  of  the 
tray. 


6.2.2  Heat  Transfer  to  the  Test  Plate 

Measurement  of  heat  transfer  rate  under  various  electrical  con¬ 
ditions  presented  more  serious  difficulty  than  mass  transfer.  In  con¬ 
trast  with  mass  transfer  measurement  where  the  total  mass  transfer  can 
be  measured  for  a  given  period  of  time,  the  heat  transfer  rate  can  only 
be  measured  instantaneously  by  reason  of  the  unsteadiness  of  the  frost¬ 
ing  process.  For  the  comparison  of  heat  transfer  rate  for  a  variety  of 
electrical  conditions  to  make  sense,  it  is  desired  to  obtain  a  heat 
transfer  rate  versus  time  curve  for  each  electric  filed  while  other 
test  conditions  being  maintained  constant.  In  other  words,  a  longer 
period  of  test  is  desired.  That  this  requirement  can  not  be  met  with 
the  test  facilities  used  in  the  present  investigation  has  already  been 
noted.  Under  such  circumstances,  it  was  decided  to  sacrifice  the  con¬ 
stancy  of  the  test  conditions.  Each  test  lasted  30  minutes.  Since  it 
usually  took  4  to  5  hours  to  cool  down  the  refrigerant  in  the  freezer, 
the  test  conditions  might  have  changed  during  this  long  interval.  It 
was  hoped,  however,  that  in  this  manner  a  qualitative  observation  might 
be  obtained  as  to  the  effects  of  electric  fields  on  heat  transfer 
associated  with  frost  formation,  when  changes  in  test  conditions  were 
not  too  drastic. 


89 


To  start  a  test  run,  refrigerant  was  circulated  through  the  back 
of  the  test  plate  to  cool  down  the  test  plate.  When  the  test  plate 
temperature  reached  steady  state,  the  field  was  applied.  Readings  of 
refrigerant  flow  rate  and  the  temperature  difference  between  the  inlet 
and  outlet  refrigerant  to  the  test  plate  were  taken  at  two  minute  inter¬ 
vals.  The  rotameter  had  been  calibrated  prior  to  use.  The  calibration 
was  plotted  as  a  milliliter  versus  cm  of  rotameter  scale  so  that  the 
flow  rate  could  be  read  directly  from  the  curve  for  each  reading  on  the 
rotameter  scale.  The  two  four- junction  the  mo  couples  installed  in  the 
inlet  and  outlet  tygon  plastic  hoses  to  the  back  of  the  test  plate  were 
put  in  series  and  connected  to  a  portable  precision  millivolt  potentiom¬ 
eter.  The  accuracy  of  measurement  was  0.03°F. 

The  total  heat  transfer  rate  was  calculated  using  the  equation 

Q  =  WcAT 

where  W  and  c  are  the  mass  flow  rate  and  specific  heat,  respectively, 
of  the  refrigerant  and  AT  is  the  inlet  and  outlet  temperature  difference. 
For  the  purpose  of  comparing  heat  transfer  data  for  various  electrical 
conditions  provided  other  conditions  Eire  identical,  it  suffices  to  com¬ 
pare  only  the  product  of  volume  flow  rate  and  the  inlet  and  outlet 
temperature  difference  for  each  case.  On  account  of  this,  heat  transfer 
through  the  edges  of  the  test  plate  and  radiation  are  disregarded. 

In  the  course  of  a  test,  it  wels  discovered  that  the  current  some¬ 
times  fluctuated  considerably.  Therefore,  it  was  necessEiry  to  Euljust 
the  current  constantly  in  order  to  maintain  constEuit  current  throughout 
the  test.  In  fact,  the  current  tended  to  increase  as  the  frost  grew. 

It  appeared  interesting  to  observe  the  VEiriation  in  time  of  the  current. 
This  was  done  by  taking  readings  from  the  milliameter  at  two-minute 
intervals  while  the  frost  was  allowed  to  form  continuously.  Typical 
curves  thus  obtained  Eire  to  be  shown  in  Section  6.3*1** 


6.2.3  Velocity  Profile  Measurement 

Measurement  of  velocity  profile  is  of  primary  importance  in  the 
experimental  investigation  of  a  flow  field.  To  help  gain  an  insight 
into  the  actual  flow  field  under  the  impingement  of  an  electric  wind, 
efforts  were  made  to  mesisure  velocity  profiles  in  certain  regions  of 
the  flow  field.  Beceuise  of  the  limited  space  between  the  wire  and  the 
plate  and  because  of  the  strong  electric  field  present  in  the  *  iiate 
vicinity  of  the  stagnation  point,  measurement  of  velocity  in  the  elec¬ 
tric  wind  and  the  stagnation  flow  regime  would  have  been  difficult  and 
was  not  attempted.  Instead,  two  sets  of  measurements  were  made  in 
regions  away  from  the  stagnation  point.  The  first  region  under  investi¬ 
gation  was  the  wall  jet  flow.  Only  the  data  at  a  point  two  inches  from 
the  stEignation  point  Eilong  the  vertical  centerline  of  the  plate  were 
taken.  More  information  in  this  region  would  have  been  available,  were 


the  test  plate  sufficiently  large  to  permit  measurements  at  points 
further  downstream. 


The  second  set  of  measurements  was  directed  to  furnish  information 
on  the  velocity  distribution  of  an  electric  wind  formed  by  a  positive 
wire-. '.'lane  corona  in  the  absence  of  the  obstruction  by  the  plane  elec¬ 
trode.  Such  measurement  was  made  possible  through  the  use  of  a  screen 
as  the  cathode.  It  was  directly  behind  the  screen  that  the  actual 
measurement  was  carried  out. 

The  velocity  distribution  was  measured  by  means  of  a  0.06-inch  di¬ 
ameter  total  pressure  probe,  constructed  of  glass,  connected  to  a  micro- 
mancmeter  using  atmospheric  pressure  for  reference.  The  micrcmancmeter 
consisted  of  a  meniscus  in  an  inclined  tube.  The  displacement  of  the 
meniscus  along  the  tube  was  read  to  0.0001  inch  with  a  Gaertner  tran- 
velUng  microscope.  By  setting  the  inclined  tube  at  an  angle  of  approxi¬ 
mately  10°  relative  to  the  horizontal,  pressure  could  be  read  to  within 
0.00002  inch  butyl  alcohol. 

For  the  measurement  of  velocity  in  the  wall  jet,  a  probe  transverse 
was  made  perpendicular  to  the  plate.  Each  of  the  traverses  consisted, 
on  the  average,  of  19  readings  at  about  0.01  inch  intervals  near  the 
wall  and  intervals  of  about  0.05  inch  in  outer  regions  of  the  flow.  In 
the  measurement  of  electric  wind  velocity,  a  probe  traverse  was  made 
perpendicular  to  the  axis  of  the  wind.  An  average  of  15  readings  at 
intervals  of  0.1  inch  was  made  for  each  traverse.  The  accuracy  of 
setting  the  probe  was  0.001  inch. 

Flows  with  mixing  layers  are  inherently  unsteady.  In  this  investi¬ 
gation,  this  unsteadiness  was  manifested  as  a  timewise  fluctuation  of 
the  pressure  signal  from  the  total  pressure  probe,  particularly  in  the 
outer  mixing  zone  and  at  higher  currents.  In  view  of  this,  two  sets  of 
data  were  recorded  for  each  profile  by  traversing  the  total  pressure 
probe  twice.  Average  value  of  the  two  was  used  in  the  velocity  calcula¬ 
tions  . 

Apart  from  errors  in  pressure  measurement,  there  is  a  source  of 
error  at  the  outer  edge  of  the  flows  where  the  longitudinal  velocity 
tends  to  zero,  while  the  transverse  inflow  velocity  tends  to  finite 
value.  As  a  result,  there  is  a  yawed  flow  relative  to  the  total  pres¬ 
sure  probe.  Hie  points  within  a  few  total  pressure  probe  diameter 
from  the  wall  are  also  subject  to  error. 


6.2.4  Current  Density  Distribution  Along  the  Plate 

As  has  been  pointed  out  in  Chapter  III,  the  velocity  distribution 
of  the  electric  wind  depends  on  the  current  density  distribution. 

Being  unable  to  calculate  the  field  intensity  and  charge  density  dis¬ 
tributions  analytical] y  for  given  boundary  conditions,  it  is 
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desirable  to  seek,  experimentally ,  possible  correlation  between  the 
electric  wind  velocity  and  the  current  density.  It  was  for  this  pur¬ 
pose  that  the  current  density  measurement  was  carried  out. 

A  segmented  plate  with  each  segment  insulated  from  its  neighboring 
segments  is  required  in  order  to  measure  the  local  current  density 
along  the  plate.  To  achieve  this  with  the  existing  setup,  the  test 
plate  was  modified  by  covering  it  an  insulating  covering  painted  with 
horizontal  strips  of  conductive  silver  paint,  each  strip  being  separated 
from  the  adjacent  ones  by  a  thin  unpainted  area  of  l/8  inch  wide.  The 
painted  strips  were  so  distributed  that  the  entire  configuration  was 
symmetrical  with  respect  to  the  corona  wire. 

To  facilitate  the  measurement  of  the  current,  each  strip  was  con¬ 
nected  to  a  switch.  Each  switch  has  two  positions,  ground  and  meter. 

The  current  to  any  particular  strip  was  read  from  the  milliamneter  by 
putting  the  appropriate  switch  in  meter  position,  while  the  other 
switches  remained  in  ground  position.  The  switches  allowed  rapid  read¬ 
ing  of  the  current  to  each  strip  which  when  divided  by  the  area  of  that 
strip  giver  the  local  current  density  at  the  location  of  that  strip. 


6.3  EXPERIMENTAL  RESULTS 

Ike  preliminary  experimental  observations  which  formed  the  basis 
of  this  investigation  have  already  been  stated  in  Section  1.2  and  will 
not  be  repeated.  Only  the  data  related  to  the  measurements  described 
in  the  preceding  section  are  here  presented.  To  facilitate  interpreta¬ 
tion,  the  data  will  be  examined  in  the  reverse  order  relative  to  the 
list  of  measurements  in  Section  6.2. 


6.3.1  Current  Density  Distribution 

The  first  series  of  test  data  to  be  investigated  is  the  current 
density  distribution  at  the  plane  electrode  of  a  wire-plane  corona. 

Hie  results  are  plotted  in  Figs.  l4  to  l£.  In  each  case  the  abscissa 
are  the  distances  from  the  horizontal  centerline  along  the  plane  elec¬ 
trode.  Positive  and  negative  values  of  x  correspond  to  lower  and  upper 
half  of  the  plane,  respectively.  The  curves  in  each  figure  were  obtained 
by  varying  the  applied  voltage  while  keeping  constant  the  distance 
between  the  wire  and  the  plane,  d. 

All  of  the  current  density  distribution  curves  exhibit  the  same 
trend.  The  maximum  current  density  occurs  at  the  centerline  as  would 
be  expected  from  the  special  electrode  configuration  used.  As  x 
increases  in  both  directions  the  current  density  drops  rapidly,  forming 
a  bell-shaped  distribution.  Essentially  the  current  in  each  case  is 
confined  to  a  region  within  about  one  inch  from  the  centerline.  The 
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Fig.  14  Currant  dansity  distribution 
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Current  density  distribution  along  the  plate,  d  »  1  inch 


symmetry  of  the  current  density  distribution  with  respect  to  the  center- 
line  seems  to  suggest  that  the  corona  discharge  may  be  regarded  to  occur 
in  a  direction  normal  to  the  plane  electrode  without  causing  appreciable 
error. 


The  following  points  may  also  be  noted  in  connection  with  the 
current  density  distribution  curves: 

1.  The  area  over  which  the  current  is  received  (to  be  refer¬ 
red  to  as  the  current -area  hereafter)  seems  not  to  be 
affected  by  increasing  voltage  V  when  d  is  kept  constant. 

But  increase  in  d  tends  to  widen  this  area. 

2.  For  constant  V,  the  current-area  increases  and  the  maxi¬ 
mum  current  density  falls  as  d  increases. 

3.  The  values  of  x  for  which  the  current  density  is  half  the 
maximum  current  density  are  approximately  constant  for 
different  voltages  but  constant  d,  and  increases  with  d. 

Figure  17  shows  the  variation  of  the  maximum  current  density  with 
voltage  for  three  values  of  a  as  plotted  fron  t’ie  data  shown  in  Figs. 

14  to  16.  The  points  corresponding  to  zero  current  density  are  adapted 
from  the  characteristic  current-voltage  curves  measured  using  the  test 
plate  as  plane  electrode.  The  latter  curves  are  shown  in  Fig.  18  for 
reference.  Upon  comparison  of  Figs.  17  and  18,  it  is  clear  that,  for 
a  given  d,  the  total  current  increases  more  rapidly  than  the  maximum 
current  density  as  voltage  increases. 

The  data  as  presented  above  provide  only  qualitative  information 
as  to  the  influences  upon  the  current  density  distribution  of  the 
applied  voltage  and  the  distance  between  the  wire  and  plane  electrodes. 
If  it  can  be  shown  that  a  similarity  exists  among  the  curves  presented 
in  Figs.  14  to  l£>,  the  current  density  distribution  for  any  V  and  d 
can  be  readily  deduced  from  the  foregoing  data,  provided  there  is  corona 
and  the  characteristic  values  of  J  and  x  are  known  for  the  values  of  V 
and  d  of  interest.  By  similarity  it  is  meant  here  that  the  current 
density  distributions  for  two  sets  of  value  of  V  and  d  differ  only  by  a 
scale  factor  in  J  and  x. 

In  an  attempt  to  seek  the  possible  similarity  between  the  curves 
shown  in  Figs.  14  and  1 £,  the  current  density  will  be  normalized  by 
Jmnv.  the  maximum  current  density.  Furthermore,  the  abscissa  variable 
x  will  be  made  dimensionless  by  X1/2,  which  corresponds  to  the  x-value 
where  J  =  i^max*  ®ie  dimensionless  current  density  distributions  so 
constructed,  Figs.  19  to  21,  show  that  the  curves  in  each  of  Figs.  14 
to  1j6  can  be  fairly  well  represented  by  a  single  curve.  From  an  inspec¬ 
tion  of  the  figures ,  it  is  seen  that  a  small  deviation  from  the  tenta¬ 
tive  similarity  profile  occurs  in  the  outer  region  of  the  curves  in  same 
occasions.  Hie  actual  reason  for  this  deviation  is  not  clear. 
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Fig*  17  Maxima  currant  density  on  the  plate  ts  voltage 


with  0,006  wira  at  eantarlina  along  9- inch 
diMntion  of  a  6x9- Inch  flat  plato 


Fig.  19  Dimensionless  currant  density  distribution,  d  •  0.5  inch 
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Maanaionlass  currant  danaity  distribution,  d  «  0.75  inch 


Fig*  21  Dimensionless  current  density  distribution,  dal  inch 


When  Figs.  19  to  21  are  plotted  together,  there  is  obtained  Fig.  22. 
Figures  20  and  21  are  seen  in  this  plot  to  coincide  with  each  other 
fairly  well.  However,  Fig.  19  does  not  seem  to  share  the  same  trend. 
Considerable  deviation  is  noted  again  in  the  outer  regions.  Despite 
this  disagreement,  a  general  conclusion  is  that  the  majority  of  the 
data  is  in  favor  of  the  existence  of  a  similarity  in  the  current  density 
distribution. 


6.3*2  Velocity  Profiles  in  a  Jet  of 
Electric  Wind 

Figures  23  to  25  show  the  velocity  profiles  of  an  electric  wind 
generated  by  a  wire-screen  corona  toj  three  distances  of  the  wire  from 
the  screen:  d  a  0.5,  O.75  and  1  inch.  Instead  of  voltage  as  in  current 
density  distribution,  the  total  current  has  been  used  as  a  parameter 
in  obtaining  the  curves  in  each  figure.  The  abscissa  in  each  case 
designate  as  in  the  previous  section  the  distances  from  the  centerline 
of  the  screen  against  which  the  corona  wire  was  placed. 

Because  of  the  fact  that  the  traversing  mechanism  on  which  the 
total  pressure  probe  was  mounted  can  only  traverse  a  distance  of  1-3/4 
inches,  measurement  of  the  wind  velocity  was  limited  to  a  region  within 
3/4  inch  from  the  centerline.  Beyond  this  range  the  data  curves  were 
extrapolated  (in  accord  with  the  measured  data  in  the  central  portion 
of  the  wind).  In  this  way  the  data  curves  as  shown  in  Figs.  23  to  25 
may  be  subject  to  inaccuracy  in  the  outer  regions. 

Ihe  presence  of  the  screen  wires  undoubtedly  affects  the  passage 
of  the  air  stream.  As  a  matter  of  fact  in  the  course  of  measurements 
it  was  found  that  the  data  point  measured  behind  a  screen  wire  did  not 
fall  on  the  same  track  as  its  neighboring  points.  Such  points  were 
omitted  in  plotting  the  velocity  profiles. 

In  some  cases  the  data  points  are  scattered,  especially  near  the 
outer  edges  of  the  flow.  In  the  latter  range,  the  scattering  of  data 
may  be  due  to  errors  incurred  by  a  yawed  flow  relative  to  the  probe  or 
to  the  decreasing  accuracy  of  the  measurements  at  the  very  low  velocities 
encountered. 

Despite  these  possible  errors,  Figs.  23  to  25  do  provide  some  use¬ 
ful  information  in  relation  to  the  velocity  distribution  of  the  electric 
wind.  A  comparison  of  the  velocity  profiles,  Figs.  23  to  25,  and  the 
current  density  distributions,  Figs.  14  to  16,  exhibits  that  these  two 
series  of  data  are  closely  related  to  each  other.  Ihe  similarity 
between  the  velocity  field  and  the  current  density  field  may  be  sunma- 
rized  as  follows: 

1.  The  velocity  profiles  exhibit  the  same  bell-shaped  dis¬ 
tribution  as  the  current  density. 


Fig.  22  Diasnaionlssa  eurrant  dsnaity  distribution 
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Velocity  distribution  of.  slnotrie 


2.  The  velocity  field  is  confined  to  a  region  about  1  inch 
from  the  centerline,  just  like  the  current  density  dis¬ 
tribution. 

3.  Increase  in  d,  the  distance  between  the  wire  and  the 
screen  tends  to  widen  velocity  field.  Similar  observa¬ 
tion  has  been  made  in  connection  with  the  current  den¬ 
sity  distribution.  One  thing  inconclusive  is  whether 
increase  in  total  current  will  affect  the  area  velocity 
field  when  d  is  kept  constant.  While  the  current  den¬ 
sity  distribution  data  tend  to  suggest  that  the  current- 
area  remains  unaffected  for  constant  d  when  voltage  and 
hence  total  current  is  increased,  extrapolation  in  the 
velocity  profiles  do  not  seem  to  encourage  a  similar 
conclusion. 

4.  In  parallel  with  the  case  of  current  density  distribu¬ 
tion,  an  inspection  of  the  velocity  profile  shows  that 
x,Af  the  value  of  x  for  which  the  velocity  is  half  the 
maximum  velocity  Vmt.  is  almost  constant  for  constant 
d  but  Increases  with  d. 

To  further  the  comparison  between  the  velocity  and  current  density 
fields,  the  velocity  profiles  are  replotted  in  dimensionless  form.  In 
analogy  with  dimensionless  current  density  distribution,  velocities  are 
normalized  by  Vmax,  and  the  coordinate  x  is  made  dimensionless  by  x%. 
Figures  26  to  28  show  the  dimensionless  velocity  profiles  as  transformed 
from  Figs.  23  to  25,  respectively.  Evidently  each  of  Figs.  23  to  25 
can  be  plotted  into  a  single  curve.  The  existence  of  the  similarity  of 
velocity  profiles  of  electric  wind  is  further  evidenced  by  Fig.  29, 
which  plots  together  the  three  curves  shown  in  Figs.  26  to  28. 

It  is  of  interest  to  compare  Figs.  22,  current  distribution,  and 
29  showing  velocity  distribution.  When  these  figures  are  plotted  toge 
together,  Fig.  30  results.  Clearly,  a  similarity  between  the  current 
density  distribution  and  the  velocity  profiles  is  convincingly  indicated 
in  this  plot. 

Figure  31  is  a  plot  of  the  variation  of  the  maximum  velocity  with 
the  total  current.  It  is  noted  from  this  plot  that  the  ma-yimiTm  velocity 
in  the  electric  wind  is  mainly  a  function  of  the  total  current  and  is 
almost  independent  of  d.  When  the  maximum  velocity  is  replotted  using 
the  square  root  of  the  current  as  abscissa  (Fig.  32),  the  data  may  be 
approximated  by  a  straight  line  which  can  be  expressed  as 


0.557^ 


(6.1) 


where  Vm*.*  is  in  ft/sec  and  I  is  in  i*a.  A  similar  relation  has  been 
derived  theoretically  in  Chapter  III  for  the  one -dimensional  model  of 
electric  wind.  For  the  two-dimensional  case,  an  analogous  equation 
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Current  (,/JLm) 


Fig*  31  Maximus  Telocity  of  electrio  wind  to  total  current 


can  also  be  obtained  if  the  relation  that  the  current  density  is  pro¬ 
portional  to  the  total  current  is  used  (Eq.  (3*28). 


6.3.3  Velocity  Profiles  in  the  Wall  Jet 

Figures  33  and  34  show  the  velocity  profiles  for  two  values  of 
d,  0.5  and  0.75  inch  respectively,  and  for  three  or  four  different 
total  currents  in  each  case.  Measurements  were  made  on  the  test  plate 
without  frost.  The  profiles  represent  the  velocity  distributions  along 
the  normal  to  the  plate  at  a  location  2  Inches  below  the  horizontal 
centerline  along  the  vortical  centerline  of  the  plate.  The  corona  wire 
was  horizontal  and  at  the  horizontal  centerline  of  the  plate.  The  test 
data  were  evaluated  under  the  assimptlon  that  static  pressure  is  inde¬ 
pendent  of  the  distance  from  the  plate  and  is  equal  to  the  ambient 
atmospheric  pressure.  A  justification  of  this  assumption  is  evidenced 
if  it  is  recalled  that  the  current  of  and  the  velocity  field  produced, 
by  a  wire-plane  corona  were  confined  primarily  to  an  area  approximately 
one  inch  form  the  centerline  at  which  the  corona  wire  was  positioned. 
Outside  this  area  the  electric  field  and  stagnation  point  effects  on 
the  static  pressure  distribution  are  practically  negligible  and  the 
ambient  pressure  stay  be  considered  to  impress  on  the  flow. 

An  Inspection  of  the  r^ocity  profiles.  Figs.  33  and  34,  indicates 
that  the  velocity  Increases  rapidly  with  y  from  zero  at  the  wall, 
reaches  its  maxlssm  value  at  approximately  y  «  0.08  to  0.1  inch,  and 
then  fails  gradually  to  null  at  y  *  0*5  to  0.55  inch.  Apparently  the 
profiles  are  of  wall  jet  type.  However,  it  is  not  obvious  whether  they 
are  lemlnar  or  turbulent.  A  Pitot  tube  can  only  meausre  the  mean  velo¬ 
city  and  is  incapable  of  detecting  the  turbulent  fluctuations.  A  pre¬ 
liminary  determination  of  the  actual  flow  pattern  (turbulent  or  laminar) 
has  to  be  judged  from  the  shape  of  the  velocity  profiles.  A  closer 
exsari  nation  of  the  foregoing  profile*  suggests  that  at  higher  currents 
(for  example  I  ■  150  m*  and  200  M*  for  d  -  0.75  inch  and  I  ■  200  ua 
for  d  ■  0.5  inch)  the  flow  is  very  likely  turbulent.  The  main  reasons 
for  this  prediction  are  that  the  maxi  mm  velocity  appears  at  lower 
value  of  y  as  I  increases  and  that  the  slope  of  the  velocity  profiles 
near  the  wall  is  steeper  for  higher  current  than  for  lower  current. 

At  low  current  the  velocity  profile  may  represent  either  im-timy  or 
more  likely  transition  flow.  When  the  flow  is  turbulent  or  laminar, 
similarity  exists  and  the  velocity  profiles  should  have  the  same  shape 
when  plotted  non-dimana  tonally.  On  the  other  hand,  if  the  data  were 
measured  in  transition  region  between  stagnation  and  fully  developed 
wall  jet  flows,  no  such  similarity  is  expected. 

Figures  35  and  36  are  the  dimensionless  profiles  plotted  from 
Figs.  33  and  34  respectively.  The  velocity  was  normalized  by  ya,  the 
y-value  at  which  u  *  ^Wr.  From  an  inspection  of  these  two  figures  it 
is  seen  that  in  the  major  portion  of  the  outer  layer  there  is  similarity 
of  the  various  profiles,  with  the  exception  at  the  outer  edge  of  the 
flow.  However,  in  the  inner  layer,  the  profiles  are  non-similar,  which 
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may  be  a  manifestation  of  different  flow  patterns  or  caused  by  the 
decreasing  accuracy  of  the  measurements  near  the  vail  (a  discussion  on 
the  measurements  near  the  vail  can  be  found  in  Dean41).  Figure  37  ie 
a  plot  of  the  following  three  profiles:  I  -  150  ua,  d  -  O.75  inch; 

I  -  100  Ma,  d  *  0.75  inch;  I  ■  200  pa,  d  ■  0.5  Inch.  For  caparison 
the  theoretical  velocity  profiles  for  laminar  (Eq.  (5*19))  and  turbu¬ 
lent1"  vail  Jets  are  also  shown  in  Fig.  37.  As  is  evident  from  this 
figure,  the  three  measured  profiles  agree  satisfactorily  with  the 
theoretical  turbulent  profile.  Further  comparison  between  the  low  cur¬ 
rent  profiles  and  the  theoretical  profiles  indicates  that  the  measured 
data  are  more  in  agreement  with  the  theoretical  turbulent  profile.  It 
thus  appears  that  for  higher  current  the  flow  at  two  inches  from  the 
centerline  is  of  the  turbulent  vail  jet  form,  whereas  at  lower  current 
the  same  location  may  fall  between  the  regions  of  transition  flow  and 
of  turbulent  vail  jet  and  the  flow  in  nearly  turbulent. 

6.3.4  Heat  Transfer  Data 

A  general  trend  is  exhibited  by  the  total  heat  transfer  rate 
versus  time  data  given  in  Figs.  «8  to  40,  which  were  obtained  using  as 
emitter  a  horizontal  single  wire  at  the  plate  centerline.  The  total 
heat  transfer  rate  is  seen  to  decrease  with  time,  but  has  the  tendency 
of  approaching  a  constant  value  as  to  Increases  further.  Similar 
observations  have  been  reported  by  others  in  relation  to  frost  forma¬ 
tion  under  different  circumstances.  Sugavara  et  al30,  in  examining  the 
frosting  process  over  a  flat  plate  under  forced  convection  conditions, 
demonstrated  that  the  total  heat  transfer  rate  decreased  with  time  and 
reached  quasi-steady  state  in  60  to  80  minutes .  In  a  study  of  frost 
formation  on  a  cylinder  surface  in  a  humid  air  cross  flow,  Chung  and 
Algren  ■  also  observed  that  the  heat  transfer  became  quasi-steady  In 
60  to  100  minutes.  Similar  trend  was  also  noted  by  Whitehurst40  when 
frost  grew  on  a  vertical  plate  under  free  convection  conditions,  but 
no  indication  was  given  as  to  the  time  when  the  heat  transfer  essen¬ 
tially  becwe  quasi-steady.  Because  each  test  run  lasted  only  30 
minutes  in  the  present  Investigation,  it  is  difficult  to  estimate  how 
long  it  actually  takes  for  the  heat  transfer  process  to  reach  quasi- 
steady  state.  In  view  of  the  above  mentioned  evidences,  it  may  gener¬ 
ally  be  concluded  that  the  Impingement  of  an  electric  wind  does  not 
change  the  characteristics  of  the  total  heat  transfer  rate  versus  time 
relation  in  connection  with  frost  formation. 

It  nay  be  noted  that  a  similar  trend  is  also  Indicated  by  the 
analysis  developed  In  Chapter  n.  If  Eq.  (2.44)  is  evaluated  at  y  ■  0, 
the  relation  between  q  and  t  (considering  absolute  value  of  q)  can  be 
written 

q  *  ao  -  ait  +  aet2  -  • •  • 


•  • 


40  0  0 
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where  ao,  and  ai  and  &z  are  positive  coefficients.  Because  of  the 
nature  of  the  solution,  the  series  on  the  right-hand  side  is  uniformly 
convergent.  When  t  is  small  q  decreases  almost  linearly  with  t,  as  t 
increases  the  rate  of  change  of  q  decreases.  As  t  tends  to  infinity, 
q  approaches  a  constant  value.  This  qualitative  description  is  in 
agreement  with  the  experimental  results. 

Figure  38,  shows  the  data  obtained  for  four  different  currents  but 
constant  d.  Although  an  inspection  of  the  data  indicates  that  increase 
in  current  results  in  higher  heat  transfer,  the  effects  is  not  signifi¬ 
cant.  Another  set  of  data  given  in  Fig.  39  shows  a  considerable  enhance¬ 
ment  of  heat  transfer  by  the  field.  With  free  convection  conditions 
(no  filed)  as  reference  state,  an  increase  of  100 %  in  total  heat  trans¬ 
fer  rate  is  obtained  at  I  *=  15 0  ua.  The  increase  in  heat  transfer  at 
I  ■  150  ua  over  the  case  where  1  =  100  m&  is  about  20^,  a  much  more 
significant  increase  than  indicated  by  Fig.  38.  While  it  is  uncertain 
what  actually  caused  the  difference  between  the  two  sets  of  data,  the 
facts  that  the  test  conditions  differ  from  one  test  run  to  another  and 
that  the  theoretical  considerations  indicate  (Chapter  VII )  the  heat 
transfer  process  is  less  sensitive  to  the  variations  in  total  current 
as  oampared  to  the  mass  transfer  process  may  have  contributed  to  this 
difference. 

Figure  40  is  a  plot  which  illustrates  the  data  for  two  different 
values  of  d  but  the  same  current.  As  is  obviously  seen,  no  difference 
can  be  noted.  It  appears  that  the  entrainment  as  d  increases  is  so 
small  that  it  can  not  produce  noticeable  influence  on  heat  transfer. 

More  data  are  necessary  before  further  comments  can  be  made  in  this 
respect. 

In  each  test  run,  readings  were  taken  as  soon  as  the  field  was 
turned  on.  The  response  of  the  inlet  and  outlet  refrigerant  tempera¬ 
ture  difference  to  the  sudden  increase  of  convection  heat  transfer  did 
not  take  place  instaneously.  It  usually  took  a  couple  of  minutes  for 
the  refrigerant  passing  through  the  copper  tubing  at  the  back  of  the 
plate  to  reach  a  steady  temperature.  Ibis  adjustment  on  the  part  of 
the  refrigerant  is  apparent  in  the  figures.  The  first  data  point  in 
each  run  is  seen  to  fall  below  the  tentative  position  as  extrapolated 
from  the  rest  of  data.  In  the  course  of  each  test  run,  the  refrigerant 
flow  rate  was  increased  a  couple  of  times  through  the  adjustment  of  the 
controlling  valve  in  order  to  maintain  as  constant  a  plate  temperature 
as  possible.  Each  time  the  adjustment  was  made,  it  took  a  few  minutes 
for  the  refrigerant  to  achieve  a  new  equilibrium  condition.  Since  the 
variation  in  time  of  the  heat  transfer  rate  is  not  drastic,  nearly  the 
same  amount  of  heat  was  transfered  to  the  refrigerant  per  unit  time  in 
the  period  of  adjustment.  Hence  at  lower  flow  rate  the  temperature 
increment  in  the  refrigerant  should  be  greater  than  at  higher  flow  rate. 
If  readings  of  the  temperature  increment  were  taken  before  new  equili¬ 
brium  state  has  been  achieved,  the  data  so  obtained  would  overestimate 
the  heat  transfer.  Ibis  is  exactly  what  happens  in  the  figures. 
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As  noted  In  Section  6.2.2,  the  current  was  observed  to  increase 
as  frost  grew  on  the  test  plate.  Figure  4l  shows  some  typical  varia¬ 
tions  of  the  current  time.  In  the  course  of  each  test  run,  the  voltage 
was  found  to  remain  approximately  constant.  Therefore,  it  is  possible 
that  the  frost  may  act  somewhat  like  a  conductor.  As  the  frost  thick¬ 
ness  increases,  the  distance  between  the  *'ire  and  uue  effective  plane 
electrode  is  shortened,  resulting  in  higher  currant . 


6.3.5  Mass  Transfer  Data 

The  results  of  mass  transfer  measurement  are  shown  in  Fig.  42 
through  48,  in  which  the  total  frost  accumulated  in  a  period  of  five 
minutes  is  plotted  against  the  current.  In  all  the  figures  the  point 
corresponding  to  zero  current  represents  mass  transfer  under  free  con¬ 
vection  conditions.  With  the  exception  of  Fig.  48,  the  data  were 
obtained  using  a  single  wire  electrode.  The  wire  was  horizontal  and 
at  the  centerline  of  the  plate,  the  centerline  being  parallel  to  the 
9-inch  edges  of  the  plate.  The  test  conditions  for  each  curve  are 
specified  by  the  test  plate  temperature  I* ,  the  specific  humidity  wa 
and  temperature  T«,  of  the  ambient  air,  and  the  distance  d  of  the  wire 
from  the  plate.  A  needle  point  was  used  as  discharge  electrode  in 
obtaining  Fig.  48. 

Owing  to  the  limited  capacity  of  the  test  setup,  only  three  or 
four  data  points  were  obtained  for  each  curve.  As  a  result  the  accuracy 
of  the  curves  plotted  on  the  basis  of  the  scare  data  points  available 
is  difficult  to  estimate.  Despite  this  disadvantage,  the  figures  do 
provide  a  general  qualitative  information  on  the  influence  of  electric 
wind  upon  frost  formation.  The  figures  unanimously  indicate  that  the 
frost  deposition  can  be  considerably  enhanced  through  the  use  of  the 
corona  discharge.  An  increase  of  over  200^  can  be  achieved  with  a 
corona  current  of  the  order  of  100  ua.  A  further  inspection  of  the 
curves  shows  that  the  rate  of  increase  of  the  frost  deposition  decreases 
with  current.  Indicative  of  the  fact  that  the  mass  transfer  will  asymp¬ 
totically  approach  a  miximum  value  as  the  current  tends  to  infinity. 

A  physically  more  reasonable  statement  of  this  is  that  beyond  certain 
value  of  the  current  no  further  appreciable  Increase  in  mass  transfer 
may  be  expected.  The  similarity  in  shape  of  the  curves  for  wire-plane 
and  point-plane  coronas  suggests  that  the  same  conclusions  also  hold 
true  for  point-plane  corona,  an  evidence  of  the  existence  of  similarity 
between  the  axially  synmetric  and  two-dimensional  flow  fields  created 
by  point-plane  and  wire-plane  coronas  repectively. 

In  Figs.  49  and  50  respectively  are  plotted  together  Figs.  42  to 
44  and  45  to  47  by  shifting  the  ordinate,  when  necessary,  such  that  the 
points  at  zero  current  coincide.  It  is  seen  in  these  figures  that  the 
data  for  different  values  of  d  can  be  fitted  approximately  by  a  single 
curve. 
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Total  frost  weight  ts  currant  for 
period  of  5  ninutee— wire  corona 
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In  Section  1.2  the  increase  of  frost  formation  under  the  action  of 
corona  discharges  has  been  attributed  to  the  impingement  of  the  electric 
wind  and  the  total  mass  transfer.  For  convenience,  let  the  velocity 
vmax  the  location  of  the  stagnation  point,  in  the  absence  of  a  solid 
plane  be  taken  as  the  characteristic  velocity  of  the  electric  wind.  A 
comparison  between  Figs.  31  and  5 0  sho  s  that  the  total  current  I  has 
quite  similar  effects  on  v  and  the  total  mass  transfer,  thus  substantiat¬ 
ing  the  existence  of  a  possible  relation  between  the  latter  two  quantities 
The  complexity  of  the  flow  field  of  an  impinging  electric  wind  seems  to 
suggest  that  such  a  relation  is  not  a  simple  one.  That  this  is  the  case 
is  further  indicated  in  Fig.  51>  which  plots  the  total  mass  transfer 
versus  v^y  on  the  basis  of  Figs.  31  and  50. 

An  explanation  of  the  behavioral  characteristics  of  the  total  mass 
transfer  data  curves  in  terms  of  I  goes  as  follows.  At  zero  current, 
there  is  no  electric  wind  and  the  mass  transfer  process  is  of  free  con¬ 
vection  form  in  nature.  When  I  differ  from  zero,  an  induced  electric 
wind  changes  the  transport  mechanism  from  free  to  forced  convection, 
resulting  in  escalation  of  the  frost  accumulation.  Since  the  total  mass 
transfer  increases  with  Vnmv  (Fig.  51)  and  v^y  increases  with  I  (Fig.  31) 
it  follows  that  the  total  mass  transfer  increases  with  I.  It  has  been 
shown  previously  that  Fig.  31  can  be  approximated  by  Eq.  (6.1).  A  dif¬ 
ferentiation  shows  that  the  rate  of  change  of  Vmay  with  respect  to  I, 
dVjjjax/dl,  is  inversely  proportional  to  ^t.  As  I  increase,  dvnwy/dl 
becomes  smaller  and  smaller.  Correspondingly,  the  rate  of  change  of  the 
total  mass  transfer  with  respect  to  I  decreases  with  I.  This  accounts 
for  the  fact  that  all  of  the  curves  shown  in  Figs.  42  through  48  tend 
to  level  off  as  I  increases. 

It  remains  to  explain  how  variation  in  vmy  affects  the  total  mass 
transfer.  To  this  end,  it  is  helpful  to  recourse  to  the  theoretical 
considerations  of  Chapter  IV  and  V.  Eqs.  (4.25)  and  (4.62)  indicate 
that  the  mass  transfer  rate  per  unit  surface  area  is  proportional  to 
</ui  in  the  stagnation  flow  regime.  The  parameter  m  is  related  to  the 
velocity  and  the  width  of  the  Impinging  stream.  For  a  uniform  jet  of 
velocity  vi  and  half  width  bj,  m  has  been  shown  to  take  the  value 
ui  =  vjit/4bj.  To  apply  thisd expression  to  the  case  of  two-dimensional 
electric  wind,  an  acceptable  selection  for  Vj  and  b-j  is  v^y  and  xvi 
respectively.  Thus,  near  the  stagnation  point,  the  mass  transfer  rate 
is  proportional  to  Jv^y.  In  the  wall  Jet  region,  Eqs.  (5.47)  and  (5.80) 
show  that  the  mass  transfer  rate  is  determined  by  two  parameters  R  and 
4.  According  to  the  definition  of  R,  it  can  b«t  deduced  that  R  is  pro¬ 
portional  to  vmy3.  If  4  is  determined  in  the  manner  suggested  in 
Section  5*3,  it  is  readily  seen  that  4  and  v^-^  are  related  linearly. 

Using  these  above  relations  between  R  and  vwlPy .  and  l  and  v^y  in 
Eq.  (5*47)  or  (5.80),  it  is  inconclusive  as  to  how  the  mass  transfer 
rate  in  a  wall  jet  varies  with  Ymav.  On  the  basis  of  physical  reason¬ 
ing,  it  is  anticipated  that  increase  in  Vm»y  will  result  in  higher  mass 
transfer  rate  in  a  wall  jet.  Bearing  this  in  mind,  it  is  found  that 
the  special  case  when  4  is  negligible  as  compared  to  x,  affords  an 


acceptable  relation  between  the  mass  transfer  coefficient  and  vTnftV| 

With  R  a  Vmfl/y3.  hjQ  is  thus  proportional  to  Vp^y3/4.  Further  discussion 
in  this  connection  will  be  taken  up  in  the  next  chapter. 

While  Figs.  49  and  50  tend  to  indicate  that  the  total  mass  transfer 
is  dominantly  controlled  by  the  total  current  and  is  not  significantly 
affected  by  d,  another  test  series  aimed  to  investigate  the  possible 
variation  of  frost  accumulation  with  d  at  fixed  current  showed  that 
such  variation  did  occur.  Typical  data  of  such  tests  are  shown  in 
Fig.  52.  A  general  trend  is  that  the  frost  accumulation  tends  to 
increase  with  d,  I  being  maintained  constant.  It  is  known  in  two-dimen¬ 
sional  laminar  free  jet  problem  that  the  volume  rate  of  discharging  per 
unit  height  of  slit  is  proportional  to  z1'3,  where  z  is  the  distance 
frcm  the  slit  along  the  axis  of  the  jet.  nils  relation  is  also  plotted 
in  Fig.  52  for  comparison.  As  can  be  seen  from  this  figure,  frost 
accumulation  increases  with  d  in  a  similar  manner  as  does  the  volume 
rate  of  flow  of  a  free  jet  with  z.  It  is  thus  deduced  that  an  electric 
wind  might  behave  like  a  free  jet,  as  the  air  stream  blows  along  fluid 
particles  from  the  surroundings  are  carried  away  with  the  wind  owing  to 
friction  on  its  boundary.  Consequently,  increase  in  d  results  in  more 
air  being  discharged  against  the  plane  electrode  and  hence  more  frost 
accumulation.  How  this  last  step  actually  takes  place  requires  further 
explanation.  In  6.3.2,  Vm»Y  has  been  found  to  remain  approximately 
constant  and  x>A  has  been  noted  to  increase  as  d  increases,  an  indication 
that  the  volume  rate  of  flow  increases  with  d.  This  means  the  ux  will 
decrease  as  d  increases  and  less  frost  will  deposit  per  unit  area  in 
the  stagnation  region.  On  the  other  hand,  it  has  been  demonstrated 
that  increase  in  d  inclines  to  widen  the  range  of  the  wind  and  hence 
the  area  of  the  stagnation  flow.  The  net  results  of  these  two  contra¬ 
dictory  effects  of  d  on  mass  transfer  is  expected  to  be  in  favor  of 
increasing  frost  deposit.  In  the  wall  jet  region,  increase  in  flow 
rate  definitely  will  increase  mass  transfer  as  is  obvious  frcm  the 
definition  of  R.  It  may  be  noted  in  passing  that  in  retrospect  the 
data  on  electric  wind  and  wall  jet  indicate  that  the  flow  rate  does 
increase  with  d.  Therefore,  the  speculation  of  the  entrainment  phenome¬ 
non  in  an  electric  wind  is  further  substantiated. 
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CHAPTER  VII 


COMPARISON  OF  THEORY  AND  EXPERIMENT 


7.1  METHODS  OF  CALCULATION 

Efforts  will  be  made  in  this  chapter  to  attempt  a  comparison  of 
theoretical  and  experimental  results.  Such  a  comparison  in  the  present 
study  presents  considerable  difficulty  owing  to  a  number  of  reasons. 
First,  no  satisfactory  theory  exists  in  the  literature  as  regards  the 
phenomenon  of  the  electric  wind.  Second,  the  electric  wind  itself  is 
highly  non-uniform  and  the  stagnation  flow  associated  with  its  normal 
impingement  on  a  plane  surface  has,  as  yet,  not  been  solved.  Hie 
theoretical  consideration  of  Chapter  IV  have  been  based  on  one-dimen¬ 
sional  model  of  electric  wind  and  the  results  are  strictly  applicable 
only  for  a  region  in  the  neighborhood  of  the  stagnation  point  small  in 
comparison  with  the  width  of  the  oncoming  stream.  Third,  the  aerodynamic 
processes  of  transition  from  stagnation  flow  to  fully  developed  wall  jet 
have  not  yet  been  explored.  Our  knowledge  concerning  transport  mecha¬ 
nisms  in  this  region  is  practically  void.  Fourth,  experimental  evidences 
indicate  turbulent  flow  prevails  in  the  wall  jet.  Lack  of  adequate 
theory  on  transport  processes  in  turbulent  wall  jet  has  forced  one  to 
consider  laminar  wall  jet  instead,  which  certainly  can  not  provide  a 
true  picture.  Fifth,  the  frost  formation  is  a  transient  process,  but 
the  theoretical  considerations  have  been  based  upon  quasi-steady  con¬ 
ditions.  Finally,  in  view  of  the  complicated  flew  field  involved,  it 
is  most  desirable  to  investigate  local  transport  coefficients.  Unfortu¬ 
nately,  in  the  present  investigation  this  was  not  accomplished.  Owing 
to  these  above  disadvantages  that  face  the  problem  of  interest,  it  is 
felt  that  the  significance  of  the  present  chapter  lies  primarily  on 
correct  theoretical  predictions  of  the  qualitative  trend  of  the  varia¬ 
tion  of  transport  processes  with  current. 

In  order  to  canpare  the  effectiveness  of  an  electric  wind  in 
inproving  the  transport  processes  relative  to  free  convection,  the  heat 
and  mass  transfer  coefficients  in  free  convection,  two-dimensional 
stagnation  flow  and  two-dimensional  laminar  wall  jet  will  be  evaluated 
using  the  following  values: 

Plate  surface  temperature ,  =  10 °F, 

Free  stream  temperature,  Too  =  70 °F, 

Plate  surface  mass  fraction  of  water  vapor,  wiw  =  0.0013, 

Free  stream  mass  fraction  of  water  viper,  wi«  =  0.01. 

The  calculations  will  be  based  upon  the  theoretical  expressions  for 
zero  normal  velocity  at  the  frost  surface. 
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and 


0.677s/Pr(Gr  +  |£  Gr’)^ 
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(0.952  +  Pr)£  L 
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where 

Gr  •  -  is  the  average  Grashof  number, 

v2 


Gr*  *  is  average  Grashof  number  for  mass  transfer, 


Pt  ■  -  is  the  coefficient  of  thermal  expansion, 
T 


(Ma  -  Hi)/Mi 
v^)  *  1 


is  the  concentration  coefficient  of  expansion 


and 

Mi,  M2  *  molecular  weight  of  water  vapor  and  dry  air,  respectively. 
All  properties  of  the  air  will  be  evaluated  at  the  mean  film  tempera¬ 
ture,  Tn  ■  40°F.  At  this  temperature,  the  values  of  k,  D,  v,  Pr,  Sc, 
gpt/v2  and  gfr/v2  are: 


V 

-  1.47  X 

I0"4ft2/sec, 

k 

-  0.0142 

Btu/hr-ft-*F 

D 

■  2.79  X 

I0“4ft2/sec, 

Pr 

-  0.714, 

Sc 

-  0.528, 

Rgt 

„2 

-  2.98  X 

10eft"3  V1 

v 
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and 

^  =  9.03  x  lcr^r3  . 

V2 

For  L  =  6  inches,  Gr  and  Gr'  are: 

Gr  =  2.2k  x  107  , 

Gr'  =  9.83  x  10s  , 

Since  Gr'  is  small  as  compared  to  Gr,  it  will  be  neglected  in  the  suc¬ 
ceeding  calculations.  Substitution  of  the  foregoing  values  into  Eq.  (7*i) 
and  (7.2)  gives  the  following  values  of  average  heat  and  mass  transfer 
coefficients  for  free  convection: 

h  *  64k 
=  50D 

Proceeding  next  to  the  calculations  of  the  heat  and  mass  transfer 
coefficients  at  the  stagnation  point  of  an  impinging  two-dimensional 
electric  wind,  Eqs.  (4.24)  and  (4.26)  will  be  evaluated,  namely, 

h  ’ 2  k  (  )*  (p*-)1/3  (7-3) 

2  '7.052V/ 

and 

\  »§  d(-HJ-)^(Sc)1/3  (7.4) 

m  2  V  7.052v  / 

The  chief  concern  here  is  the  determination  of  the  parameter  ui.  As 
noted  earlier  (Section  6.3«5)j  ui  is  related  to  the  velocity  vj  and  the 
half  width  bj  of  a  uniform  jet  by  ui  =  nvj/4bj.  In  order  to  utilize 
this  relation  to  estimate  the  value  of  ux  for  a  two-dimensional  non- 
uniform  electric  wind,  certain  modification  has  to  be  excerised.  It 
has  been  suggested  in  the  preceding  chapter  to  replace  Vj  and  bj  by 
Vnutv  and  x%,  respectively,  of  the  actual  impinging  wind.  This  same 
procedure  will  again  be  adopted  here  in  carrying  out  the  computations. 

To  be  specific,  the  case  where  v^y  =7-8  ft/sec  and  x>A  ■  0.28  inch, 
corresponding  to  I  =  200  ua  and  d  =  0.5  inch,  will  be  used  in  a  sample 
calculation.  For  this  particular  case  ui  is  found  to  be  ux  =  262  sec"1 
Substituting  this  value  of  ui  and  the  required  properties  of  the  air  at 
a  mean  temperature  of  40°F  into  Eqs.  (7.3)  and  (7.4),  h  and  h_  are  found 
to  be  h  *  675k  and  =  6l2D. 

Evaluation  of  the  heat  and  mass  transfer  coefficients  for  laminar 
wall  jet  will  be  based  upon  Eqs.  (5 .4l)  and  (5.48)  with  0x(O)  =  0.273 
and  cp’(0)  =  O.258,  namely, 
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h(x)  =  0.273k(5R/32v3)^(x+ir3/4 


(7.i>) 


and 

hjx)  =  0.258D(5R/32v3)^(x+i)“3/4  (7.6) 

The  main  task  now  is  the  determination  of  the  two  constants  R  and  £. 
Arguing  that  the  volume  flow  rate  as  well  as  the  maximum  velocity  of 
the  impinging  jet  remains  practically  unchanged  after  being  deflected 
by  the  plane  surface,  Glauert37*30  has  suggested  to  estimate  R  and  l 
from  the  conditions  of  the  impinging  stream.  From  the  definition  of 
R  (that  is  Eq.  (5«7))j  be  has  deduced  Rq.  (5.22)  for  R,  namely, 

R  =  £  (typical  velocity) (volume  flow/unit  length)2  (7*7) 

The  constant  length  l  may  be  obtained  by  equating  the  maximum  wall  jet 
velocity  at  x=0  as  given  by  Eq.  (5 .21)  to  the  maximum  initial  velocity, 
namely, 

/  =  (5R/2v)(0.315/vnax)2  (7.8) 

where  0.315  is  the  maximum  value  of  f* .  Eqs.  (7-7)  and  (7.8)  have  not 
been  subject  to  experimental  verification  previously,  and  their  validity 
should  not  be  taken  for  granted.  In  the  present  program  two  other 
methods  of  determining  R  and  i  are  also  considered.  The  transport  coef¬ 
ficients  calculated  using  the  three  methods  of  estimating  R  and  i  will 
be  compared  with  the  case  of  free  convection  as  well  as  experimental 
results.  Conclusion  will  be  made  as  to  which  method  yields  the  most 
satisfactory  values  of  h  and 

The  determination  of  R  by  means  of  Eq.  (7*7)  depends  on  a  suitably 
selected  typical  velocity  of  the  impinging  jet  of  air  stream.  Since 
the  volume  flow  rate  per  unit  length  of  the  jet  is  constant,  a  typical 
velocity  of  the  jet  which  will  give  an  average  value  of  R  across  the 
jet  is  the  average  of  the  Jet.  The  flow  rate  per  unit  length  and  the 
average  velocity  of  the  jet  of  electric  wind  are  determined  from  experi¬ 
mentally  measured  velocity  profiles  shown  in  the  foregoing  chapter. 

For  the  case  pf  0  =  200  ua  and  d  =  0.5  inch,  the  flow  rate  per  unit 
length  is  0.225  ft2/sec  and  the  average  velocity  is  2.7  ft/sec.  Using 
these  values  in  Eq.  (7.7) >  it  is  found  R  =  0.0684  ft5/sec3»  The  value 
of  i  determined  by  Eq.  (7*8)  is  1.89  ft.  Hie  heat  and  mass  transfer 
coefficients  at  x=2  inches  for  this  case  are  h  =  4l.7k  and  »  39«7D. 

The  second  method  of  determining  R  and  /  is  similar  to  the  one 
described  above,  except  that  the  wall  jet  data  are  used  in  lieu  of  the 
electric  wind  data.  For  I  =  20C  ua  end  d  =  0.5  inch,  the  flow  rate 
per  unit  length  as  measured  from  Fig.  33  is  0.133  ft2/sec  and  the 
average  velocity  is  2.96  ft/sec.  The  values  of  R  and  l  thus  found  are 
0.0262  ft5/sec3  and  O.818  ft,  respectively.  The  corresponding  transport 
coefficients  are  h  =  52.5k  and  hgj  =  49. 5D. 
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The  third  methos  Is  based  on  the  assumption  that  l  is  small  as 
compared  to  x  ■  2  Inches  and  can  be  neglected.  And  R  Is  obtained  by 
substituting  the  measured  «***•<  «»">  wall  jet  velocity  In  the  second  of 
Eq.  (5.21),  putting  x  =  0.167  ft  (that  is  x  =  2  inches)  and  f*  =  0.315> 
and  solving  for  R,  namely, 

R  «  (7-9) 

For  the  example  considered  above,  R  =  0.00445  ft5/sec3>  h  =  128k  and 

hm  -  1211). 

With  the  transport  coefficients  for  stagnation  and  wall  jet  flows 
determined  as  above,  the  total  heat  and  mass  transfer  rates  will  now  be 
estimated.  Since  no  information  is  available  concerning  the  transport 
processes  in  the  transition  region,  this  region  will  be  neglected.  It 
will  be  assumed  that  the  stagnation  flow  covers  an  area  of  4xA  and  the 
average  transport  coefficients  in  this  region  is  one  half  the  sum  of 
the  values  at  the  stagnation  point  and  at  x  ■  2  inches  for  all  jet,  the 
latter  are  considered  as  the  average  transport  coefficients  outside  the 
stagnation  region. 


7.2  DISCUSSION  OF  THEORETICAL  PREDICTIONS 

The  heat  and  mass  transfer  coefficients  at  the  two-dimensional 
stagnation  point  ore  listed  in  Table  1.  As  expected  h  and  increases 
with  I  for  constant  d  and  decreases  with  increasing  d  for  constant  I. 

A  discussion  of  this  has  been  given  in  Section  6.3.5  and  no  further 
conments  will  be  made. 


Table  1  -  Transport  coefficients  at  two-dimensional  stagnation  point 
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For  the  two-dimensional  laminar  wall  jet,  the  local  heat  and  mass 
transfer  coefficients  evaluated  using  the  three  methods  show  a  certain 
peculiar  behavior.  Tables  2  and  3  give  the  values  of  h  and  deter¬ 
mined  by  using  the  electric  wind  and  the  wall  jet  data,  respectively. 

An  inspection  of  these  tables  indicated  that  no  unique  relation  exists 
between  the  transport  coefficients  and  the  total  current.  An  increase  in 
I  does  not  always  result  in  higher  values  of  h  and  1^,  in  constrast  to 
expectation.  It  is  all  the  more  surprising  to  also  note  that  the  trans¬ 
port  coefficients  for  wall  jets  are  at  most  of  the  same  magnitude  as 
those  for  free  convection.  According  to  the  theory  on  free  convection28* 
the  maximum  velocity  at  a  position  x  from  the  upper  edge  of  the  vertical 
plate  is  given  by 

w.3v[sei^]\» 

At  x  ■  6  inches,  iwy  =  0.272  ft/sec,  which  is  the  maximum  attainable 
velocity  along  the  test  plate  under  free  convection  conditions  and  is 
much  smaller  than  the  measured  mj^Hminn  velocity  in  the  wall  jet.  It 
is  rather  mysterious  that  higher  velocity  does  not  have  the  effect  of 
enhancing  transport  processes.  Furthermore  this  theoretical  result 
contradicts  the  experimental  finding  that  frost  was  found  to  form  more 
rapidly  even  in  regions  far  away  from  the  stagnation  region  (i.e.  in 
the  wall  jet  region)  under  the  action  of  electric  wind  than  in  the 
absence  of  field.  Apparently,  Eqs.  (7*7)  and  (7.8)  do  not  seem  to 
provide  correct  information  for  the  wall  Jet. 
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Table  2  -  Transport  coefficients  for  two-dimensional  laminar 
wall  jet  —  based  on  electric  wind  data 


n 
1 1 


I 

(ua) 

d 

(in) 

vmax 

(ft/sec) 

V.F.R.W 

(ft2/sec) 

A.V.W 
(ft /sec) 

R 

(ft*/s«c3) 

i 

(ft) 

h/k 

(ft"1) 

hm/D 

(ft’1) 

200 

0.5 

7.8 

0.225 

2.7 

0.0684 

1.89 

4l 

39.7 

150 

0.5 

7.4 

0.196 

2.35 

0.0451 

1.39 

4j  , 

44.4 

100 

0.5 

6.1 

0.162 

1.95 

0.0255 

1.23 

46.6 

44.1 

50 

0.5 

4.4 

0.121 

1.45 

0.0132 

1.28 

39.3 

37.2 

150 

0.75 

7.2 

0.275 

2.36 

0.0891 

2.9 

31.6 

29.9 

100 

0.75 

6.1 

0.233 

2.15 

0.0584 

2.65 

30.6 

28.9 

50 

0.75 

4.2 

0.154 

1.54 

0.0183 

1.75 

31.2 

29.5 

[, 

1 

I 

0 


0 


u 


0 

B 


(a)  V.F.R.  =  volume  flow  rate 

(b)  A.V.  ■  average  velocity 


Table  3  -  Transport  coefficients  of  two-dimensional  laminar  wall 
jet  —  based  on  wall  Jet  data 


(Ma) 

(in) 

(ft/sec) 

(ft2/sec) 

(ft/sec) 

(ft5/sec5) 

(«) 

(ft-1) 

(ft-1) 

200 

0.5 

6.7 

0.133 

2.96 

0.0262 

0.818 

52.5 

49.5 

150 

0.5 

5.9 

0.112 

2.49 

0.0156 

0.59 

56.1 

53.0 

100 

0.5 

4.7 

0.09 

2.0 

0.008 

0.465 

54.3 

51.3 

50 

0.5 

3.2 

0.057 

1.26 

0.002 

0.167 

62.3 

58.8 

150 

0.75 

6.4 

0.15 

3.34 

0.038 

1.38 

40.8 

38.6 

100 

0.75 

5.6 

0.12 

2.67 

0.019 

0.863 

47.0 

44.4 

50 

0.75 

3.6 

0.078 

1.74 

0.053 

O.528 

45.7 

43.2 

Tabulated  in  Table  4  are  the  transport  coefficients  for  laminar  wall 
Jets  as  predicted  by  the  third  method.  It  is  noted  that  h  and  increase 
with  I  and  are  approximately  twice  the  magnitude  of  the  corresponding 
average  free  convection  transport  coefficients.  Thus,  it  is  concluded 
that  the  values  of  h  and  1%  determined  by  this  method  are  more  reasonable. 


Table  4  -  Transport  coefficients  of  two-dimensional  .laminar  wall 
Jet  —  based  on  Eq.  (7. 9) 


I 

(M*) 

d 

(in) 

Umax 

(ft/sec) 

R 

(ft5/sec3) 

h/k 

(ft"1) 

bm/b 

(n-1) 

200 

0.5 

6.7 

0.00455 

128 

121 

150 

0.5 

5.9 

0.00344 

119.5 

113 

100 

0.5 

4.7 

0.00214 

106 

100 

50 

0.5 

3.2 

0.00101 

89.8 

84.9 

150 

0.75 

6.4 

0.00405 

125 

118 

100 

0.75 

5.6 

0.00311 

116.5 

no 

50 

0.75 

3.6 

0.00128 

93.5 

88.3 

The  total  heat  and  mass  transfer  rates  given  in  Table  5  are  pre¬ 
dicted  on  the  basis  of  Tables  1  and  4.  The  ratio  of  the  predicted 
total  mass  transfer  rate  with  field  and  the  corresponding  value  for 
free  convection  (no  field)  is  plotted  in  Fig.  53  as  function  of  current. 
Since  the  effect  of  free  convection  is  not  considered  in  the  analysis, 
the  predicted  mass  transfer  rate  vanishes  at  zero  current.  For  com¬ 
parison,  Fig.  50  is  also  shown  in  Fig.  53 »  using  the  value  of  total 
mass  transfer  at  zero  current  as  a  reference.  The  theoretical  prediction 
is  observed  to  agree  quite  well  with  expe’  r  wUit.  However,  it  should  be 
cautioned  that  the  present  investigation  is  aimed  only  to  provide  means 
of  comparing  qualitatively  the  simplified  theoretical  considerations 
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and  the  experiment,  and  no  quant  tive  accuracy  is  claimed.  It  is 
worth-while  to  emphasize  that  wh*  the  laminar  wall  jet  theory  tends 
to  underestimate  the  mass  (as  well  as  heat)  transfer  in  the  wall  Jet 
region  (which  is  actually  turbulent),  the  assumption  that  the  situation 
is  quasi-steady  at  the  test  plate  conditions  inclines  to  overestimate 
the  total  mass  (and  heat)  transfer.  Another  point  worth  noting  is  the 
selection  of  the  stagnation  flow  area  and  the  average  transport  coef¬ 
ficient  in  the  wall  jet  region  is  rather  arbitrary,  for  no  information 
on  the  position  of  transition  is  available. 

Table  5  -  Total  heat  and  mass  transfer  rates 


(Ma) 

(in) 

^tOt/PDfrleo^w) 

(ft) 

200 

0.5 

67.1 

62.7 

150 

0.5 

63.4 

59.3 

100 

0.5 

56.8 

52.9 

50 

0.5 

48.3 

44.9 

150 

0.75 

67.3 

62.7 

100 

0.75 

62.7 

58.2 

50 

0.75 

50.9 

47.3 

Free 

convection 

24 

18.75 

A  similar  comparison  between  theory  and  experiment  for  total  heat 
transfer  is  given  in  Fig.  54,  the  experimental  data  being  taken  from 
Fig.  40.  The  correlation  between  theory  and  experiment  is  poor  in  this 
case.  The  discrepancy  may  be  attributed  to  the  fact  that  the  measured 
total  heat  transfer  rate  represents  the  total  heat  gain  to  the  refri¬ 
gerant  through  convection,  radiation  and  conduction.  While  the  theoret¬ 
ical  curve  considers  convective  heat  transfer  only,  the  uncorrected 
data  actually  represent  the  ratio  Q/(Q)free  conv>  where  Q  and 

(Q)free  corrv  8X6  the  total  heat  transfer  rates  with  and  without  field 
respectively.  If  the  radiative  and  conductive  heat  transfer  is  inde¬ 
pendent  of  the  electric  field,  then 


9tot 


(^tot)f: 


ree  conv 


_ Q, 

(Q)free 


conv 


Ifce  conductive  heat  transfer,  which  includes  heat  gain  from  the  environ¬ 
ment  through  the  edged  of  the  test  plate  and  the  insulating  material, 
is  estimated  to  be  small  as  compared  with  the  radiative  heat  transfer. 
After  making  correction  for  radiative  heat  transfer,  the  measured  data 
are  seen  to  be  in  closer  agreement  with  the  theory. 


(  1 I 
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To  further  the  comparison  between  theory  and  experiment,  a  numeri¬ 
cal  value  of  the  magnitude  of  the  total  mass  transfer  will  be  given  and 
compared  with  the  measured  frost  deposit.  For  the  case  I  =  200  ua  and 
d  =  0.5  inch,  the  total  mass  transfer  for  a  period  of  five  minutes  is 

“tot  =  6  2 . 7pD  (wi  oTwiw )  £q 

At  a  mean  temperature  of  4o°F,  the  density  of  the  air  is  0.0795  lbm/fta. 
Using  this  value  of  p  and  the  value  of  (wiw-wiw)  given  in  the  previous 
section,  the  total  mass  transfer  to  the  test  plate  is  calculated  to  be 
1.65  0a.  In  the  case  of  heat  transfer,  the  magnitude  of  the  measured 
convective  heat  transfer  rate  is  found  to  be  approximately  2.5  times 
greater  than  the  value  estimated  by  the  theory.  The  reason  for  this 
discrepancy  is  rather  obscure  and  further  research  is  required  in  order 
to  made  a  more  conclusive  comparison  between  theory  and  experiment. 

On  the  ground  of  the  close  agreement  between  theory  and  experiment 
as  indicated  in  Fig.  53  for  the  total  mass  transfer,  a  number  of  con¬ 
clusions  may  be  drawn.  First,  within  the  velocity  range  of  the  experi¬ 
ments,  the  assumption  of  laminar  flow  in  the  wall  jet  region  provides 
a  reasonable  prediction  on  mass  transfer.  Second,  the  deviation  of  the 
actually  transient  frosting  process  from  the  quasi-steady  assumption 
within  an  interval  of  5  minutes  is  not  significant  under  the  test  con¬ 
ditions.  Finally,  the  method  of  calculation  adopted  in  estimating  the 
total  mass  transfer  is  very  satisfactory.  An  example  of  the  distribu¬ 
tion  of  the  mass  transfer  coefficient  over  the  entire  test  plate  used 
in  the  present  calculation  is  displayed  in  Fig.  55  by  the  dashed  line. 
Judging  from  the  velocity  distribution  of  the  impinging  electric  wind, 
it  is  expected  that  the  transfort  coefficients  may  also  assume  a  bell¬ 
shaped  distribution.  The  predicted  distribution  based  on  the  simplified 
model  is  shown  in  the  same  figure  by  the  curve  made  up  of  alternate 
long  and  short  dashes.  For  comparison,  the  mass  transfer  coefficients 
as  predicted  by  laminar  wall  jet  alone  is  also  indicated  by  Fig.  55  by 
solid  lines. 

Fi  sure  56  shows  the  relative  effectiveness  of  an  electric  wind  in 
enhancing  the  heat  and  mass  transfer  as  predicted  by  the  theory  for 
the  test  plate  under  consideration.  A  general  conclusion  is  an  electric 
wind  is  more  effective  in  increasing  mass  transfer  than  heat  transfer. 
For  example,  at  I  =  200  ua  the  increase  in  mass  transfer  relative  to 
free  convection  conditions  is  233$;  whereas  the  increase  in  heat  trans¬ 
fer  is  only  l80“$.  Of  course  the  above  percentage  increase  are  valid 
only  for  the  particular  test  plate  considered  in  the  present  investiga¬ 
tion. 


CHAPTER  VIII 


CONCLUSIONS  AND  RECOMMENDATIONS 


8.1  SUMMARY  AND  CONCLUSIONS 

An  experimental  program  was  undertaken  to  investigate  the  behavior 
of  frost  formation  in  an  electric  field.  To  facilitate  comparison  of 
the  effectiveness  relative  to  free  convection  conditions  of  an  electric 
field  in  improving  the  heat  and  mass  transfer  processes  associated  with 
frost  formation,  the  test  plate  was  installed  vertically  and  the  experi¬ 
ments  were  carried  out  under  ambient  air  conditions. 

The  experimental  study  included  two  phases.  The  preliminary 
exploratory  tests  were  aimed  to  explore  the  possible  effects  of  electric 
fields  upon  the  frosting  phenomenon.  The  following  conclusions  may  be 
drawn  from  the  results  of  this  test  sequence: 

1.  A  uniform  electric  field  does  not  appear  to  have  any 
observable  effects  on  frost  formation  at  the  initial 
stages  of  the  process. 

2.  A  non-uniform  field  characterized  by  a  corona  discharge 
can  enhance  considerably  the  frost  formation. 

3.  The  effects  of  corona  discharges  on  frosting  can  be 
attributed  to  the  impingement  of  electric  wind  induced 
by  the  corona,  changing  the  flow  field  from  free  con¬ 
vection  to  forced  convection. 

On  the  basis  of  the  preliminary  experimental  findings,  a  second 
experimental  sequence,  using  single  wire  electrode  configuration,  was 
designed  to  study  the  flow  field  generated  by  the  impinging  electric 
wind  and  to  provide  further  information  on  the  electric  wind  effects 
upon  heat  and  mass  transfer  in  connection  with  frost  formation.  The 
mass  transfer  data  indicated  that  at  a  total  current  of  between  100  pa 
and  200  na,  the  total  frost  deposit  for  five  minutes  could  be  increased 
up  to  200$  for  the  particular  test  plate  used  in  this  investigation, 
as  compared  to  the  case  of  no  field.  The  rate  of  increase  in  total 
frost  accumulation  was  found  to  decrease  with  increading  current  and  to 
approach  zero  asymptotically  as  the  current  tended  to  infinity.  The 
heat  transfer  data  also  showed  significant  increase  in  total  heat  trans¬ 
fer  rate  in  the  presence  of  electric  wind  and  an  increase  of  up  to 
100$  at  150  m&  was  measured  for  the  6x9  inches  test  plate.  However, 
the  data  were  insufficient  to  indicate  a  similar  trend  as  exhibited  by 
the  mass  transfer  data. 

From  velocity  measurements  behind  a  wire-screen  collector,  the 
electric  wind  was  found  to  be  confined  to  within  approximately  one  inch 
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from  the  stagnation  point  and  the  velocity  assume  a  bell-shaped  distri¬ 
bution.  Outside  this  region,  the  velocity  was  essentially  zero.  When 
plotted  non-dimens ionally,  the  data  favored  the  existence  of  a  similar¬ 
ity  not  only  among  the  velocity  profiles  but  also  between  the  current 
density  distribution  and  the  velocity  profile  for  various  electric  con¬ 
ditions.  At  a  distance  two  inches  from  the  stagnation  point  along  the 
test  plate,  the  velocity  profile  normal  to  the  plate  was  found  to  be  of 
a  turbulent  wall  jet  type  at  high  current.  At  lower  current  the  profile 
was  observed  to  deviate  slightly  from  the  turbulent  one. 

The  above  velocity  measurements  substantiate  the  postulate  that, 
in  analogy  with  the  impingement  of  a  free  jet,  the  flow  field  produced 
by  an  impinging  electric  wind  may  be  divided  into  four  regions.  Prior 
to  impingement,  the  flow  is  controlled  by  the  electric  field  and  may  be 
called  the  electric  wind  proper.  After  impingement,  the  wind  is  spread 
out  over  the  plate.  In  the  vicinity  of  the  stagnation  point,  a  stagna¬ 
tion  flow  is  expected  and  its  properties  will  be  determined  by  the 
static  pressure  distribution  on  the  plate.  At  distances  from  stagnation 
point  sufficiently  large  to  approximate  a  zero  static  pressure  gradient 
at  the  plate,  the  flow  will  assume  the  form  of  a  wall  jet.  The  regions 
of  stagnation  flow  and  wall  jet  are  separated  by  a  transition  zone. 

An  analytical  program  was  conducted  to  study  the  heat  and  mass 
transfer  processes  in  the  stagnation  and  wall  jet  flow  regimes  of  an 
impinging  electric  wind.  Evaluation  of  the  heat  and  mass  transfer 
rates  at  the  stagnation  point  was  carried  out  by  approximating  the 
electric  wind  by  a  uniform  jet  of  velocity  VmflY.  the  maximum  velocity 
of  the  actual  wind,  and  half  width  xvj,  the  value  of  x  at  which  the 
velocity  of  the  wind  is  ^v^y .  Analysis  for  the  wall  jet  region  was 
based  upon  laminar  flow.  Comparison  between  theory  and  experiment  was 
accomplished  through  the  introduction  of  a  simplified  model  for  calcu¬ 
lating  the  total  heat  and  mass  transfer  to  the  test  plate.  According 
to  this  model,  the  stagnation  flow  region  extends  to  2x%  from  the  stag¬ 
nation  point  and  the  transport  coefficients  between  the  stagnation 
point  and  the  wall  jet  follow  a  linear  relation.  In  the  wall  jet  region, 
the  average  transport  coefficients  are  assumed  to  take  the  values  at 
a  point  two  inches  from  the  stagnation  point. 

Theoretically  the  transport  coefficients  for  wall  jets  depend  on 
two  parameters  R  and  l.  Three  different  methods  were  proposed  to  deter¬ 
mine  these  two  constants.  On  the  ground  of  physical  arguments,  the 
transport  coefficients  in  a  wall  jet  are  expected  to  increase  with  the 
maximum  wall  jet  velocity.  The  only  method  which  gives  results  in 
agreement  with  this  expectation  is  the  one  which  assumes  l  negligible 
and  R  determinable  from  the  maximum  wall  jet  velocity. 

The  total  mass  transfer  versus  current  relation  predicted  by  the 
theory  was  found  to  agree  satisfactorily  with  the  experimental  results. 
Validity  of  the  theory  was  further  verified  by  the  agreement  between 
the  absolute  magnitude  of  the  predicted  frost  deposit  and  the  measured 
data.  It  is  therefore  beleived  that  the  present  theoretical  model 
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should  also  be  adequate  to  predict  the  heat  transfer.  However,  the  com¬ 
parison  between  experimental  and  theoretical  heat  transfer  results  was 
found  to  be  less  satisfactory.  The  causes  of  this  discrepancy  demand 
further  investigation. 

The  theoretical  predictions  on  heat  and  mass  transfer  were  based 
on  the  assumption  of  negligible  norma]  velocity  at  the  frost  surface. 

The  assumption  appears  to  be  a  good  one  under  the  test  conditions.  With 
a  view  of  furthering  the  study  to  include  the  electric  wind  effects  upon 
other  mass  transfer  processes  such  as  drying  or  other  evaporation  pro¬ 
cesses  where  the  assumption  is  no  longer  valid,  an  analysis  was  made  of 
the  heat  and  mass  transfer  in  the  presence  of  non-negligible  normal 
velocity  at  the  surface  for  plane  stagnation  flows  as  well  as  plane 
laminar  wall  jet,  using  an  integral  method.  No  explicit  expressions  for 
heat  and  mass  transfer  coefficients  are  possible  for  this  case  and, 
therefore,  numerical  or  graphical,  solutions  must  be  employed.  No 
examples  are  given  in  this  connection’. 

A  simplified  model  of  heat  conduction  in  a  frost  layer  was  also 
analyzed  using  a  power  series  solution  introduced  by  Portnov.  As  expected, 
the  controlling  factor  responsible  for  the  unsteadiness  of  the  transport 
mechanisms  in  the  frost  layer  was  found  to  be  the  mass  flux  from  the 
surroundings  to  the  wall.  The  heat  flux,  on  the  other  hand,  was  found 
to  have  the  effect  of  counteracting  the  unsteadiness  caused  by  the  mass 
flux. 


8.2  KECOMMEDATIONS  FOR  FURTHER  RESEARCH 

The  impinging  flow  associated  with  an  electric  wind  is  of  particular 
interest  because  of  its  close  analogy  to  a  free  jet  impinging  on  a  sur¬ 
face.  In  recent  years  the  so-called  impingement  flow  has  been  increas¬ 
ingly  utilized  to  enhance  the  intensity  of  heat  and  mass  transfer  pro¬ 
cesses  in  various  industries.  At  the  present  time  no  satisfactory 
theory  on  the  impingement  flow  has  yet  been  developed.  A  study  of  the 
flow  and  the  heat  and  mass  transfer  processes  in  connection  with  the 
impingement  of  an  electric  wind  may  provide  useful  information  on  the 
impingement  flow. 

The  present  investigation  was  preliminary  and  incomplete.  More 
extensive  and  r^stematic  research  remains  to  be  done  to  furnish  a 
clearer  picture  and  to  explore  potential  applications  of  the  phenomena. 
The  following  are  some  suggestions  for  further  research: 

a.  Experimental 

Since  the  heat  and  mass  transfer  coefficients  depend  on  the  flow 
field,  the  investigation  of  the  flow  pattern  in  various  regions  of  the 
flow  field  is  of  primary  importance.  The  first  region  to  be  studied  is 
the  electric  wind  itself.  It  appears  that  the  flow  development  of  an 
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electric  wind  along  its  own  axis  has  never  been  investigated  before  and 
it  would  be  a  valuable  contribution  to  the  understanding  of  the  aero¬ 
dynamic  processes  of  an  electric  wind  if  such  measurements  can  be  accom¬ 
plished.  Two  cases  have  to  be  considered.  The  first  is  to  investigate 
the  variation  of  the  velocity  profile  with  distance  from  the  emitting 
electrode  in  the  absence  of  the  obstruction  by  plane  electrode.  Ifce 
other  case  is  to  examine  the  influence  of  the  solid  collecting  electrode 
on  the  development  of  the  wind.  It  is  also  desirable  to  determine 
whether  the  flow  is  turbulent  or  laminar  in  the  electric  wind. 

Systematic  measurement  of  the  velocity  profile  at  various  distances 
from  the  stagnation  point  along  the  plane  surface  is  indispensable  in 
understanding  the  transition  from  stagnation  flow  to  wall  jet.  It  helps 
locate  the  position  of  transition  and  is  also  essential  in  determining 
the  effects  on  the  stagnation  flow  of  the  ambient  still  fluid  and  of 
the  non-uniforaity  of  the  impinging  stream.  Measurement  of  pressure 
distribution  along  the  plane  electrode  is  helpful  in  providing  further 
insight  in  the  stagnation  flow  generated  by  a  non-uniform  impinging 
stream.  Such  measurement  has  been  made  by  Velkoff  for  positive  wire- 
plane  corona  and  by  Chattock  for  point-plane  corona,  both  positive  and 
negative. 

In  view  of  the  complex  nature  of  the  flow  field  under  consideration, 
measurement  of  local  heat  and  mass  transfer  coefficients  is  highly 
desirable  and  should  form  the  core  of  further  research  programs.  Other 
measurements  which  are  significant  and  should  be  attempted  in  future 
research  include  temperature  and  concentration  profiles  and  frost 
density  distribution  along  the  plate. 

b .  Theoretical 


Although  experimental  evidence  suggested  the  existence  of  a  simi¬ 
larity  between  the  current  density  distribution  and  the  velocity  pro¬ 
file  of  an  electric  wind,  it  should  be  emphasized  that  this  tentative 
similarity  needs  farther  analytical  verification.  Further  analytical 
research  should,  therefore,  be  directed  toward  establishment  of  the 
observed  similarity.  It  is  also  necessary  to  determine  the  role  played 
by  viscosity  and  entrainment  in  an  electric  wind.  Attention  should  be 
called  to  a  possible  analogy  between  an  electric  wind  an  thermal  plume 
generated  by  a  line  or  point  heat  source. 

Further  analytical  work  is  needed  in  determining  the  viscous  flow 
and  the  associated  heat  and  mass  transfer  processes  in  the  stagnation 
region  of  an  Impinging  non-uniform  stream.  According  to  the  stagnation 
flow  theory  of  a  uniform  stream,  the  boundary  layer  thickness  and  hence 
the  heat  and  mass  transfer  coefficients  are  constant  in  this  region. 
Experimental  evidences,  however,  indicated  that  the  local  heat  transfer 
coefficients  also  showed  a  bell-shaped  distribution  in  the  vicinity  of 
stagnation  point4.  Evidently,  this  discrepancy  is  due  to  the  non¬ 
uniformity  of  the  impinging  stream. 
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Two  other  problems  deserving  theoretical  considerations  are  the 
instability  of  a  laminar  wall  jet  and  the  transition  from  stagnation 
flow  to  wall  jet.  To  the  author's  knowledge,  the  first  problem  has  not 
even  been  mentioned  in  the  literature.  A  prompt  question  that  one 
intends  to  ask  is:  Is  the  instability  theory  of  a  free  jet  applicable 
to  a  wall  jet?  As  to  the  second  problem,  preliminary  prediction  of  the 
position  has  been  made  by  Brady  and  Ludwig10.  However,  the  problem  is 
still  far  from  being  comple te.ly  solved. 


APPENDIX  A 


PHOTOGRAPHS  OF  THE  GROWTH  OF  FROST  AND  FROST  THICKNESS  DISTRIBUTION 


Shown  in  Figs.  57  and  58  are  the  photographs  of  the  growth  of  the 
frost  strip  under  the  action  of  an  electric  wind.  Although  several 
corona  wires  were  mounted  on  the  electrode  support,  only  one  of  them 
(the  third  one  from  the  top)  was  active.  The  active  corona  wire  was 
located  along  and  at  a  distance  of  2  Inch  from  the  horizontal  center- 
line  of  the  test  plate.  The  pictures  were  taken  successively  at  4  min¬ 
ute  intervals  with  the  first  one  taken  one  minute  after  the  application 
of  the  field.  In  the  figures  the  letter  t  designates  the  time  when  the 
picture  was  taken  (the  field  being  applied  at  t  «*  0).  The  total  current 
for  this  test  run  was  about  150  pa  and  the  test  plate  temperature  was 
about  10 °F.  The  specific  humidity  of  the  ambient  air  was  not  recorded. 
The  vertical  ridges  mentioned  in  Section  1.2  can  also  been  seen  in  the 
pictures . 

Measurement  of  the  frost  thickness  distribution  along  the  test 
plate  was  made  with  a  depth  gauge  micrometer  mounted  on  a  sliding 
member,  the  latter  including  the  micrometer.  Measurements  were  made 
along  the  vertical  centerline  of  the  test  plate  at  \  inch  intervals  by 
positioning  the  sliding  member  along  the  channel.  Two  readings  were 
taken  at  each  location,  one  at  the  frost  surface  and  the  other  at  the 
test  plate  surface.  The  difference  of  the  two  readings  is  the  frost 
thickness  at  that  location.  To  facilitate  locating  the  position  of 
the  frost  surface,  a  pointed  cap  was  attached  to  the  end  of  the  measur¬ 
ing  rod  of  the  micrometer.  When  the  tip  of  the  modified  measuring  rod 
approached  the  frost  surface,  a  flash  light  wi »  turned  on  to  illuminate 
the  frost  surface  near  the  rod  tip.  The  rod  was  advanced  until  it 
reached  the  test  plate  surface  and  then  a  second  reading  was  taken.  In 
order  to  make  the  measurements  the  field  was  turned  off  after  the  frost 
had  been  allowed  to  grow  for  a  predetermined  period  of  time,  the  dis¬ 
charging  electrode  removed,  and  then  the  measuring  device  put  in  posi¬ 
tion. 


This  method  of  measuring  the  frost  thickness  is  not  precise,  for 
it  is  very  difficult  to  determine  the  exact  position  of  the  frost  sur¬ 
face.  Furthermore,  it  usually  took  about  ten  minutes  to  complete  the 
entire  readings  for  each  test  run.  To  prevent  the  frost  from  melting 
it  was  necessary  to  keep  circulating  the  refrigerant  through  the  tubing 
of  the  test  plate.  Hence  in  the  course  of  measurement,  the  frost  thick¬ 
ness  in  the  portion  yet  to  be  measured  might  have  changed  during  this 
period  of  time.  In  addition,  the  variation  in  time  of  the  frost  thick¬ 
ness  distribution  can  not  be  measured  with  the  present  device.  Because 
of  these  disadvantages  no  extensive  measurements  were  attempted,  nor 
was  it  intended  to  provide  quantitative  information.  Rather  it  was  only 
intended  to  give  seme  idea  of  the  pattern  of  frost  thickness  distribu¬ 
tion  at  an  arbitrarily  selected  instant  under  arbitrary  test  conditions. 


161 


■ 

Given  in  Fig.  59  are  frost  thickness  distributions  for  three  test  runs, 
each  being  measured  approximately  ten  minutes  after  the  field  was 
applied.  The  total  current  in  each  case  was  about  120  ua  and  the  test 
plate  temperature  was  about  10 °F.  No  specific  humidity  of  ambient  air 
was  recorded.  Positive  x  in  the  figure  corresponds  to  the  lower  half 
and  negative  x  the  upper  half  of  the  test  plate.  Again  a  single  hori¬ 
zontal  wire  located  at  the  centerline  of  the  test  plate  was  used  as 
discharging  electrode. 
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Frost  thickness  (in)  Frost  thickness  (in)  Frost  thickness  (in) 


APPENDIX  B 

CALCULATION  OF  WATER  VAPOR  CLUSTER  SIZE 


The  calculation  will  be  carried  out  for  the  case  where  the  current 
is  150  ua  and  the  total  frost  deposit  on  the  test  plate  is  1.5  gm  for 
five  minutes . 

Consider  each  ion  being  singly  charged  and  each  water  vapor  cluster 
containing  only  one  ion.  The  total  number  of  ions  reaching  the  test 
plate  in  5  minutes  for  a  current  of  150  pa  is 

%  ,  XOxjT*  xJOQ  .  2>3i  x  1017 
1.6  x  icr19 

where  1.6  x  1CH19  coulomb  is  the  electronic  charge.  The  total  number 
of  water  molecules  condensed  as  frost  on  the  test  plate  corresponding 
to  a  total  frost  deposit  of  1.5  ©a  is 

Bv  =  x  6.023  x  1023  =  5.02  x  1022 
18 

where  6.023  x  1023  is  Avogadro's  number.  If  each  water  vapor  cluster 
has  only  one  ion,  then  each  cluster  must  consist  of  the  following 
number  of  water  molecules. 

^  =  1.78  x  10s 

Ni 
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APPENDIX  C 


THE  SIX  COEFFICIENTS  OF  Y,  <pi  AND  92 

CP10  =  920  =  Tw  ,  Yi  =  0 

911  -  921  =  0  ,  911  +  921  “  |  (C3-C4TW) 

[Cl  -  (9lO+92o)J 

912  =  922  =0  ,  Y3  =  0 

913  =  923  =  0 

2 

913  +  923  = - C4Y2(9ll-»92l) 

k 

Y4  =  -  ^  Y2(9n+92i) 

914  ®  924  =0  ,  Ys  “  0 

925  “  925  "  0 

915  +  925  =  -  J  [l8Y2(9l3-Kfe3)  +  12Y4(cpii+<te1)] 

+  ~  (Y2)2(9l3+923)  I 

Y<S  =  -  ^l8Y2( 913^92 3)  +  12Y4(9ii+92i)j 

In  tejms  of  the  given  quantities  Tw,  T^,  wi  .  wif,  L_,  h,  hjn, 
P,  Pf  and  b,  the  following  expressions  are  obtained. 

Ya  =  £  ^  (wi.-v^) 

Y4  "  "  32  (p?  S  )2  +  *<*-*>]  (wioo-wiw) 

Ye  "  (^)3|  I 

•[  pLfhm(wiw-wiw)  +  hCT^-T*)] 

+  l  b  [  PLfV^co^iw)  +  h<V^>]2  |  (wico'Wlw) 


910  =  920  =  T„ 

911  =  921  =  i  ^^^(wi^-wiw)  +  h(Tw-Tw)J 

913  =  923  =  -  \  ^  ^  (p^h^h)  (wioo-w^) 

*  [pLfhm^ieo-^w)  + 

915  =  *»  ■  <-&&&■  |(^)a  (v^-wO^Vh) 

•  [  PLfMwi„-»i„)  ♦  h(T.-T„)]  [§  +  ±  (^“  ) 

(wi«>“Wiw)J  +  ^(^p£02  [pLfhm(wico“wiw) 

+  h(Tw-Tw)  ]  2  (vloo-wiw)  J 
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13  ABSTRACT 


An  experimental  and  analytical  program  was  undertaken  to  investigate  the 
behavior  of  frost  formation  under  an  impinging  electric  winf  induced  by  a  positive  wir« 
plane  corona  in  a  nonuniform  electric  field.  A  vertical  test  plate  precooled  to  sub¬ 
freezing  temperatures  was  exposed  to  the  ambient  air  and  the  effects  of  electric  fields 
were  studied  through  the  application  of  a  high  voltage  between  the  grounded  test  plate 
and  a  corona  wire.  Considerable  increases  in  heat  and  mess  transfer  were  found.  The 
changes  were  attributed  to  the  alteration  of  flow  field  from  free  convection  to  forced 
convection. 

\A  simple  theoretical  model  was  developed  for  predicting  the  effects  of  electric 
fields  upon  the  heat  and  mass  transfer  processes  associated  with  frost  formation.  Two 
regions  of  the  flow  field  created  by  the  impinging  electric  wind  are  of  primary  import¬ 
ance  in  heat  and  mass  transfer,  namely,  the  stagnation  flow  and  the  wall  jet  regimes. 
In  the  stagnation  flow,  a  linear  relation  for  the  local  heat  and  mass  transfer  coeffi¬ 
cients  was  assumed  between  the  stagnation  point  and  the  wall  jet  region.  In  the  wall 
Jet,  the  analysis  was  based  on  laminar  flow.  The  mass  transfer  predicted  by  this  mode] 
agrees  quite  satisfactorily  with  experimental  results,  but  discrepancy  between  theory 
and  experiment  was  observed  in  heat  transfer. 

*An  analysis  was  also  made  of  a  simplified  model  of  transient  heat  conduction  in  a 
frost  layer.  Using  a  method  due  to  Portnov  the  position  and  the  temperature  of  the 
progressing  frost-air  interface  was  obtained  in  a  series  ex}  ision  in  powers  of  \t. 

The  coefficients  up  to  the  power  n  =  6  were  given. 
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